SIC  FILE  COPY  ADA038215 


HEAVE  AND  ROLL  RESPONSE  OF  FREE  FLOATING  BODIES 
OF  CYLINDRICAL  SHAPE 


H.  O.  Berteaux 
R.  A.  Goldsmith 
W.  E.  Schott,  III 


February  1977 


TECHNICAL  REPORT 


Prepared  for  the  Office  of  Naval  Research 
under  Contract  N00014-75-C-1064 ; NR  294-044 
and  from  the  NOAA  Data  Buoy  Office. 

Approved  for  public  release;  distribution 
unlimited. 


UNCLASSIFIED  2/77 


*tCuR»TV  Cl.  *SS»  P»C  *Tio»4  O*"  This  PACE  (Whmri  Dat  m Frttarad) 


REPORT  DOCUMENTATION  PAGE 


I GOVT  ACCESSION  NO 


J^EAVE  AND  ^DLL  RESPONSE  OF  _FREE  FLOATING  BODIES 
OF  .CYLINDRICAL  SHAPE,  ~ ~ ~ 

t—  y 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


J recipient's  catalog  Humm 


s -REPORT  « PIBI&0  COVENed 


^technical 


ljJcfoStIl 


Ngfe'bl4-75-C-\^J^ 


.-nw/.l  - - 

JJ.  O.  ^erteaux,  Jl.  //oldsmi th/and 
JJ.  J-.^chott,  III  ^ i 


/:  £'-■ 


* PCnronuiNC  ORGANIZATION  namC  AnO  ADDRESS 


Woods  Hole  Oceanographic  Institution 
Woods  Hole,  MA  02543 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


NR  294-044 


I’  C ON  TROLLING  OEFtCE  NAME  AMO  ADDRESS 


Naval  Ocean  Research  and  Development  Activity 
Bay  St.  Louis,  Mississippi  39520 
ATTN : NORDA  400 


*ys report  nrrr — . 

(J/j/Feb aamw-  911  j 

f/C  n r Diml  -J 


112 


Ta  monitoring  agency  NAME  • AOORESSCir  Olltoront  from  Control  I Inf  OIUco) 


IS  SECURITY  CLASS,  fof  IM,  roporfj 

Unclassified 


ls«  OECL  ASSIFICATION  ;dornoraoing 
SCHEDULE 


■ • Distribution  STATEMENT  fof  thlm  Roporl) 

Approved  for  public  release;  distribution  unlimited. 


p o c- 

-V- 


IT  OiSTRiBuTioh  statement  rot  t>o  obotroct  ontmrmd  In  Bloc*  yt,  " Milormi  from  Roport)  1 i 

, ; » I2  1911  i 


\ V\  ' 

\ 

1 \ U 1-iLU  - 1 “ 


_ J t 

v.,i 


7V 


*8  SUPPLCMCnTARY  notes 


19  KEY  WOBOS  'Coorinu*  on  ravmraa  aid*  It  nacaaaary  and  Identity  by  block  numb**) 


1 . Spar  buoys 

2 . Buoy  Dynamics 

\.  Heave  and  Roll  of  Spar  Buoys 


20  ABSTRACT  'Conttnua  on  ravara*  mid*  It  nacaaaary  and  Identify  by  block  numb**) 

The  following  report  describes  a computer  solution  to  help  predict  the 
heave  and  roll  response  of  free  floating  bodies  of  cylindrical  shape  when 
excited  by  random  seas  with  known  spectra. 

The  basic  concepts  of  harmonic  analysis  and  statistics  used  in  the 
method  are  first  briefly  reviewed.  The  report  then  presents  a detailed 
derivation  of  the  linear  heave  and  roll  response  amplitude  operators  • ^Cont. ) 


DD  1 jan*7J  1473  EDITION  OF  1 NOV  ES  IS  OBSOLETE 

S/N  0 I02-01*-  MO' 


UNCLASSIFIED  2/77 


OCC 


• ECURITY  CLARIFICATION  OP  THIS  RAGE  fRtiw.  Dto  Bntmod) 


SECURITY  CLASSIFICATION  OF  THIS  PAOEf*?i«n  Dal* 


WHOI-77-12 


HEAVE  AND  ROLL  RESPONSE  OF  FREE  FLOATING  BODIES 
OF  CYLINDRICAL  SHAPE 


by 


H.  O.  Berteaux 
R.  A.  Goldsmith 
W.  E.  Schott,  III 


WOODS  HOLE  OCEANOGRAPHIC  INSTITUTION 
Woods  Hole,  Massachusetts  02543 

1 


February  1977 


TECHNICAL  REPORT 


Prepared  for  the  Office  of  Naval  Research  under 
Contract  N00014-7S-C-1064  ? NR  294-004  and  from  the 
NOAA  Data  Buoy  Office. 

Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government. 

In  citing  this  manuscript  in  a bibliography , 
the  reference  should  be  followed  by  the  phrase: 
UNPUBLISHED  MANUSCRIPT. 


Approved  for  public  release;  distribution 
unlimited. 


Approved  for  Distribution 


Melvin  A.  Rosenfeld,  Acting  Chairman 


Department  of  Ocean  Engineering 


ABSTRACT 


The  following  report  describes  a computer  solution  to  help  predict 
the  heave  and  roll  response  of  free  floating  bodies  of  cylindrical  shape 
when  excited  by  random  seas  with  known  spectra. 

The  basic  concepts  of  harmonic  analysis  and  statistics  used  in  the 
method  are  first  briefly  reviewed.  The  report  then  presents  a detailed 
derivation  of  the  linear  heave  and  roll  response  amplitude  operators, 
that  is  the  expressions  of  the  vertical  and  angular  displacements  produced 
by  a simple  harmonic  wave  of  one  foot  amplitude. 

The  second  part  of  the  report  reviews  the  computation  procedure 
and  the  program's  logic.  It  gives  a detailed  set  of  instructions  for  the 
program  users,  reviews  the  program's  capabilities  and  limitations,  and 
presents  three  case  studies. 

The  heave  and  roll  response  programs  are  written  for  use  with 
XEROX  SIGMA  7 computers.  Program  listings  are  given  in  the  appendix. 


K .'iu  Ufim  pf' 
1st!  S-iilts*  □ 


1 

i 

i 


Acknowledgments 


The  authors  wish  to  express  their  gratitude  to  Mr.  R.  Walden 
for  his  critical  review  of  the  manuscript. 

The  assistance  of  Mrs.  A.  Henry  and  of  the  Graphic.  Arts 
Department  in  typing  and  preparing  this  report  is  also  gratefully 
acknowledged. 

The  work  reported  herein  received  support  from  the  Office  of 
Naval  Research  (Contract  No.  N00014-75-C-1064;  NR294-044)  and 
from  the  NOAA  Data  Buoy  Office. 


TABLE  OF  CONTENTS 


Page  Number 

1.0  PROBLEM  STATEMENT 1 

2.0  THEORETICAL  BACKGROUND 

2.  1 Statistical  Response  of  Floating  Bodies  to 

Ocean  Wave  Excitation  1 

2.2  Derivation  of  Heave  and  Roll  Response 

Amplitude  Operators  (RAO).. 5 

2.2.1  Heave  response  5 

2.  2.  1.  1 Initial  conditions 5 

2.  2.  1. 2 General  equation  of  heave 

motion 7 

2.  2.  1.  3 Expression  of  the  forces 

applied  to  the  buoy 7 

2.  2.  1.  4 Expression  of  the  differential 

equation  of  heave  . . 14 

2.  2.  1.  5 Expression  of  the  heave  RAO 14 

2.  2.  1.6  Phase  relationship  between 

heave  and  wave 16 

2.2.2  Roll  response 17 

2.  2.  2.  1 Initial  conditions  17 

2.  2.  2.  2 General  equation  of  roll  motion  ...  18 

2.  2,  2.  3 Expression  of  the  moments 


applied  to  the  buoy 


18 


Page  Number 


2.  2.  2.  4 Expression  of  the  added 

moment  of  inertia  1^ 26 

2.  2.  2.  5 Expression  of  the  differential 

equation  of  roll 27 

2.  2.  2.  6 Expression  of  the  roll  response 

amplitude  operator 28 

2.  2. 2.  7 Phase  relationship  between 

roll  and  wave. 29 

3.0  COMPUTER  PROGRAMS 

3.  1 Heave  Computer  Program  (HERAO) 29 

3.1.1  Program  logic 29 

3.1.2  Program  input 31 

3.  1.  3 Program  output 40 

3.2  Roll  Computer  Program  (ROLLRAO)  41 

3.2.  1 Program  logic 41 

3.2.2  Program  input 42 

3.  2.  3 Program  output 52 

4.0  CASE  STUDIES 

4.  1 Heave  and  Roll  Response  of  r.  Small  Flat  Cylinder  ...  53 

4.  1.  1 Program  input 54 

4.1.2  Program  output  56 

4.  2 Heave  and  Roll  Response  of  a Ballasted 

"Telephone  Pole"..... 57 


Page  Number 


4.  2.  1 Program  input 57 

4.2.2  Program  output 61 

4.  3 Heave  and  Roll  Response  of  a Complex 

Shape  W.H.O.I.  Soar  Buoy 61 

4.  3.  1 Program  input  63 

4.  3.  2 Program  output  ...........  66 

5.0  CONCLUSIONS  AND  LIMITATIONS 67 

6.0  REFERENCES 69 

7.0  APPENDICES  70 

I.  Expression  of  Linearized  Damping  Coefficient....  70 


II.  Evaluation  of  the  Coefficient  "B"  of  Damping 

Moment 72 

III.  Evaluation  of  the  Coefficient  "D"  of  Wave 

Drag  Moment  75 

IV.  Evaluation  of  the  Coefficient  "P"  of  Wave 

Inertia  Moment..... 77 

V.  Evaluation  of  the  Coefficient  "Ip."  of 


Added  Moment  of  Inertia  B0 

VI,  Computation  Method  for  Coefficients  "B", 

"D",  "P" 81 

VII.  Heave  Program  Listing 83 

VIII.  Roll  Program  Lir  ting 91 


1 


.0  PROBLEM  STATEMENT 

The  heave  and  roll  motion  of  a cylindrical  body  of  constant  cross 
section  when  excited  by  a simple  harmonic  wave  is  a relatively  straight- 
forward problem.  However,  very  few  buoys  can  be  realistically  modeled 
as  a pillbox  or  a telephone  pole.  Most  spar  buoys  are  made  of  circular 
cylinders  of  varying  diameters  (see  Fig.  #1).  Some  spar  buoys  extend 
to  considerable  depths  below  the  water  level.  Furthermore,  most  sea- 
ways are  not  made  of  regular  harmonic  waves  of  single  frequency  and 
amplitude  and  in  general  irregularity  and  randomness  of  the  sea  surface 
will  prevail. 

The  objective  of  this  report  is  to  present  a method  which  can  be 
used  to  compute  reasonable  expectations  of  vertical  and  angular  displace- 
ment that  a complex  shape  buoy  will  experience  when  free  floating  in  a 
random  stationary  seaway. 

The  computer  solution  presented  in  this  report  was  originally 
derived  to  investigate  the  dynamic  behavior  of  specific  spar  buoys  used 
by  the  Woods  Hole  Oceanographic  Institution.  This  solution  is  here 
presented  in  a generalized  form,  with  the  hope  that  it  becomes  a con- 
structive addition  to  the  solutions  already  in  the  literature. 

0 THEORETICAL  BACKGROUND 

2.1  Statistical  Response  of  Floating  Bodies  to  Ocean  Waves  Excitation 
Readers  unfamiliar  with  the  probabilistic  theory  of  ship 
and  buoy  dynamics  should  resort  to  References  I,  II,  and  III 
for  a theoretical  introduction  to  the  subject. 


Fig.  No 
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Basic  concepts  borrowed  from  this  theory  and  used  in  the 
formulation  of  the  heave  and  roll  computer  programs  described  in 
this  report  are  hereafter  summarized. 

If  the  probability  density  function^6^VJ  of  the  wave  ampli- 
tudes for  a given  seaway  can  be  explicitly  expressed,  then 

the  expectation  of  certain  values  of  wave  amplitudes  can  be 
directly  computed. 

For  example: 

- The  most  probable  amplitude  is  the  value  of  ^ 
for  which  , 

£ M * 0 


The  average  amplitude  0(  is  given  by 


(2. 1.  1) 


ty  = Jxpf'x) 


(2.  1.2) 


w 

- The  average  of  a fraction  f of 

wave  amplitudes  larger  than  a given  amplitude  ^ a 


can  be  obtained 
oo 


fry  c 1 J (2.1.3) 

fpC*)  ^ OCo 


fa' 


etc. 


When  certain  restrictive  conditions  prevail,  wave  amplitudes 
have  been  found  to  follow  a Rayleigh  distribution  given  by: 


bf¥)  = Ajt 

__  %>■ 

where  is  the  mean  square  value  of  the  wave 

amplitudes. 


(2.  1.4) 
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This  probability  density  function  has  been  used  to  compute 
the  expectation  of  particular  wave  amplitude  means  and  maxima. 
Results  of  these  computations  are  found  to  be  proportional  to  the 


root  mean  square^#  of  the  wave  amplitudes.  Those  retained 
in  this  study  are  summarized  in  Tables  land  2,  "Value  of  expected 
means"  and  "Value  of  expected  maxima.  " 

Now  let\^  the  spectral  density  function  of  the  wave 

amplitudes  of  the  given  sea  way.  Assuming  the  seaway  to  be 
stationary,  then  the  integral  over  all  positive  frequency  ranges 
of  is  equal  to  the  mean  square  value  of  the  wave  amplitudes, 

i.  e. 

oo 

'^r  = (2.1.5) 

o 

This  result  can  be  used  to  compute  X2  and  . The 

value  of  & thus  obtained  can  in  turn  be  used  to  compute  the  ex- 
pected wave  amplitude  means  and  maxima  listed  in  Tables  1 and  2. 


From  the  definition  of  the  mean  square  value 

j-  rV 

/V-*=o  A/  '—I  L 

and  the  result  (2.  1.  5)  it  is  clear  that  the  quantity 


o 

is  proportional  to  the  amplitude  of  the  elementary  component 


(2.  1.  6) 


wave  of  the  spectrum  with  frequency  6*2^  . 

Now  if  is  the  expression  of  the  linear  response  of  a 
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Table  1 


Wave  Amplitude  Means 


Fraction,  j' , of 
Largest  Ampli- 
tudes Considered 

Mean  Values 
ij  -r  J7l 

mam 

2.  359 

1.800 

KB9  • 

1.416 

0.  50 

1.256 

1.00 

0.  886 

Table  2 


Expected  Maximum  Amplitudes 


Number 

of 

Waves 

50 

100 

500 

1, 

000 

limn 


Maximum  Wave 
Amplitudes 

max  * ▼ » 


2.  12 
2.  28 
2.  61 

2.  78 

3.  13 
3.47 
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free  floating  body  to  a simple  harmonic  wave  of  unit  amplitude 
and  frequency  Od  , then  clearly  the  quantity 


dti)  O 


ywiAto-; 


is  proportional  to  the  amplitude  of  the  body  response  to  the 
elementary  component  wave  of  frequency  . 


It  thus  follows  that 


Y(CJ)  'ffco* ) cl(X> 


is  r jportional  to  the  amount  of  the  response  mean  square 
value  contained  in  the  frequency  band  oUd  centered  at  . 

The  response  mean  square  value  A1  ^ is  therefore  given  by: 


The  response  /fij)  of  the  body  being  linear,  the  probability 
density  function  of  the  response  will  also  follow  the  probability 
density  function  of  the  wave  amplitudes.  Thus  the  results  tabu- 
lated in  Tables  1 and  2 can  be  used  again,  together  with  expression 
(2.  1.  9)  to  compute  statistical  means  and  maxima  of  body  response 
amplitude. 

For  example,  the  average  of  the  one  third  highest  response 
amplitudes  will  be  given  by  _____ 

7 1 = wei/fc 


— jjYru) 


(2.  1.  10) 
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Empirical  formulation  of  wave  amplitude  spectra  used  in 
the  computer  program  are: 


- 

JfuA*.  16.876  £ V*** 


ft  - sec  (Pierson  Moskowitz)  2.  1.  11 


where  V is  the  wind  speed  (knots). 


-/ oso 


Jru»=  Ststfe**** 

pu>s 


ft  - sec  (Bretschneider)  2.1.12 


- <T3  O 


\fra>)  = ^5^  e.G**» 
Ts*a>* 


ft  - sec  (I.S.S.C.) 


2. 1. 13 


In  formula  (2.  1.  12)  and  (2.  1.  13)  i8  the  significant  wave  height 
(feet)  and  is  the  significant  wave  period  (seconds). 

Z*  2 ££,?lvatlon  of  the  Heave  and  Roll  Response  Amplitude  Operators  (RAO) 
2.2.1  Heave  response 

2.2.  1.  1 Initial  conditions 


T 

A 

i 


Fig.  No.  2 
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w 


i 


i 


f 


Let  us  consider  a simple  harmonic,  deep  sea  wave,  as 
shown  on  Fig.  No.  2.  The  coordinates  of  a point  on  the  surface 
of  this  wave  are  given  by: 


J ta  At**  (<*>£  - 

g r ‘/Vo  (cot  — ^ S') 


(2.2.  1) 
(2.2.2) 


If  we  select  to  observe  this  wave  at  S * @ * then  the  parametric 

equation  of  the  water  particle  motion  around  this  point  become: 

a Cd  t 

£ s s4  C<m  cot 

where  4 i.  the  wave  amplitude,  and  Cd  the  wave  angular 
frequency. 

We  also  know  from  the  simple  harmonic  wave  theory  (Ref.  I, 
pp.  14-27)  that  the  parametric  equations  of  water  particles  at  any 

depth  7 below  the  mean  water  level  would  then  be 

- 

^ *c  A*'**  Ct)  t 

CvO  Cjt 

where  /(  is  the  wave  number,  q s J* { for  deep  water  waves, 

<f 

^ being  the  gravity  acceleration. 

The  vertical  components  of  water  particle  velocity  and  ac- 
celeration would  in  turn  be  given  by: 

£ = — Q CO  t 

2 e _ cot 

In  this  case,  at  time  t*  0 + £ , the  amplitude  of  the  vertical 
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f 


displacement  of  the  water  particles  start  to  decrease  from  their 
maximum  value,  the  water  particles  vertical  velocity  component 
starts  to  increase  and  is  in  the  downwaids  direction,  and  the 
water  particles  vertical  acceleration  component  starts  decrees- 
ing  and  is  also  in  the  downwards  direction. 

As  shown  on  Fig.  No.  3,  let  ft  be  the  distance  from  the 
still  water  surface  to  the  buoy  water  line.  At  time  (f  c 0+  £ the 
buoy  is  assumed  to  move  downwards,  that  is  the  distance  ft  is 
increasing. 

2.2. 1.2  General  equation  of  heave  motion 

The  equation  of  heave  motion  will  be  obtained  from: 


2.2.  1.  3 


\ •* 

= [4*\-¥AAa' Wy/K 

sum  of  the  vertical  forces  applied  to 


(2.2. 


the  buoy, 

- mass  of  the  buoy 

added  mass  of  the  buoy  due  to  the  water  entrained 
in  the  vertical  direction 

- virtual  mass  of  the  buoy  (in  the  vertical 


3) 


direction). 

Expression  of  the  forces  applied  to  the  buoy 

The  vertical  forces  applied  to  the  buoy  are: 

- Its  weight  " 

- The  resultant  " J 3,1  of  the  pressure  forces  exerted 

by  water  particles  on  the  top  and  the  bottom  plates  of 


Fig.  No.  3 


the  watertight  compartments  of  the  buoy 

- The  damping  force  "J]"  resulting  from  the  water 
opposing  the  buoy  vertical  motion 

- The  friction  force  "(S"  exerted  by  the  water  particles 
vertical  velocity  on  the  buoy 

- The  inertial  force  "Jf  " exerted  by  the  water  particles 
vertical  acceleration  on  the  buoy. 

Forces  in  the  direction  of  increasing  (downwards)  will 

be  considered  positive. 

The  expression  of  these  forces  is  obtained  as  follows: 

“ Pressure  force  "1  M. 

With  the  initial  conditions  assumed,  the  pressure  jo  at  a depth 

Z is  eiven  by; 

jo  = + <L  CoS  Lot  J 

To  help  find  a general  expression  for  the  resultant  ^ , let  us 
consider  the  spar  buoy  shown  in  Fig.  No.  3. 

At  the  bottom  of  the  buoyancy  tank 

Z = x+hz. 

where  depth  of  the  bottom  plate  below  buoy  water  line 

Assuming  , the  upwards  pressure  force  7^  on  the  tank 

bottom  is  thus  given  by 

where  = area  of  bottom  plate  subjected  to  water  pressure 
(the  entire  area  of  the  plate  in  this  case). 
Similarly  the  pressure  force  ~P  on  the  top  plate  of  the  buoyancy 
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tank  is  given  by 


~kJ), 


£ = ft  ft+l),  + 4e.  c*s  cot jST 

where  \Sy  ~ area  of  top  plate  subjected  to  water  pressure. 

If  the  spar  mast  has  a cross  section  and  is 
watertight  then  obviously  Sy  ® 

The  resultant  _/  will  be  the  difference  between  the  bottom  pres- 
sure force  and  the  top  pressure  force.  Being  in  the  upwards 

J1  . T>  /p  'P  \ . 

direction,  J,  *-//*-  iy  l.  e. 

ffos-SrJx  + it$T)-/4casco£(fjt.  _ Srzj  j 

The  constant  terms  in  the  expression  of  the  pressure  force  must 
equal  the  buoy  static  weight  . This  can  be  easily  established. 

Noting  that 

where  //  = length  of  the  buoyancy  tank,  the  constant  terms 

( saiL  - sr6,J 

can  be  written: 

JS  fa  (6^)-  Sr 6,7  = ((Sq  - Sr)/},  + SqMJ 

which  obviously  is  the  sum  of  the  weight  of  the  water  displaced  by 
the  immersed  portion  of  the  mast  and  by  the  buoyancy  tank  under 
equilibrium  conditions,  and  therefore  is  equal  to  the  buoy  weight. 

The  sum  of  the  weight  force  and  the  pressure  force  can  then  be  in 


general  expressed  by: 

"P, 


(2.2.4) 


where 


is  the  surface  at  a depth 


hi 


subjected  to  the 
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pressure,  sS^  ia  positive  if  the  pressure  exerted  upon  it 
is  in  the  upwards  direction,  and  vice  versa  is  nega- 

tive if  the  preasure  exerted  upon  it  is  in  the  downwards 
direction. 

This  expression  can  further  be  simplified  and  written: 

7>+  W=  - + (A/t  CoS  cot 


where 


c =JS  (se-s')^fSSe- 

is  the  heave  restoring  force  constant  and  ^ . 


(2.2.5) 
(2.2.  6) 


(2.2.7) 


<■  c 

- Damping  force  'iD". 

The  damping  force  exerted  by  the  water  on  a buoy  component 
" L"  will  be  assumed  to  be  directly  proportional  to  the  buoy 


speed  ^ . It  will  therefore  be  of  the  form 

J>C  ~ -~£>c  £ 

where  is  the  linearized  coefficient  of  damping 


(2.2.8) 


am 


associated  with  buoy  heave  motion. 

It  can  be  shown  (see  Appendix  I)  that  the  general  expression  of 

linearized  damping  coefficients  "dj1  for  periodic  motion  of 

plitude  X and  frequency  20  is  of  the  form 

d x Jl  (2.2.9) 

3K  O 

<3 

where^*  = water  mass  density  = 2 slugs/ft3 
Cp  = conventional  drag  coefficient 
vT  = area  normal  to  the  flow. 


r 

» 
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In  order  to  keep  the  differential  equation  of  heave  motion  linear, 
an  arbitrary  constant  value  of  average  heave  motion  must 

be  selected  to  compute  the  linearized  damping  coefficients  l?l.  . 
The  value  of  X^,  selected  is  left  as  an  input  for  the  program 
users.  One  can  use,  for  example,  a reasonable  fraction  of  the 
average  wave  amplitude  for  the  sea  state  considered  in  a given 
study. 

The  expression  of  then  becomes 

i>c  = Xj,  c ^4 


(2.2. 10) 


The  total  damping  force  is  thus  finally 

j>.  £4  . X, 

L l 

or  simply 

D*  -3<k 


with 


(2.2.11) 


i L 

- Wave  induced  drag  force 

The  drag  force  resulting  from  the  water  particle  impinging 

* 

with  a velocity  on  a buoy  component  " is  also  assumed  to 
be  linearly  proportional  to  £ . It  therefore  will  be  expressed  by 

» cirL  (2.2.  12> 

where  Q is  the  linearized  coefficient  of  drag  associated 
with  water  particle  velocity. 


Following  previous  reasoning  the  expression  of  CTc  will  be 
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given  by 

Ci  s J Q>i  ^ ^ (2.2.13) 

where  is  now  the  arbitrary  average  value  of  wave 
amplitudes  retained  for  the  particular  study. 

A comment  should  also  now  be  made  regarding  the  water  particle 

* . 
velocity  V?  . It  will  be  recalled  that  the  expression  of  h is 
U ~«Z.  C 

to  £. 

In  the  case  of  a plate  or  a cylinder  of  small  height  placed  at  a 
distance  fi  below  the  buoy  W<L.  and  if  /V  , then  2.-4. 

and  the  speed  of  the  water  particles  acting  on  this  plate  is  well 
established.  On  the  other  hand,  if  the  cylinder  is  one  of  con- 
siderable height,  as  for  example  the  buoyancy  tank  shown  in 
Fig.  No.  3,  then  the  speeds  at  the  top  and  at  the  bottom  must  be 
somehow  averaged  and  replaced  by  a unique  equivalent  speed. 


For  simplicity,  one  could  consider  this  averaged  speed  to  be  the 
speed  at  the  depth  of  the  cylinder  midpoint.  For  the  buoyancy 
tank  previously  mentioned  this  speed  would  then  be: 

_ /?ii/  e *■'  st>r«  ait 


More  appropriate  values  of  equivalent  depths  could  also  be  devised. 
With  these  remarks  in  mind,  the  expression  of  the  friction  force 
becomes 


Gr  9 jjTj  ss  - ^<0  U)  6 J>  (J;  t 


> c c 

being  the  true  or  the  equivalent  depth  of  the 

n,  9 

component  . 
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More  simply  written. 


£ - - 00  JU+1  Cd  t 


with 


-v-  «>£j|/6^eA  „ 

C • • 

c 

- Inertial  force  "L 

The  inertial  force  I-  produced  by  the  water  particle  accelera 
• « 

i 

tion  ^ on  a component  "c"  of  the  buoy  is  given  by 

■£  - -*ii  i 

/ t 

where  is  the  added  mass  of  the  component  and 


is  given  by 


with  = added  mass  coefficient  of  component 

= volume  of  the  component  (ft^). 

The  values  of  and  1/^  depending  of  course  on  the  dimen- 

» 

sions  and  shape  of  the  component  are  left  as  an  input  for 

the  program  users. 

The  remarks  on  the  averaged  value  of  the  water  particle  speed 
also  apply  for  the  water  particle  acceleration. 

The  expression  of  the  inertial  force  "J" " is  therefore  given  by 

/»  JT  = ’ CytCdt  . €.  (2.2. 


1 •=.  — Q/f  co  ^ OJ  t 

& m Z^l « e~^ 


- 14  - 


k** 


2.  2. 1.4  Expression  of  the  differential  equation  of  heave 

Using  <X  yields: 

c. 

_ Cfi  + MA cos<^>  £ - 3o(  - A/A ol>£ - <pA «>  c**su±£  * % 

or , 

Cix  + 5%  4 % =■  A j (£f-k)lcpja>su>t--  c/tosi+icot j 


This  expression  can  be  further  reduced  to: 

ii  — / i 

C*  4 B4u-*t¥  ft  a 4 CoS  (uit-t^Tj 

where  ^ , the  exciting  force  is  given  by 

E * V J(M-  ?«?)  *+ 

and  <T  , the  phase  angle  between  the  wave  and  the  force  is 
given  by 

<r-  a/6o 


(2.  2. 


(2.  2 


M-  u>lQ 

2.2. 1.5  Expression  of  the  heave  response  amplitude  operator 

Let  us  assume  that  a particular  solution  of  the  heave 
equation  is  given  by: 

/x  - Xo  c°s  + ty) 

where  is  the  phase  angle  between  the  exciting  force  and  the 
heave  response.  Then, 

•Jf  « Xo  ( 6l0  Coo  OJ  t ^ ) 

o CO  ^ ’ /Qjm  00  £ Ckxs  -f  CtrOuj  t \jj  J 

# * y 

a - (Xo  Co  L ^Coo  C 0<£  Corp  ojs  — lUtyt  (jj ^ 


(2.2. 


Introducing  these  values  of  * •#  --X  in  the  equation  of 


18) 


.19) 


20) 
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heave  motion  and  ignoring  for  the  moment  the  phase  angle  Q* 
yield: 

CX  ( tot  U)t/9*4to<pj 

_ 5 X°  ^ ( *£4  **  to  £ £<3t>  ^ y-  <^vT  o<j  £ 

- 'HlyX*  <4X  ( too  cot  t/s  - if***  tot  &'u  j/J  * 'FCoStot 


Thus, 


Xo  j ( C - Wy  e*ro  (js  _ \ — l 

/,  j(-  C i - Bex)  toSlfj  SS  0 


From  the  second  result. 


tl&As!  (f)  ~ 


Therefore 


S-4M 


1»* 


/ = 


- £c>  jg> 

C - 

- °o  & 

\j(C-/*¥  to  L)  \ (<o  3j 

C - **fv«jz' 


(2.2.21) 


Introducing  these  values  of£/*<  ^ and  C&,s  ^ in  the 


first  result,  yields: 


X,  = 


The  expression  of  the  heave  response  is  thus  finally  given  by 


- lb  - 


* = ^/06  + n.a,) 


(2.  2.22) 


The  response  amplitude  operator  being  the  ratio  of  the 
heave  amplitude  ^ by  the  wave  amplitude  is  thus  in  turn 
expressed  by 

-R.A.O.  = — = 

A 


As  previously  established,  the  phase  angle  0^  between 
wave  and  exciting  force  is  given  by 

_/  (0 


C^=  (qm~ 


, _ u>  * 

.*■  c<  e 7*- 


(2.2.24) 


oyZbktizfr' 

«-  c 

The  phase  angle  ft  between  the  exciting  force  and  the 
heave  response  is  in  turn  given  by 

f*  ~ ^ j’J'Sc  - 

The  phase  angle  <^*  between  wave  and  heave  response 

ft 


(2.2.25) 


is  finally  given  by 


(2.2.26) 


T 
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2.2.2  Roll  Response 
2.  2.  2.1  Initial  conditions 

Let  us  consider  again  the  simple  harmonic,  deep  sea  wave 
shown  in  Fig.  No.  2. 

The  slope  of  this  wave  is  given  by: 

a A-***  [ U)t  - ) 

If  we  again  select  to  observe  this  wave  at  5 s 0 , the  ex- 
pression of  the  slope  becomes 

s Cofc 


The  horizontal  components  of  the  velocity  and  acceleration 
of  a water  particle  at  a depth  Z will  in  turn  be  given  by: 

£ X 5 CO  /I &-  GoC  cot 

c * 

3 = ~GQ  yi  e u)£. 


At  the  time  0 + £ the  magnitude  of  the  horizontal 

particle  velocity  is  maximum  and  is  positive  (i.e.,  in  the  direc 
tion  of  increasing^  ),  the  magnitude  of  the  horizontal  particle 
acceleration  is  minimum  and  in  the  opposite  direction,  and  the 
magnitude  of  the  slope  is  minimum  and  starts  to  increase. 

Let  © be  the  angle  of  roll,  measured  from  the  vertical 
in  a clockwise  direction.  At  time  $+£,,  the  buoy  will  be 
assumed  to  roll  in  this  direction,  i.e.,  the  angle  of  roll  is 
increasing. 


These  initial  conditions  are  depicted  in  Fig.  No.  4.  Ro- 
tation of  the  buoy  is  assumed  to  take  place  around  the  buoy  center 
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of  gravity. 

2.  2.  2.  2 General  equation  of  roll  motion 

The  equation  of  roll  motion  is  given  by: 


(2.2.27) 


where  2*i  - sum  of  the  moments  applied  to  the  buoy 

4 

c r 

/ = moment  of  inertial  with  respect  to  c.g.  of  buoy 

T - added  moment  of  inertia  due  to  entrained  water, 
also  with  respect  to  buoy  c.g. 

- virtual  moment  of  inertia  = £ + Tp: 

2.  2.  2,  3 Expression  of  the  moments  applied  to  the  buoy 
Moments  applied  to  the  buoy  are: 

- righting  moment  caused  by  displacement  of  center  of 
buoyancy,  45* 

- damping  moment  due  to  buoy  motion  in  the  water 

- friction  moment  due  to  drag  forces  induced  on  the  buoy 

by  horizontal  water  particle  velocity  . 

- inertia  moment  due  to  inertia  forces  induced  on  the 
buoy  by  horizontal  water  particle  acceleration 

Clockwise  capsizing  moments  will  be  considered  positive, 
and  vice  versa  counterclockwise  righting  moments  will  be  con- 
sidered negative. 

The  expression  of  these  moments  can  be  derived  as 
follows: 
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- Righting  moment  - % 

The  righting  moment  opposes  buoy  motion.  Its  value  is 
given  by: 

^ -4j  - _ A/p (6-  (/fstmcutj  (2. 

where  (^/  - buoy  weight 

r - distance  from  buoy  center  of  gravity  to 
buoy  metacenter. 

- Damping  moment  - 2Z» 


The  drag  forces  due  to  buoy  motion  alone  oppose  the  roll  both 
above  and  below  the  buoy  center  of  gravity.  Therefore,  the 
damping  moment  is  negative.  Its  expression  is  derived  as 
follows: 


VV  V 


2.  28) 


Fig.  No.  5 
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Let  ud  consider  an  elementary  buoy  section  at  a distance  A*  from 
the  buoy  c.g.  (See  Fig.  No.  5) 

The  elementary  damping  force  on  this  element  will  be  assumed 
to  be  of  the  form:  • 

d£m  ifrjr'&cn  Q 

or,  for  small  angles  of  roll  G 

where  b a linearized  damping  coefficient  again  given  by: 


j ? = fluid  density  = 2 slugs/ft^ 

O = drag  coefficient  for  cylinders,  normal  flow 

s^’/4/=  area  across  the  flow  = 4( f*J  d t*  wi  th  the 

/ 

cylinder  diameter  at  distance  P* 

K/rJ  = amplitude  of  cyclic  motion  at  distance  r* 
w - r9  (in  order  to  keep  the  equation  of  motion  linear 
an  arbitrary  constant  value  of  0 must  be  selected, 
say  (B m G ), 

60  = angular  frequency  of  cyclic  motion,  which  under  steady 
state  conditions  should  equal  the  frequency  of  the  ex- 
citing wave. 

The  expression  of  the  damping  force  thus  becomes:  # 

JFb  = /£  j>  O & * Jfrj  t*  dr  j r 0 
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or 


c/%  = c/ti)  d(r)  r Vr  Q 


where 


(2.2.29) 


The  moment  of  this  elementary  force  is: 

-~r  /5  = - ^ r^c/r 

and  the  total  damping  moment  is  found  from 


or 


where 


( fro  qio 


(2.2.  30) 


(2.  2.  31 : 


Appendix  II  outlines  a method  for  computing  these  integrals. 

- Wave  drag  moment  - 

Drag  forces  due  to  water  particle  velocity  will  tend  to  capsize 
the  buoy  or  to  upright  it  depending  on  their  point  of  application 
with  respect  to  the  c.g. 

The  resulting  moment  will  thus  be  positive  above  the  c.g.  and 
negative  below  the  c.g. 
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Fig.  No.  6 


Consider  again  an  elementary  buoy  section  of  area  J/r/  Jr 
at  a distance  from  the  buoy  center  of  gravity  (Fig.  No.  6). 
The  elementary  drag  force  due  to  the  water  horizontal  velocity 
on  this  elementary  section  will  be  assumed  to  be  of  the  form: 

JF_  = / 

where  the  linearized  damping  coefficient,  will  be  expressed 

by: 

c(r/  = J-f  Q X/rl  u) 

3 ^ 4 

'H  in  this  case  is  the  amplitude  of  the  water  particle  cyclic 
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i 


motion  and  is  therefore  given  by 

){(rj  = y~>  e"  ^ 

Here  again,  in  order  to  maintain  linearity  in  the  expression  of 
the  roll  RAO,  an  arbitrary  constant  amplitude  must  be 

selected.  One  could,  for  example,  select  the  average  amplitude 
4 of  the  waves  in  the  particular  sea  state. 

With  these  remarks  in  mind,  the  expression  of  can  be 


written: 


or 


~ -4* 


Cfr'j  - _v_  f Cb  4 (r)  G- 


, - 
c(r)  — h)  <4 (r)  e.  d'r 

The  expression  of  the  elementary  drag  force  4F  thu! 


where 


(2.2.  32) 


becomes 


S /S/4 CcrOcot  dr 


The  moment  of  this  elementary  force  is  in  turn: 


» r <dF  » co  Ctfo  cot  F 6 c Jr 


Noting  that  the  drag  forces  have  a tendency  to  capsize  the  buoy 
when  applied  above  the  buoy  c.g,,  and  to  upright  the  buoy  when 
applied  below  the  c.g.  , the  expression  of  the  wave  drag  moment 
becomes: 
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or, 

r J)^iO  CtOOji  (2.2.33) 


where 


/ /'■r  -ife  #'fe  -i*i  > 

dfy  i (2.2.34) 

A method  for  the  evaluation  of  the  coefficient  is  outlined 


in  Appendix  III. 

- Wave  inertia  moment  - 


The  elementary  inertia  force  dj.  due  to  the  water  particle  hori- 


zontal acceleration  acting  on  an  elementary  buoy  section  of 

£ 

Jl/  s A d(r*)  dr 
4 


volume 


is  of  the  form: 


dl  = C„j  JV  f 


Cro 


0 


or,  for  small  angles  ® ^ .. 

JL  = dfnjJr 

where  = coefficient  of  added  mass  for  cylinders.  This 

elementary  force  can  be  more  simply  written: 

#• 


Jl  * f I d(rf  Jr 
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(2.2.35) 


Fig.  No.  7 


The  moment  with  respect  to  the  c.g.  of  this  elementary  force  it 


in  turn  given  by 


r'w r, IJr  z. 
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Nothing  again  (Fig.  No.  7)  that  inertia  forces  due  to  the  wave 
action  have  a tendency  to  upright  the  buoy  above  the  c.g.  and  to 
capsize  it  below  the  c.g.  , the  expression  of  the  total  moment 
will  be  given  by: 


or, 

cuk 

*~i'  (2.2.37) 


A method  for  the  evaluation  of  the  coefficient  ^ is  outlined 
in  Appendix  IV. 

2.  2.  2.  4 Expression  of  the  added  moment  of  inertia  Ip 

The  added  mass  of  an  elementary  buoy  section  of  volume 
Jv,  located  at  a distance  t*  (rom  the  buoy  c.g.  is  given  by 


<d*t 1 = C*  p c/V 


r 

*1  = added  mass  coefficient  = 1 for  cylinders. 

The  moment  of  inertia  of  this  elementary  mass  with  respect 
to  the  c. g.  is:  » . f 


d[p  = A-V*,  = tydV 


and  the  total  moment  of  inertia  is 


4 -M  -fffvVv 


(2.2.  38) 
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Thus 


71.  - moment  of  inertia  of  the  water  displaced  by  the  buoy 
r 

with  respect  to  the  buoy  c.g. 

J^p  can  be  evaluated  following  the  method  outlined  in  Appendix  V. 

2.  2. 2.  5 Expression  of  the  differential  equation  of  roll 

Summing  the  moments  and  applying  the  angular  form  of 
Newton's  law  yield: 

* • 


The  resulting  equation  of  motion  is  then: 

lyO-+BO-t  CO  = ^ 'j,  cot  +J)  Cl)  J 

where  C - is  the  roll  restoring  constant. 

The  equation  of  roll  motion  can  be  further  reduced  to: 

C0+  0 ~ W &S(u£+r)  (2.2. 39) 

where  A , the  exciting  torque  due  to  wave  action,is 
given  by: 


M * J j/p^  P*')  *■■+  (pu>)  *■ 

and  the  phase  angle  between  wave  and  resulting  torque  is 
in  turn  given  by: 

, Q?  + Tfc*- 
- to*'  * 


(2.  2.40) 


(2.2.41) 


J>co 
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2.2.2.  6 Expression  of  the  roll  response  amplitude  operator 

Assuming  again  that  a particular  solution  of  the  differential 
equation  of  roll  is  given  by; 

9 s 0e  Cao  ('cofr  + ftj 

where  (^A  is  the  phase  angle  between  roll  and  the  external 
torque  M and  introducing  this  value  of  9 and  the  values  of 

o o m 

its  first  and  second  derivatives  (9  and  Q in  the  equation  of 
motion,  will  yield 


(9= 


M 


V(C-Zvtoy 


and 


r 1 


~ coB 


C-  Ivco*- 


(2.2.42) 


by: 


The  expression  of  the  roll  response  will  then  be  given 


Q. 


A 


\/(fL-  (Pcf 


\Z(c-rvu)l)T-+  (cob) 


(2.2.43) 


The  response  amplitude  operator  being  the  ratio  of  the 
roll  amplitude  by  the  wave  amplitude  ^ , will  thus  be  given 
by: 


KAO.  = 


(2.2.44) 
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2.  2.  2.  7 Phage  relationship  between  roll  and  wave 

As  previously  established,  the  phase  angle  0^  between 
wave  and  exciting  torque  is  given  by: 


/l  . -I 


(T  = 


U) 


J>00 


(2.2.45) 


The  phase  angle  between  the  exciting  torque  and  the 
roll  response  is  in  turn  given  by: 


- ttz* 


-/ 


60S 


(2.2.46) 


C-  rvaf 

Finally  the  phase  angle  between  wave  and  roll  response 
will  be  the  sum  of  the  two,  i.e. , 


^5  = 'T-f  Kf) 


(2.2.47) 


3.0  COMPUTER  PROGRAMS 


3.  1 


3.  1.  1 


< 


I 


Heave  Computer  Program.  (HERAO) 

Program  logic 

The  operations  performed  by  the  heave  computer  program 
can  be  summarized  as  follows: 

- It  computes  the  heave  Response  Amplitude  Operator* 
using  formula  (2.2.23)  for  decreasing  values  of  the 
wave  angular  frequency  (j  • 
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The  recurrence  formula  used  to  change  the  value  of  the 
angular  frequency  between  two  consecutive  computations 
of  the  RAO  is 

cot.  = — &n: 

"In- 1 *♦*  At 

where  AT,  the  change  in  wave  period  is  constant. 

The  value  of  AT  to  use  for  a particular  set  of  compu- 
tations is  left  as  a program  input,  and  so  is  the  range 
of  variation  of  wave  periods  to  be  considered. 

- It  computes  the  phase  angles  between  force  and  wave, 
heave  and  force,  and  heave  and  wave  using  formulas 
(2.2.24),  (2.2.25),  and  (2.2.26),  for  the  same  set  of 


angular  frequencies 


- It  computes  the  wave  amplitudes  spectral  density  using 
one  of  the  spectral  density  formulas  (2.  1.  11),  (2.  1.  12), 
or  (2.  1.  13)  for  the  same  set The  choice  of  spectral 
density  formula  is  left  as  an  input. 

- It  computes  the  heave  response  spectral  density  KCco) 
using  formula  (2.  1.8)  and  the  computed  values  of  the 


RAO  for  the  selected  set 


H 


It  computes  the  root  mean  square  values  of  the  wave 
amplitudes  and  of  the  heave  response  amplitudes  by 
taking,  as  suggested  by  formulas  (2.  1.  5)  and  (2. 1.  9), 
the  square  root  of  the  area  under  the  wave  and  heave 
amplitudes  spectral  density  curves  established  for  the 
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- It  vises  the  statistical  results  of  Tables  I and  II  and  the 
two  root  mean  square  values  of  wave  amplitudes  and 
heave  response  to  compute  the  corresponding  expecta- 
tions of  wave  and  heave  means  and  maxima. 

3.1.2  Program  input 

The  program  is  designed  to  handle  any  number  of  cases  in 
consecutive  order.  All  input  is  format  free.  Values  for  any 
parameters  are  entered  in  a continuous  string,  separated  by 
commas.  The  program  is  designed  to  run  either  in  the  batch 
mode  or  interactively  from  an  on-line  remote  terminal.  The 
method  of  input  is  the  same  for  either  case.  As  the  interactive 
mode  is  also  self-explanatory  and  types  user  prompts,  the 
input  will  be  discussed  for  the  batch  mode.  All  depths  are 
considered  as  positive  downwards.  An  equivalent  depth  is  an 
average,  or  more  accurately  an  effective,  depth  at  which  a 
body  or  surface  is  located. 

Input  data  must  be  provided  on  the  following  cards: 

Card  1 --  Number  of  pressure  surfaces. 

NP  an  integer  value,  starting  in  column  1,  used  to 
specify  the  number  of  horizontal  pressure 
surfaces. 

Card  2 --  Pressure  depth,  surface  area. 

There  will  be  as  many  card  2's  as  specified  on 
card  1.  Each  will  contain  the  following  informa- 


tion. 
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DEPTHP  a value  specifying  the  "equivalent  depth" 
of  the  ith  pressure  surface,  in  feet, 

AREA  a value  specifying  the  area  of  the  pressure  sur- 
face, in  square  feet,  A negative  value  is  entered 
for  a surface  that  has  exerted  upon  it  a downward 
force.  A positive  value  is  entered  for  surface 
subjected  to  an  upward  force. 

Card  3 --  Number  of  inertial  components. 

NE  an  integer  value,  starting  in  column  1,  specifying 
the  number  of  inertial  components  which  comprise 
the  buoy. 

Card  4 --  Depth,  added  mass  coefficient,  volume. 

There  will  be  as  many  card  4's  as  specified  on 
card  3.  Each  card  will  contain  the  following  infor- 
mation. 

DEPTHI  a value  specifying  the  equivalent  depth  of  the  ith 
inertial  component,  in  feet. 

ADDMSC  a value  specifying  the  added  mass  coefficient 
for  the  ith  inertial  component. 

VOLUME  a value  specifying  the  volume  of  the  ith  inertial 
component,  in  cubic  feet. 

Card  5 --  Number  of  drag  surfaces. 

ND  an  integer  value,  starting  in  column  1,  specifying 
the  number  of  drag  surfaces  of  the  buoy  body. 
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Card  b --  Drag  depth,  linearized  damping  coefficient, 
linearized  wave  drag  coefficient. 

There  will  be  as  many  card  6's  as  there  are  drag 
surfaces  specified  on  card  5.  Each  card  will  con- 
tain the  following  information. 

DEPTHD  a value  specifying  the  equivalent  depth,  in 
feet,  of  the  ith  drag  surface. 

DAMPC  a value  specifying  the  linearized  damping 
coefficient  of  the  ith  drag  surface,  in  lbs 
force/ (ft /sec)/(  rad/ sec). 

WDRAGC  a value  specifying  the  linearized  wave  drag 
coefficient  of  the  ith  drag  surface,  in  lbs 
force /(ft /sec) /(rad/ sec). 

Card  7 --  Cross  sectional  area  at  water  surface. 

CAREAWL  a value,  starting  in  column  1,  specifying  the 
cross  sectional  area  at  the  water  line,  in  square 
feet.  For  the  purposes  of  this  analysis  this  area 
is  assumed  to  be  constant  over  the  range  of  motion 
at  the  water  line. 

Card  8 - - Virtual  mass. 

VIRTMASS  a value  specifying  the  virtual  mass  of  the 
body,  in  slugs. 

Card  9 (3F.  0)  --  Starting,  ending,  increments  of  wave  periods. 

TIME1  a value  specifying  the  lowest  wave  period  to  be 


studied,  in  seconds. 
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TIME2  a value  specifying  the  highest  wave  period  to  be 
studied,  in  seconds.  This  time  should  be  an  in- 
tegral multiple  of  the  incremental  time  (TIMEDEL) 
greater  than  TIME1. 

TIMEDEL  a value  specifying  the  incremental  wave  period, 
in  seconds,  used  in  the  analysis  from  TIME1  to 
TIME2. 

Card  10  (I,F.O,F.O)  --Amplitude  spectrum  selection,  para- 
meters. 

This  card  has  a general  form  as  follows. 

ISEASEL,  PARAMA,  PARAMB 
The  necessity  and  meaning  of  the  parameters  will 
depend  on  the  amplitude  spectrum  (ISEASEL) 
selected.  In  reality  these  are  double  height  formulas 
which  are  converted  internally  to  give  the  amplitude 
spectrum.  The  following  options  are  available. 
ISEASEL  = 1 Pierson-Moskowitz  formula. 

PARAMA  = wind  speed,  in  Knots 
no  PARAMB 

ISEASEL  = 2 Bretschneider  formula. 

PARAMA  = significant  wave  height,  in  feet 
PARAMB  = significant  wave  period,  in  seconds 
ISEASEL  = 3 International  Ship  Structure  Congress 
PARAMA  = significant  wave  height,  in  feet 
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PARAMB  = significant  wave  period,  in  feet 
Card  11  (A  1,  IX,  A1,1X,  2E.0)  --  Selection  of  listing,  line  printer 
plot,  plot  scale  minimum,  plot  scale  maximum. 

ILIST  enter  a Y (for  yes)  in  column  1 if  you  desire  a 
listing  of  the  various  output  parameters.  Any 
other  character  in  column  1 will  not  produce  a 
listing. 

IPLOT  enter  a Y (for  yes)  in  column  3 if  you  desire  a 
line  printer  plot  of  the  RAO.  Any  other  char- 
acter in  column  3 will  not  produce  a plot.  The 
line  printer  will  plot  a point  at  each  selected  wave 
period.  A check  is  made  on  the  length  of  the  plot 
for  the  following  criteria. 

N - T2  - T1  + TIMEDEL  < 2SQ 
TIMEDEL 

Note:  This  limitation  is  computer  dependent. 

RAOMIN  if  a plot  is  desired,  you  may  enter,  beginning 

in  column  5,  the  minimum  value  for  the  RAO 

scale.  If  left  blank,  RAOMIN  = 0. 
f 

RAOMAX  if  a plot  was  selected,  the  maximum  value  of 
the  RAO  scale  may  be  entered  following  RAOMIN 

(separated  by  a comma).  Under  the  current 

i 

version  RAOMAX  > 5.  For  best  results  RAOMAX- 
RAOMIN  should  be  an  integral  multiple  of  5,  If 


*** 
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it  is  not,  the  program  adjusts  it  to  be  so.  The 
current  default  is  RAOMAX  = 5. 

Card  12  (Al)  --  Another  case? 

IEND  If  you  wish  to  run  another  case,  enter  a Y or 

YES  beginning  in  column  1.  Any  character  other 
than  a Y in  column  1 will  cause  the  program  to 
terminate. 

The  sequence  of  card  types  1 through  12  is  repeated  for 
each  additional  case  desired. 

There  is  a special  entry  mode  for  additional  cases.  Be- 
cause the  buoy  configuration  may  be  quite  complex,  it  is  un- 
desireable  to  enter  all  the  descriptive  parameters  if  all  that  is 
changing  is  the  wind  speed  for  the  sea  state.  Another  alternative 
is  that  all  the  parameters  may  remain  constant  except  the  inertia 
terms.  A special  input  code  of  -1  will  allow  the  user  to  keep  in 
effect  the  values  last  entered  for  any  of  the  parameters.  This 
input  code  may  be  used  for  any  of  the  following  input  cards. 

Card  1 
NI  = -1 

Use  the  pressure  parameters  from  the  previous  case. 

Do  not  input  any  type  2 cards. 

Card  3 
NI  = -1 

Use  the  inertia  parameters  from  the  previous  case. 

Do  not  input  any  type  4 cards. 


h 
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Card  5 
ND  = -1 

Use  the  drag  parameters  from  the  previous  case. 

Do  not  input  any  type  6 cards.  ^ 

f 

Card  7 

CAREAWL  = -1 

Use  the  previously  entered  value  of  the  cross  sectional 
area. 

Card  8 

VIRTMASS  = -1 

Use  the  previously  entered  value  of  the  virtual  mass. 
Card  9 

TIME  1 = -1 

Use  the  previously  entered  time  range  and  increments. 
Card  10 

ISEASEL  = -1 

Use  the  previously  entered  amplitude  spectrum  for  the 
sea  state. 

Usage  modes 

As  previously  mentioned,  the  program  can  be  used  in  either 
a batch  or  an  interactive  mode. 

The  control  card  sequence  necessary  to  compile,  load,  and 
run  the  HERAO  program  in  a batch  mode  is  as  follows: 


! JOB  aaa,  uuu 

! LIMIT  (TIME,  3),  (CORE,  10) 
! FORTRAN  LS,  GO 


f!tP#SP !*¥« 


- 38  - 


FORTRAN  source  deck  of  program  HERAO 


IMETASYM  SI,  LO,  GO 


source  deck  for  subprogram  LAMTERM 


! LOAD  (GO),  (UNSAT),  (3)),  (MAP),  (LDEF),  (LMN,  HERAOR), 
(PERM) 

! RUN  (LMN,  HERAOR) 

(DATA 


data  cards  for  each  case 


! EOD 

To  run  a subsequent  job  utilizing  the  existing  load  module, 
only  the  following  cards  need  to  be  submitted: 

! JOB  aaa,  uuu 

! LIMIT  (TIME,  2),  (CORE,  10) 

! RUN  (LMN,  HERAOR) 

! DATA 


data  cards  for  each  case 

! EOD 

To  use  the  run  module  in  an  interactive  mode  from  a terminal, 
simply  log  on  and  enter,  to  a ! prompt,  the  following: 

START  HERAOR 

where  HERAOR  is  the  name  of  a previously  created  load  module. 

From  this  point,  the  operator  simply  responds  to  the  pro- 
gram prompts  as  if  you  were  punching  up  the  cards.  The  only 
difference  is  that  the  operator  does  not  need  to  start  in  column  X, 
but  should  start  as  though  the  head  were  already  positioned  cor- 
rectly. It  is. 
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Restrictions 


The  user  should  not  enter  a wave  period  of  zero  (0.0)  or 
less.  While  the  program  will  still  run,  the  integral  of  the  wave 
amplitudes  from  an  angular  frequency  of  oo  will  be  unreasonable. 
This  in  turn  will  cause  the  resulting  wave  statistics  to  be  in 
error.  All  other  parameters  computed  should  be  satisfactory. 

The  method  of  integration  used  is  that  of  a trapezoidal  approx- 
imation. The  user  must  therefore  exercise  some  care  in 
selecting  the  time  period  increment.  Too  large  an  increment 
may  cause  the  peak  of  the  RAO  or  wave  amplitude  spectrums  to 
be  "smoother",  resulting  in  lower  values  for  the  integrals  of  the 
he,  e response  and  amplitude  spectrum. 

Subprograms  required 

IAMTERM  a metasymbol  subprogram  which  checks  to  see  if 

the  program  is  being  run  in  batch  or  from  an 
on-line  terminal, 

computes  the  double  height  density  spectrum  for 
the  sea  state  according  to  one  of  several  empirical 
formulas;  internal. 

initializes  the  line  printer  plot  routine;  internal, 
executes  the  line  printer  plot  routine;  internal. 


SEAS PEC 


PLOTINJT 

PLOTHEAV 


LPPLOT  (PLOT1,  PLOT2,  PLOT3,  PLOT4,  PLOT5,  PLOT7) 


a subprogram  which  helps  create  and  list  a line 
printer  plot;  fromW.H.O.I.  account  3 library. 
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Program  output 

The  output  of  the  program  is  comprised  of  four  basic  parts. 
These  are: 

1.  Summary  of  input  parameters. 

2.  Summary  of  RAO,  phases,  and  amplitude  spectrum, 
all  given  as  a function  of  time  and  frequency. 

3.  Tabular  summaries  of  wave  and  heave  response 
statistical  properties. 

4.  Line  printer  plot  of  the  RAO. 

The  summary  of  the  input  parameters  is  only  given  for  a 
run  made  in  the  batch  mode.  For  an  on-line  hard  copy  terminal, 
the  users  entries  constitute  the  input  summary. 

The  listing  of  the  RAO,  phases,  and  other  information  is 
optional,  as  specified  in  column  1 of  input  card  type  11.  The  list 
has  the  same  format  whether  in  the  batch  or  on-line  mode.  Note 
that  the  amplitude  spectrum  is  output  for  the  sea  state. 

Tabular  summaries  of  the  wave  and  heave  response  statis- 
tical properties  are  always  output  and  are  the  same  regardless  of 
the  mode  of  operation. 

The  line  printer  plot  is  also  optional,  as  specified  in  column 
3 of  input  card  type  11. 

Typical  program  outputs  for  the  batch  and  terminal  modes 


are  shown  under  "Case  Studies",  Section  4. 
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Errors  and  diagnostics 

♦♦♦NUMBER  OF  ENTRIES  IS  GREATER  THAN  ARRAY 
SIZE  ALLOWS 

nrij  nn^ 

THE  PROGRAM  TERMINATES 

The  input  for  the  number  of  components  describing  the  buoy- 
configuration  exceeds  the  array  size  allocated.  Currently  nn2  - 20. 

♦♦♦THE  PLOT  BUFFER  IS  NOT  LARGE  ENGOUGH  FOR 
THE  PERIOD  RANGE  SPECIFIED 
THE  PLOT  IS  SUPPRESSED 

The  number  of  wave  periods  analyzed  must  meet  constraints 
described  in  the  input  section,  card  11. 

Timing 

The  program's  execution  time  is  a function  of  the  buoy  con- 
figuration and  the  number  of  wave  periods  analyzed.  In  any  case, 
the  execution  time  normally  is  negligible,  being  about  3 seconds 
(0.05  minutes)  per  case. 

3.  2 Roll  Computer  Program  (ROLLRAQ) 

3.  2.  1 Program  logic 

The  operations  performed  by  the  roll  computer  program 
are  similar  to  those  performed  by  the  heave  program.  They 
include: 

Computation  of  the  roll  RAO,  using  formula  (2.2.44) 
over  the  set  M previously  defined.  The  roll  RAO 
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is  expressed  in  units  of  degrees  of  roll  per  foot  of 
wave  amplitude.  The  recurrence  formula  to  change 
the  value  of  the  angular  frequency  between  two  consecu- 
tive  computations  of  the  RAO  is  again 

X-f+'AT 

where  Zi  I * the  change  in  the  wave  period  is  a constant 
set  by  the  user. 

Computation  of  the  phase  angles  between  external  torque 
and  wave,  roll  response  and  torque,  and  roll  and  wave 
using  formulas  (2,2,45,  46,  47)  for  the  same  set  {H}- 

- Computation  of  the  wave  amplitude  and  roll  response 
spectral  densities  S(cj)  and  K&>)  and  of  the  root  mean 
square  values  of  wave  amplitudes  and  roll  amplitudes. 
The  choice  of  spectral  density  formula  is  left  as  a pro- 
gram input. 

- Finally,  computations  of  expectations  of  means  and 
maxima  of  wave  and  roll  amplitudes  with  the  help  of 
the  statistical  parameters  shown  in  Tables  I and  II. 

3.2.2  Program  input 

The  program  is  designed  to  handle  any  number  of  cases. 
Almost  all  input  is  format  free.  Values  for  any  parameter  are 
entered  in  a continuous  string,  separated  by  commas  or  blanks. 
The  program  is  designed  to  run  either  in  the  batch  mode  or  in- 
teractively from  an  on-line  remote  terminal.  As  the  interactive 
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mode  is  also  self-explanatory  and  types  user  prompts,  the 
input  will  be  discussed  for  the  batch  mode.  All  depths  are 
considered  as  positive  downwards. 

Card  1 (3F,  0)  --  Period  range  of  time. 

TIME1  a value  specifying  the  lowest  wave  period  to 
be  studied,  in  seconds.  As  the  wave  velocity 
expression  contains  an  exponent  with  wave  fre- 
quency, the  user  is  cautioned  against  usu.g  a 
starting  period  of  less  than  1.0  seconds. 

TIME2  a value  specifying  the  highest  wave  period  to  be 
studied,  in  seconds.  This  value  of  time  should 
be  an  integral  multiple  of  the  incremental  time 
(TIMEDEL)  greater  than  the  value  of  TIME1, 

TIMEDEL  a value  specifying  the  incremental  wave 

period,  in  seconds,  used  in  the  analysis  from 
TIME1  to  TIME2. 

Card  2 (I,F.0,F.0)  --Amplitude  spectrum  selection,  parameters. 
This  card  has  a general  form  as  follows. 

ISEASEL,  PA  RAMA,  PARAMB 
The  necessity  and  meaning  of  the  parameters  will  depend 
on  the  amplitude  spectrum  (ISEASEL)  selected.  In 
reality,  these  are  double  height  formulas  which  are 
converted  internally  to  produce  the  amplitude  spectrum, 
the  following  options  are  available. 
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ISEASEL  = 1 Pierson -Mo  skowitz  formula 
PARAMA  wind  speed,  in  knots, 
no  PA  RAM  B 

ISEASEL  = 2 Brets chneider  formula. 

PARAMA  significant  wavs  height,  in  feet. 

PARAMB  significant  wave  period,  in  seconds. 

ISEASEL  * 3 International  Ship  Structure  Congress 
PARAMA  significant  wave  height,  in  feet. 

PARAMB  significant  wave  period,  in  seconds. 

Card_3  (F.  0)  - Radius  of  buoy  at  water  surface  plane. 

RWL  a value  specifying  the  outer  radius  of  the 

buoy  at  the  surface  (still  water  assumed),  in  feet. 
Card_4  (F. 0)  . . Depth  to  keel. 

DEPTHK  a value  specifying  the  depth  to  the  keel  (bottom) 
of  the  buoy,  in  feet. 

~rd  ~ (F,0)  ““  Average  wave  amplitude. 

AVERGAMP  a value  soer-i fvina 

specifying  the  average  expected 

wave  amplitude,  in  feet. 

~rd  6 (Fa0>  — Average  roll  constant. 

THETABAR  a value  specifying  the  average  expected 

roll,  in  degrees. 

— -d  7 (I)  ~ Number  of  buoy  components. 

NP  a value  specifying  the  number  of  buoy  com- 

ponents. This  may  be  set  to  zero  in  subsequent 


cases. 


Cards  8 ...  (7  + NP)  (I,  6F.  0,  3A4, A2)  --  Component  specification. 
ISHAPE  an  integer  value  used  as  a code  to  specify 
the  component  shape. 

1)  hollow  cylinder 

2)  solid  cylinder 

3 ) solid  disc 

4)  right  triangular  plate 

WIDTH  a value  specifying  the  width  (diameter  or 
base)  of  the  component,  in  feet. 

HEIGHT  a value  specifying  the  height  of  the  component, 
in  feet. 

THICK  a value  specifying  the  thickness  of  the  com- 
ponent, in  inches.  Entering  THICK  = -1  for 
ISHAPE  = 1 will  generate  a solid  (THICK  = 
WIDTH/2.0).  For  ISHAPE  = 2 or  3,  also 
enter  a -1  as  it  is  ignored. 

DENSITY  a value  specifying  the  density  of  the  com- 
ponent, in  pounds  mass  per  cubic  foot 
(lbsm/ft^). 

DISTCG  a value  specifying  the  vertical  distance  from 
the  buoy  keel  to  the  component  center  of 
gravity,  in  feet.  The  component  c.g.  is 
vertical  vector  only. 
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FRACNORM  a value  specifying  the  fractional  propor- 
tion of  the  component  area  normal  to  the  roll 
motion.  For  cylinders  this  entry  = 1.0.  For 
triangular  plates  it  will  vary  from  0.0  (oriented 
in  line  with  roll)  to  1.0  (area  normal  to  roll 
motion). 

ICOMMENT  a character  string  used  to  describe  the 
component. 

Note:  For  the  purpose  of  visually  inspecting  data  cards  used 
in  batch  input,  it  may  be  desirable  to  "format"  the 
data.  The  recommended  format  is  (15,  6F10.  0,  IX, 
3A4.A2). 

Card  9 (I)  --  Redefined  part  code. 

N an  integer  value  specifying  the  number  (index) 

of  the  component  to  be  redefined.  This  allows 
the  user  to  change  the  dimensions  of  a com- 
ponent^) or  to  add  new  components.  To 
add  a new  component,  the  specified  value  of 
N must  be  one  greater  than  the  current  maxi- 
mum number  of  parts  defined  for  the  buoy. 

There  may  be  as  many  redefinition  pairs  of 
cards  as  desired.  To  terminate  the  sequence 


enter  a value  for  N = -1. 
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Card  10  (I,  6F.  0,  3A4,  A2)  --  Redefinition  specification. 

This  card  has  the  same  format  as  card  7,  the  component 
specification.  Any  number  of  pairs  of  this  and  the  pre- 
vious card  may  be  entered  as  desired. 

Card  11  (Al,  IX,  Al,  IX,  2F.  0)  --  Selection  of  listing,  line 
printer  plot,  plot  scale  minimum  and  maximum. 
ILIST  enter  a Y (for  yes)  in  column  1 if  you  desire 

a listing  of  the  various  output  parameters. 
Any  other  character  in  column  1 will  cause 
the  listing  to  be  suppressed. 

IPLOT  enter  a Y (for  yes)  in  column  3 if  you  desire 

a line  printer  plot  of  the  roll  RAO.  Any  other 
character  in  column  3 will  cause  the  line 

t 

printer  plot  to  be  suppressed.  The  line 
printer  plots  a point  at  each  selected  wave 
period.  A check  is  made  on  the  length  of  the 
plot  for  the  following  criteria. 


T2  - T1  + TIMEDEL  < 100 
TIMEDEL  “ 

RAOMIN  if  a plot  is  desired,  you  may  enter,  begin- 
ning in  column  5,  the  minimum  value  for  the 
RAO  scale.  If  left  blank  RAOMIN  = 0.0. 

RAOMAX  if  a plot  was  selected,  the  maximum  value 
of  the  RAO  scale  may  be  entered  following 
the  RAOMIN  value  (separated  by  a blank  or 
comma).  For  the  best  results  RAOMAX  - 


f 


RAOMIN  should  be  an  integral  multiple  of 
5.  If  it  is  not,  the  program  adjusts  it  to  be 
so.  The  current  default  is  RAOMAX  = 10,0. 


Card  12  (Al)  --  Another  case? 


IEND 


If  you  wish  to  run  another  case,  enter  a Y or 


YES  beginning  in  column  1.  Any  character 
other  than  a Y in  column  1 will  cause  the  pro- 
gram to  terminate. 

The  sequence  of  card  types  1 through  12  is  repeated  for 
each  additional  case  desired. 

There  is  a special  entry  mode  for  additional  cases.  Be- 
cause the  buoy  configuration  may  be  quite  complex,  it  is  undesir- 
able to  enter  all  the  descriptive  parameters  to  investigate  the 
effect  of  a different  wind  speed  or  sea  state.  Another  alternative 
is  that  the  design  of  the  buoy  remains  the  same  except  for  the 
size  of  the  counterweight.  A special  input  code  of  -1  and  the 
component  redefinition  options  allow  the  user  to  keep  in  effect 
the  values  last  entered  for  any  of  the  parameters  or  components. 

Cards  1-7  A -1  will  maintain  those  values  previously 
entered.  For  card  type  7,  this  simply  uses 
the  same  buoy  and  no  new  components  are 
defined  (no  type  8 cards). 

Cards  9-10  As  many  of  these  pairs  as  desired  may  be 
entered  to  define  a new  buoy  configuration. 


The  entry  sequence  is  terminated  by  a -1 
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for  card  type  9 and  no  card  type  10. 

- * Hia&g 

As  previously  mentioned,  this  program  is  meant  to  be 
used  in  either  a batch  or  an  interactive  mode.  The  following 
control  card  sequence  is  a complete  list  of  steps  necessary  to 
compile,  load,  and  execute  the  program  ROLLRAO  in  the 
batch  mode. 

! JOG  aaa,  uuuu 

! LIMIT  (TIME,  3),  (CORE,  20) 

! FORTRAN  LS,GO 

Fortran  source  deck  of  program  ROLLRAO 
JMETASYM  SI,LO,GO 

source  deck  of  subprogram  IAMTERM 
! LOAD  (GO),  (UNSAT,  (3)1  (MAP),  (LDEF),  (LMN,  ROLLR),  (PERM) 
! RUN  (LMN,  ROLLR) 

(DATA 

data  cards  for  all  cases 
! EOD 

To  run  a subsequent  job  utilizing  the  existing  load  module, 
only  the  following  cards  need  to  be  submitted: 

! JOB  aaa,  uuuu 

! LIMIT  (TIME,  2),  (CORE,  20) 

(RUN  (LMN,  ROLLR) 

! DATA 

data  cards  for  all  cases 
! EOD 

To  run  the  load  module  from  a terminal,  simply  log  on  and 
enter,  to  a ! prompt,  the  following: 

START  ROLLR 

From  this  point  on  you  simply  respond  to  the  program 
prompts  as  if  you  were  punching  the  data  cards.  The  only  dif- 


ference is  that  you  do  not  need  to  start  in  column  1,  but  should 


- 50  - 


start  as  if  you  are  already  there.  (There  is  also  a slight 
difference  in  the  manner  in  which  the  buoy  components  are 
described  and  entered.  This  is  explained  to  the  user  on  line.  ) 
Restrictions 

The  user  should  exercise  care  in  entering  wave  periods  of 
less  than  1.0  seconds.  The  water  velocity  terms  contain  a 
natural  exponent  of  the  wave  number  times  displacement  from 
the  center  of  gravity  and  this  computation  may  exceed  the 
machine  computational  capability. 

The  method  of  integration  used  is  that  of  a trapezoidal 
approximation  over  frequency.  The  user  must  therefore  exer- 
cise additional  care  in  selecting  the  time  period  increment.  Too 
large  an  increment  may  cause  the  peak  of  the  RAO  or  wave 
amplitude  spectrums  to  be  "smoothed",  resulting  in  lower 
values  for  the  integrals  of  the  amplitude  spectrum  and  roll 
response.  At  the  other  extreme,  the  time  period  of  0.0  will 
cause  an  "infinite"  angular  frequency  and  the  resulting  wave 
statistics  will  be  in  error. 

Also  note  that  the  entire  program  is  executed  in  single 
precision.  For  normal  buoy  configurations  this  mode  is 
adequate.  However,  in  certain  cases  such  as  the  case  of  a 
small  flat  cylinder  (case  study  no.  1,  Section  4.  1),  the  response 
amplitude  operator  exhibited  signs  of  instability  in  the  numerical 
computation.  This  can  be  overcome  by  using  double  precision 
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computations. 

Subprograms  Required 

IAMTERM  a Metasymbol  subprogram  which  checks  to 
see  if  the  program  is  being  run  in  batch  or 
from  an  on-line  terminal. 

SEASPEC  computes  the  amplitude  density  spectrum 

for  the  sea  state  according  to  one  of  several 
empirical  formulae;  internal. 

PLOTINIT  initializes  the  line  printer  plot  routine; 
internal. 

PLOTROLL  executes  the  line  printer  plot  routine;  internal. 

TINPUT  prompts  and  inputs  data  from  a user  on-line. 

BIN  PUT  reads  and  summarizes  data  entered  in  batch 

mode. 

BODYVOL  computes  the  volume  of  a component. 

BODYMI  computes  the  shape  dependent  contribution  of 

a component's  moment  of  inertia. 

DISPLACE  computes  those  parameters  associated  with 
the  displacement  of  the  buoy's  components. 

BUOYDAMP  (WATERDAMP)  computes  the  moments  of 

damping  for  the  buoy  and  water  drag  forces. 

WATRINRT  computes  the  inertia  moment  contribution 


from  the  water  wave  particle  acceleration 
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LPPLOT  (PLOT1,  PLOT2,  PLOT3,  PLOT4,  PLOT5,  PLOT7) 
a subprogram  which  creates  and  teats  a 
line  printer  plot;  W.H.O.I.  account  3 
library. 

3.  2.  3 Program  output 

The  output  of  the  program  is  comprised  of  four  basic 
parts.  These  are: 

1.  Summary  of  input  parameters  and  physical  properties 
of  the  buoy. 

2.  Summary  of  RAO,  phases,  and  amplitude  spectrum, 
all  given  as  a function  of  time  frequency. 

3.  Tabular  summaries  of  wave  and  roll  response  statis- 
tical properties. 

4.  Line  printer  plot  of  the  roll  RAO. 

The  summary  of  the  input  parameters  is  only  given  for  a 
run  made  in  the  batch  mode.  For  an  on-line  hard  copy  terminal, 
the  user  entries  constitute  the  input  summary. 

The  listing  of  the  RAO,  phases  and  other  information  is 
optional,  as  specified  in  column  1 on  input  card  type  11.  The 
list  has  the  same  format  whether  in  batch  or  on-line  mode.  Note 
that  the  amplitude  spectrum  is  output  for  the  sea  state.  This 
is  derived  directly  from  the  selected  double  height  spectrum. 


Tabular  summaries  of  the  wave  and  roll  response  statis- 
tical properties  are  always  output  and  are  the  same  regardless 
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of  the  mode  of  operation. 

The  line  printer  plot  is  also  optional,  as  specified  in 
column  3 of  the  input  card  11. 

Timing 

The  program's  execution  time  is  a function  of  the  buoy 
configuration  and  the  number  of  wave  periods  analyzed.  Nor- 
mally the  execution  time  is  much  less  than  a minute  per  case. 
Errors  and  Diagnostics 

A number  of  checks  are  made  on  the  input  parameters  and 
the  stability  of  the  buoy.  The  messages  indicate  the  nature  of 
the  error  and  take  appropriate  action  depending  on  the  execution 
mode. 

4.  0 CASE  STUDIES 

To  illustrate  the  use  of  the  computer  solutions  the  following 
three  case  studies  are  presented. 

a.  Heave  and  roll  response  of  a flat  cylinder  of  small 
dimensions. 

b.  Heave  and  roll  response  of  a ballasted  "telephone  pole.  " 

c.  Heave  and  roll  response  of  a complex  shape  W.H.O.I. 
spar  buoy. 

4.  1 Heave  and  Roll  Response  of  a Small  Flat  Cylinder 

This  case  study  has  relatively  little  practical  value.  Every 
one  knows  that  a thin  slice  of  pulpwood,  if  thrown  in  the  sea, 
will  essentially  follow  the  heave  and  slope  of  the  waves,  both 
in  magnitude  and  phase.  It  is  included  here  mainly  for  the 
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purpose  of  program  result  verification. 

The  small  cylinder  considered  has  a diameter  of  one 
foot.  The  ratio  of  its  height  to  its  diameter  is  1:3.  Its 
density  is  3/4  that  of  sea  water.  The  heave  and  roll  motion 
of  this  cylinder  will  be  studied  using  a Pierson-Moskowitz 
spectrum  with  a 20  knots  wind. 

4.  1.  1 Program  input 

Buoy  draft.  The  buoy  draft  is  obviously  0.  25'. 

Pressure  surface.  The  pressure  force  is  exerted  ex- 
clusively on  the  cylinder  lower  face.  The  pressure  surface  is 
therefore  at  0.25'  from  the  surface.  Its  area  is 
IX  / 4 = 0.  785  sq.  ft.  The  pressure  force  being  upwards  «S 
is  positive. 

Inertial  component.  The  added  mass  effect  will  be  con- 
considered  to  take  place  at  the  lower  face  of  the  cylinder,  i.e. 
at  0.25'  from  the  surface.  It  will  be  assumed  to  take  place 
only  half  of  the  time,  during  the  downwards  part  of  the  heave 
cycle.  It  will  be  estimated  to  be  the  same  as  the  one  produced 
by  a flat  plate  of  same  radius  as  the  cylinder.  Thus  the  averaged 
added  mass  will  be: 

* ’ m0 ’ 33  ,lug' 

Assuming  the  added  mass  coefficient  to  be  equal  to  one,  the 


averaged  added  mass  coefficient  and  volume  will  be 
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Cm  - ^ <^1  = O.So 
Ml  _ -®^3  * ^.33  ft. 

3 

Damping/drag  surface.  Damping  and  wave  drag  will  also 

be  considered  to  take  place  at  the  lower  face  of  the  cylinder,  i.e 

at  0.25'  from  the  surface.  These  effects  will  again  be  assumed 

to  take  place  only  half  the  time. 

A choice  must  now  be  made  for  the  arbitrary  values  of 

average  heave  and  average  wave  amplitude  . Assum 

ing  that  the  heave  equals  the  wave  amplitude,  a fair  assumption 

in  this  case,  and  selecting  the  average  wave  amplitude  for  winds 

of  20  knots  to  be  V will  yield  /X  ^ — Xc  a 

Using  this  value  of  ^ anc*  *n  exPress*ons  (2.2.10) 

and  (2.  2.  13),  and  a drag  coefficient  « 0._9  , the  value  of 

I ' 1 

the  linearized  damping  and  wave  drag  coefficients  fc>  and  C is 
found  to  be:  , 

6'  . c!  = /.&  Ik-'  j 

//-  SeJ 

/ 

To  account  for  the  time  average  only  half  of  these  coefficients 
value,  i.e.  0.9,  should  be  used  as  program  input. 

Cross  sectional  area  at  surface.  This  area  equals  rL  /4 
or  0.  785  sq.  ft. 

Buoy  virtual  mass  '^*1  r 

/ft i - buoy  mass  = ~ V = 0.  392  slug. 

1 * * 4 31* 

t 

- added  mass  = 

Thus,  buoy  virtual  mass 


0,  333  slug 
= 0. 725  slug 
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In  addition  to  theoe  parameters,  the  user  must  also 
specify  the  arbitrary  roll  constant  (S)  and  wave  amplitude 
constant^  . In  this  case  the  average  angle  roll  <3  will 
be  assumed  to  equal  the  slope  of  the  average  wave  in  20  knots 
wind.  Assuming  an  average  amplitude  d of  3*  and  an  average 
wave  length  of  120',  the  average  wave  slope  is  found  to 
be: 


O Zr>s4  * JL  = 0.  157  radians 

“ ZT 

A summary  of  the  data  input  is  shown  in  the  Data  Coding 
Form  Fig.  No.  8 and  Fig.  No.  9. 

4.1.2  Program  output 

Heave  and  roll  response  amplitude  operators  are  depicted 
in  Fig.  Nos.  10  and  11. 

Computed  expectations  of  average  and  maximum  values 
of  wave  amplitudes  and  of  heave  and  roll  motions  are  also 
obtained  with  the  help  of  the  computer  programs.  Correspond- 
ing values  of  wave  slopes  are  calculated  independently  using 


^ j fj 


where  ^ is  the  applicable  Raleigh  constant 
and  51  > 

being  the  wave  number  = and 

being  the  wave  amplitude  spectral  density. 
These  results  are  summarized  in  Table  Nos.  3 and  4. 


The  heave  RAO,  with  a value  of  one  over  most  of  the 
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Fig.  No.  10 
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RMS  pF  WAVE  SpECT«U*  ■ 9.683  FEET 

PROS ARLF  AMPLITUDE 

*F  wAVr  ■ 1.897  FEET 


EXPECTFD 

WAVE 

MAXIMUM 

amplitude 


5*6«9 

6.118 

7*003 

7.459 

8.399 

9.311 


FRACT I ON  OF 
LARGEST 
AMPLITUDES 

considfrfd 

. DIO 
• 100 
.339 

.500 

1*000 


AVFPAQt 

*AV& 

AMPLITUDE 

6*330 

4.830 

3.799 

3.370 

2.377 


NUMBER 
OF  WAVES 

50 

100 

500 

1000 

10000 

100000 


RmS  OF  RESPONSE  SPECTRUM  • 


P.684  FEET 


PROBABLE  AMPLITUDE 

OF 

HEAVE  RESPONSF 

■ 

1 .897 

FEET 

FRACTION  of 

AVERAGE 

I 

expected 

LABQFST 

HEAVE 

! 

HEAVE 

AMPLITUDES 

amplitude 

I 

NUMBER 

MAXIMUM 

considered 

RESPONSE 

I 

f 

OF  WAVES 

amplitude 

*010 

6*331 

I 

I 

50 

5*690 

*100 

4.831 

I 

too 

6.119 

.333 

3.800 

I 

500 

7*005 

•500 

3.371 

I 

1000 

7.461 

1*000 

2.378 

I 

10000 

8.400 

I 

1DOOOO 

9.313 

Table  No.  3 


Heave  Response  of  Flat  Small  Cylinder 
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Amplitudes  of  Roll  are  in  Degrees 


Fraction  of 
Largest 
Amplitudes 
Considered 


.010 
. 100 
. 333 
. 500 
1.000 


Ave  rage 
Roll 

Amplitude 

Response 

25.070 
19.  129 
15.048 
13. 348 
9.416 


I 

I 

I Number 
I of  Waves 
I 

I 50 

I 100 

I 500 

I 1000 

I 10000 

I 100000 


Expected 

Roll 

Maximum 

Amplitude 

22. 530 
24. 230 
27. 737 
29. 544 
33.264 
36.877 


Amplitudes  of  Wave  Slope  are  in  Degrees 


Fraction  of 

I 

Largest 

I 

Amplitudes 

Average 

I 

Number 

Considered 

Slo  pe 

I 

of  Wave 

.010 

23. 905 

I 

I 

50 

. 100 

18.  240 

I 

100 

. 333 

14.  349 

I 

500 

.500 

12. 728 

I 

1000 

1.000 

8.978 

I 

10000 

I 

100000 

Expected 

Maximum 

Slope 

21.483 
23. 104 
26.448 
28. 171 
31.718 
35.  163 


Table  No.  4 


Roil  Response  of  Small  Flat  Cylinder 
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wave  periods  considered,  clearly  indicates  that  the  small 
cylinder  is  a perfect  wave  follower.  Averages  and  expected 
maximum  values  of  heave  when  compared  to  corresponding 
values  of  wave  amplitude  confirm  this  expected  result. 

The  roll  RAO,  being  in  degree  of  roll  per  foot  of  wave 
rather  than  per  degree  of  slope  does  not  immediately  correlate 
roll  and  wave  slope.  It  shows  simply  that  roll  is  large  at 
small  periods,  and  tends  to  zero  as  the  wave  period  increases, 
which  is  of  course  precisely  what  the  slope  of  the  wave  does. 

The  computed  statistical  averages  however  do  confirm  that 
for  practical  purposes  the  roll  angle  is  strongly  correlated  to 
wave  slope.  It  thus  appears  that  this  first  study  case  is  a good 
test  of  the  program  validity. 

4.  2 Heave  and  Roll  Response  of  a Ballasted  "Telephone  Pole" 

We  next  consider  a cylindrical  body  made  of  two  cylinders 
of  same  diameter,  but  of  different  lengths  and  densities,  as 
shown  on  Fig.  No.  12. 

The  density  of  the  large  cylinder  is  16  lbs/cu.ft.  The 
density  of  the  small  cylinder  is  384  lbs/cu.ft.  As  in  the  previous 
case,  the  heave  and  roll  response  of  this  body  will  be  studied 
using  a Pierson-Moskowitz  spectrum  with  a 20  knots  wind. 

4.  2.  1 Program  input 

In  order  to  provide  the  necessary  input  data  to  the  heave 
and  roll  programs,  the  following  computations  must  be  first 


performed. 


Fig.  No.  12 


Buoy  draft. 

weight  of  small  cylinder  = K xlxlx2x384  = 2413  lbs 
weight  of  long  cylinder  = ft  xlxlx22xl6  = 1105  lbs 

Total  Weight  = 3518  lbs 

Draft  = iii® = 17.  5' 

)f  xlxlx64 

Number,  depth,  and  area  of  pressure  surfaces.  In  thi a 
simple  case  there  is  again  only  one  pressure  surface,  namely 
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the  pole  lower  face.  It  is  located  17.  5'  from  the  surface.  Its 
area  " is 

\ 

* 

s n x /v  / j - 3.  14  sq.  ft. 

The  pressure  force  being  upwards,  is  positive. 

Number,  depth,  added  mass  coefficient  and  volume  of 
inertial  components.  In  heave  motion,  the  added  mass  effect 
will  be  considered  to  be  essentially  produced  by  the  pole  lower 
end.  There  will  thus  be  only  one  inertial  component,  acting  at 
17.  5'  from  the  surface.  The  added  mass  effect  will  be  estimated 
to  be  the  same  as  the  one  produced  by  a sphere  of  sami  radius 
as  the  pole  (long  cylinder  approximation),  but  acting  only  half 
of  the  time.  Thus  the  averaged  added  mass  will  be  expressed 

* H -i  (Ok  * ik  = *•  °& 

The  corresponding  added  mass  coefficient  and  volume 
are  therefore 

= ■£ 

VcL  >=  r 4.  (3  Cu^t 

Number,  depth,  damping  and  wave  drag  coefficients  of 
drag  surfaces.  Damping  and  wave  drag  effects  arc  assumed  to 
be  also  essentially  produced  by  the  pole  lower  face.  Thus  there 
will  be  only  one  drag  surface  located  at  17.  51.  These  effects 
will  be  assumed  to  take  place  only  half  of  the  time. 
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The  average  heave  will  be  assumed  to  be  half  of 

the  average  wave  amplitude  and  the  later  will  be  assumed 

to  be  3*.  Using  these  values  of  ^ and  X^in  expressions 
(2.2.10)  and  (2.2.  13)  together  with  Cl  =0.9,  yield 

t>'  = 5.6  lk,J_  I nod. 

' - 77 


/ 7 a ft-  'Wc 

c,  •«  (><s  r * / 
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Again,  to  account  for  the  fact  that  damping  effects  are 
assumed  to  occur  only  half  of  the  time,  only  half  of  the  co- 
efficients values,  i.e.  1.8  and  3.6  respectively,  should  be  used 
as  program  input. 

Cross  sectional  area  at  surface. 

n+t«i  « . 


Buoy  virtual  mass 

s .3«5/<S 

3*.* 

**}  =r 


109  slugs 
2.08  slugs 
111. 08  slugs 


In  addition  to  these  computed  parameters,  the  program 
user  must  select  the  arbitrary  average  roll  constant  Q and 
wave  amplitude  ft  . For  this  example  <9  is  set  equal  to 
5°  and  ft  equal  to  3.0  feet. 

Obviously  the  number  of  buoy  parts  is  two.  Their  charac- 
teristics are  summarized  as  follows: 


Part 

No. 

Name 

Sha  pe 

Width 

(ft) 

Height 

(ft) 

Thick 

(ft) 

Density 

(lbs/cu.ft) 

H 

1 

Upper 

Cyl. 

Solid 

Cyl. 

2 

22 

16 

13 

2 

Lower 

Solid 

Cyl. 

2 

2 

! 

384 

1 

1 

All  pertinent  input  data  are  listed  in  the  data  coding  form 
shown  in  Fig.  Nos.  13  and  14. 

4.2.2  Program  output 

Values  of  the  heave  and  roll  response  amplitude  operators 
of  the  telephone  pole  for  the  prescribed  period  interval  (50  seconds) 
and  increment  (1  second)  are  presented  in  the  typical  computer 
printouts  shown  in  Fig.  No.  15  and  Fig.  No.  16,  The  response 
amplitude  operators  are  also  graphically  represented  in  Fig. 

Nos.  17  and  18.  Expected  average  and  maximum  values  of 
wave  amplitudes  and  of  heave  and  roll  responses  are  as  tabu- 
lated in  Tables  5 and  6, 

4.  3 Heave  and  Roll  Response  of  a Complex  Shape  W.H.O.I. 

Spar  Buoy 

Fig.  No.  19  shows  the  dimension  and  shape  of  the  spar 
buoy  to  be  studied  next. 

Heave  and  roll  response  will  be  again  established  using  a 
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Pierson-Moskowitz  spectrum  with  a wind  of  20  knots. 

This  relatively  complex  shape  buoy  is  made  of  the  follow- 
ing parts: 

- A reserve  buoyancy  cylinder  3 ft.  high,  2.  5 ft.  in 
diameter,  made  of  2 lbs/cu.ft.  polyurethane  foam. 

- A spar  24  ft.  long  made  of  8"  O.  D by  1/4"  thick  wall 
aluminum  tubing. 

- A spar  base  plate,  made  of  4 ft.  diameter  by  1/2" 
thick  aluminum  plate. 

- A 3'.0"  diameter  by  8'.0"  long  buoyancy  tank  made  of 
3/l6"  steel  plates.  The  buoyancy  tank  is  filled  with 

4 lbs/cu.ft.  foam. 

- A 3'.0"  diameter  by  8'.0"  water  ballast  tank  made  of 
1/8"  steel  plates.  The  tank  is  filled  with  sea  water. 
The  bottom  plate  is  4'.0"  in  diameter. 

- A 10'.  0"  long  stem  made  of  6 5/8"  O.  D.  schedule  40 
steel  pipe  filled  with  sea  water. 

- A 4'.0"  diameter  by  1/2"  thick  damping  plate. 

- A counterweight  cylinder  2,5  ft.  in  diameter  by  0.848 
ft.high,  made  of  cast  iron. 

The  physical  parameters  of  the  buoy  main  components  are 


I 


summarized  in  Table  No.  7. 


Table  No 


Total  = 7773.84 
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4.  3.  1 Program  input 

Computations  and  considerations  in  the  support  of  pro- 
gram input  are  as  follows. 

Buoy  draft.  The  weight  of  the  water  displaced  by  the 
buoy  equals  the  weight  of  the  buoy. 

Weight  of  water  displaced  by  the  stem 

JL  (0  - 55Z.J  \ /Oy  _ 15 Z.  /£ 

T 

Weight  of  water  displaced  by  the  ballast  tank 

_/ix3x3y6V£4  c $r/e.  n 

H 

Weight  of  water  displaced  by  the  buoyancy  tank 

J3  x 3 x 3k  6 *64  =x3  £/3 , // 

& 

Weight  of  water  displaced  by  the  immersed  portion  of 
the  spar  of  length  " 

Jt  +(0.  66 6 J 4.  y CJf  r 22. 

Solving  for 

4 = 7/73.84  - +Jx3Xre  „}  g tz.o ' 

The  buoy  draft  is  therefore  38'. 

Number,  depth  and  area  of  pressure  surfaces.  There 
are  two  pressure  surfaces  to  consider,  namely  the  top  and  the 
bottom  of  the  foam  filled  buoyancy  tank.  The  area  of  the  first 
pressure  surface  St  is  given  by 

'-C  = K (/ft  2-  ffz  J where 
*7^  is  the  radius  of  the  tank,  Kj  =1.5' 
is  the  radius  of  the  spar,  - 0.  33' 

Thus  = ?T(1.  52  - 0.  332)  - 6.71  sq.  ft. 
is  located  1 2 ' below  the  surface.  The  pressure  force  acting 
on  *3^  being  downwards,  is  negative.  The  area  of  the  second 

• • --  «■  ■?*  ■**&!*: ' mmw******™ 
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pressure  surface  v5 £ is  in  turn  given  by 
SL  = = 7.  06  sq.  ft. 


It  is  located  20*  below  the  surface.  The  pressure  force  is 
upwards  and  thus  A is  positive. 

Number,  depth,  added  mass  coefficients,  and  volume  of 
inertial  components.  For  computing  the  heave  response,  two 
distinct  added  mass  effects  must  be  accounted  for:  the  added 
mass  due  to  the  water  entrained  by  the  top  and  bottom  plates 
of  the  buoyancy  and  water  ballast  tanks,  and  the  added  mass  due 
to  the  water  entrained  by  the  damping  plate. 

The  first  added  mass  effect  will  be  assumed  to  be  the 
same  as  the  one  produced  by  a sphere  with  a diameter  equal 
to  the  diameter  (4')  of  the  plates  located  at  the  top  of  the 
buoyancy  tank  and  the  bottom  of  the  water  ballast  tank.  The 
equivalent  depth  will  be  selected  midway  between  the  two 
plates,  i.e.  20'  below  the  surface.  The  added  mass  coeffi- 
cient for  a sphere  is  J.  . The  volume  VOL.  of  this  first 
inertial  component  is 

VoLt  — A ft  — 33. 57  CM-fe- 

The  second  added  mass  effect  will  be  considered  to 
take  place  at  the  buoy  keel,  i.e.  at  38'  below  the  surface.  The 
formula  for  the  added  mass  of  a circular  plate  of  radius 
being 


Ai> 1 


^ «ai?  « ^oL 
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An  arbitrary  added  mass  coefficient  of  1 will  yield  a volume  1 10L± 

i/oL^  = -§a.3=  2h  33 

Both  effects  in  this  case  are  happening  all  the  time. 

Number,  depth,  damping  and  wave  drag  coefficients  of 
drag  surfaces.  Damping  and  wave  drag  will  be  assumed  to  be 
produced  mainly  by  or  on  the  upper  and  lower  faces  of  the 
buoyancy  tank  and  water  ballast  tank  and  by  the  damping  plate 
at  the  buoy  lower  end.  There  will  thus  be  two  drag  surfaces, 
one  assumed  to  be  located  half  way  between  the  two  ends  of 
the  tanks  at  an  equivalent  depth  of  20',  and  the  other  at  the 
buoy  keel  38*  below  the  surface.  Assuming  ^ X,-_ 

and  X =3*  the  corresponding  damping  and  wave  drag 

c 

coefficients  given  by  expressions  (2.  2.  10)  and  (2.  2.  13)  are 


found  to  be 


/>y  sr  A y Z#  V'9  V X 

3/t. 

K ~jLxZy/'ZynM\  / S = 

C 'm  £ yl  y 5 
; *1 

c’«  4 xl  ylZ  3 

1 3/l  ' 


« /M* 
= Sf.£0 
- 28.60 
a 36.  40 


Cross  sectional  area  at  surface. 


Kx0-3%z  = 0.3^Z,  -'t 


f 
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■j  Buoy  virtual  mass. 

& 

Added  mass  of  first  inertial  component 

V = ~ 33‘3  slu®8 

' 2 3Z.Z.  3/t 

added  mass  of  second  inertial  component 
/>ff’  — &.x  x /z  ) ss  42.4  slugs 

*■  *3  3*  a. 

Buoy  mass 

/)yj  — 77JS.  _ 237.  85  slugs 

32  2.  t I 

Virtual  mass  ^ sr  /%  + /*>  + /Hi  s.  313.  55  slugs 

v ^ Z> 

As  in  the  preceding  case  study,  the  program  user  must 
also  provide  an  arbitrary  value  of  average  buoy  roll  © and 
wave  amplitude./^  . In  this  case  Q and^  are  selected  to 
be  5 ° and  2.5  feet  respectively. 

All  pertinent  data  are  listed  in  the  data  coding  form 
shown  in  Fig.  Nos.  20  and  21. 

4.  3.  2 Program  output 

The  heave  and  roll  response  amplitude  operators  are 
graphically  represented  in  Fig.  No.  22  and  Fig,  No.  23.  The 
expected  average  and  maximum  values  of  wave  amplitudes  and 
of  heave  and  roll  motion  are  summarized  in  Table  Nos.  8 and  9. 

As  a point  of  interest,  Table  No.  10  presents  a succinct 
performance  comparison  of  the  three  buoy  types  when  sub- 


mitted to  the  same  random  excitation. 


NO  l * A t 0 * V 
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RMS  8P  WAVE 
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PRObArLE  amplitude 
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CpNSlDpRgD 
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AVERAUE  1 

*AVS  I 

amplitude  i 

NUMB|P 

OP  ^a^es 

ejected 

WAVE 

Maximum 

amplitude 
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•100 

•331 

•500 

1*000 
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2*377  I 

I 

so 

100 

500 

1000 

10000 

100000 
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amplitudes 
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AVEHaQe  I 
weave  I 

rM«»bn«°e  ' 
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op  WA&ES 

E*PECteD 

WEAVE 

MaX|MUM 

AMPLITUDE 

,010 

•100 

< *333 

•*00 

1*000 
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1*404  { 

1 • 1*5  I 
1*01*  I 
•*16  I 
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100 
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1000 
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Table  No.  8 

. 

Heave  Response  of  Complex  Shape  Buoy 

RMS  OF  WAVE 

SPECTRUM  • 

8*683  FEET 

probable  amplitude 

OF  WAVE 

a 

1.897  FEET 

Fraction,  jf 

EXPECTED 

LARGEST 

average 

WAVE 

amplitudes 

WAVE 

I NUMBER 

MAXIMUM 

CONSIDERED 

amplitudf 

I OF  WAVE8 

AMPLITUDE 

• 010 

6*330 

! 50 

5*689 

• 100 

4.830 

I 100 

6*118 

• 333 

3*799 

! 500 

7.003 

• 500 

3.370 

I 1000 

7.459 

1.000 

8*377 

I 10000 

8.399 

? 100000 

9.311 

RmS  8f  RESPONSE  SPECTRUM  • 

7.489 

DEG 

PROBABLE  AMPLITUDE 

OF  ROLL  RESPONSE  ■ 

5.295 

DEG 

amplitudes  OF 

fraction  of 
largest 

AMPLITUDES 

CONSIDERED 

ROLL  ARfc  IN  DEGREES 

avehase  I 
ROLL  I 

AMPLITUDE  1 NUMBER 

RESPONSE  I OF  WAVES 

EXPECTED 

roll 

MAXIMUM 

AMPLITUDE 

.010 

17.667  I 

50 

15*877 

.100 

13.481  I 

100 

17*076 

.333 

10*605  I 

500 

19*547 

• 500 

9.407  I 

1000 

80*820 

1*000 

6*636  I 

10000 

83*448 
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100000 
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Table  No.  9 Roll  Response  of  Complex  Shape  Buoy 


Average  Heave 

Significant  Heave 

Average  Wave  Amplitude 

Average  Wave  Amplitude 

1.000 

3.961 

1.  625 

6.591 

0.301 

2.  791 
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5.0  CONCLUSIONS  AND  LIMITATIONS 

The  theoretical  introduction  and  the  case  studies  presented  point 
out  the  positive  aspects  of  the  computer  solution  as  well  as  some  of  its 
limitations.  In  the  formulation  of  the  equations  of  heave  and  roll  motion, 
an  attempt  has  been  made  to  account  for  the  effects  of  the  water  particle 
velocity  and  acceleration.  The  depth  dependence  of  these  effects  has 
been  included.  The  model  thus  obtained  is  more  realistic  than  simpler 
models  which  consider  only  buoy  displacement  and  wave  slope  as  the 
predominant  exciting  forces. 

As  illustrated  in  case  study  number  3,  the  response  of  buoys  of 
relatively  complex  shape  can  be  easily  studied.  The  heave  and  roll 
response  amplitude  operators  can  be  used  to  compute  the  response  of 
the  buoys  to  waves  of  known  or  specified  amplitude  and  frequencies.  In 
addition  to  this  time  domain  approach,  specified  spectral  densities  can 
be  used  to  derive  certain  statistical  expectancies  of  buoy  heave  and  roll 
amplitudes.  Parametric  studies  of  buoy  performance  can  thus  readily 
be  made. 

On  the  other  hand,  to  satisfy  the  condition  of  linearity,  certain 
assumptions  are  made  which  introduce  in  the  solution  a degree  of 
arbitrariness  difficult  to  evaluate.  Certainly  the  initial  choice  of  the 
average  values  of  wave  amplitude,  buoy  heave,  and  buoy  roll  angle 
used  to  compute  the  linearized  coefficients  of  drag  and  inertia  will 
reflect  on  the  accuracy  of  the  solution.  To  improve  this  accuracy  an 
iterative  procedure  can  be  followed  which  replaces  the  initial  assumed 
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values  by  computed  ones  until  sufficient  agreement  is  achieved. 

Experimentally  verified  values  of  linear  equivalents  of  inertia 
and  viscous  effects  would  greatly  help  validate  or  improve  the  com- 
puter solution  described  in  this  report.  The  assumption  of  small 
roll  angles,  also  required  by  the  condition  of  linearity,  further 
limits  the  use  of  this  program. 

Energy  dissipation  by  wave  radiation  is  not  considered.  This 
factor  could  be  important  in  large  disk  buoys.  Finally  the  effects 
that  mooring  lines  and  tether  lines  would  have  on  the  buoy  response 
have  not  been  included,  thus  restricting  this  solution  to  free  floating 
buoys. 

Despite  these  limitations,  the  rationale  used  in  the  derivation 
of  the  solution  and  the  program  input  flexibility  make  the  computer 
solution  useful  as  well  as  practical. 
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7.0  APPENDICES 


APPENDIX  I 


Expression  of  Linearized  Damping  Coefficient 


locity  V i( 


When  the  drag  force  on  a body  moving  with  a velocity  r is  assumed  to  be 
linearly  proportional  to  the  velocity,  the  expression  of  the  force  is  simply 

J>t  = dV 

In  most  cases,  however,  drag  forces  are  expressed  in  terms  of  V , 


using  the  familiar  formula 


JW  fc^r 


where  f>  - water  mass  density 
Op  = drag  coefficient 

~ body  area  across  the  flow  (blunt  bodies). 

If  the  motion  of  the  body  is  periodic,  with  amplitude  /\a  and  frequency 
CU  - §JL. , i.  e.  if  for  example 

ft  — CO  ^ 

then  the  amount  of  energy  dissipated  per  cycle  by  the  damping  force 
is  given  by  'T  ‘ • 

= J\  JC  = fpt  Vdt  - ^ 

° 0 o 

The  amount  of  energy  dissipated  per  cycle  by  the  linear  damping  force  is 


in  turn  given  by  7- 


T 

/ 2. 


t S Jp(  Jt  r Jz>f  Vjt  «,  t/  / CO  todrj 


t = ncL  X 


/V 
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Assuming  both  forces  to  dissipate  the  same  amount  of  energy  will  yield 
the  expression  of  the  linearized  damping  coefficient  o[  , namely 

d J:  f Cj  u> 
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APPENDIX  II 


Evaluation  of  the  Coefficient  "B"  of  Damping  Moment 
"B"  has  been  previously  defined  as: 


Let  us  consider  a buoy  made 
of  different  cylinders  as  shown 
on  the  sketch.  (Fig.  No.  24) 
Let  cj(  , , be  the  values  of 

and  lft  and  be  the 
corresponding  limits  of  the 
variable  i . The  integration 
of  the  first  integral 


s 


r.-o 


Fig.  No.  24 


yields: 


This  result  in  turn  yields  to  the  recurrence  formula 

% *h  - 7- ! r'** 

n*°  J 
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The  expression  of  the  coefficient  "B"  is  therefore 


o/gu 

4 


(pn^) 


where 


+ Y~  /•  ■ 

<L<7  w /*. 


*’6/ 


C*  6> 
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APPENDIX  III 


Evaluation  of  the  Coefficient  "D11  of  Wave  Drag  Moment 
The  coefficient  "D"  has  been  previously  defined  as 

.iJft  K-te  .!Xz, 


t]  )/?e 


!?*o 


- J4y  fi 


<L  o!rl_ 


The  integrals  can  be  readily  evaluated  with  the  help  of  the  following  argu- 
ment. For  small  angles  Q , the  projection  of  A*  on  the  vertical  is 
approximately  equal  to  A*. 


Thus,  from 

A * 0 
/ 

to 

n = J 

/ 

;z 

s.  -7 

and  from 

r{*0 

to 

ut 

H 

/ 

-z 

SV 

Introducing  these 

values  of 

in 

the  integrals 

yields 

-ZfS  . K‘S  if?  - ZKrl. 

X)  -/&*.  ) r;drt  _ Jd(n.) n e «/q.  < 

*Cx0  rjo  ‘ 


Considering  the  same  buoy  geometry  as  in  Appendix  II,  and  evaluating  the 
first  integral  over  the  domain  of  variation  of  A*  yield: 

fa [*  - e (z^r. 

This  result  leads  to  the  recurrence  formula: 

^ , dl 

MW  e dyH,  =4  U A (in.  -/)  ~e  (^i. 


= ^ ~)  ~e  (Z*L  -/)] 


L ~ 
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\0 


Similarly,  the  evaluation  of  the  second  integral  over  the  domain  of  variation 


The  recurrence  formula  thus  is 
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APPENDIX  IV 


Evaluation  of  the  Coefficient  "P"  of  Wave  Inertia  Moment 


'P"  has  been  previously  defined  as 


( rr-o  k*o 


Using  the  change  of  variables  previously  discussed,  the  integrals  in  the 


bracket  can  be  written: 


i t -*M 


? , 4p+n) 

+ M(fi)  >ze-  dfL 


^ + J V '(n^. 7je  <j\  l 

° 0 J 


-ML 


Noting  that 


r **  , , 

-/'e  ^ = fr 


Then,  over  the  intervals 


O*  Kg  A,  M/&  cf(n)  ac/, 

6/$*?  d(>i)  - 


the  evaluation  of  the  integrals  yield: 


1)  over  the  first  interval 


I 
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2)  over  the  second  interval, 


ji  r 

*L  e.  ( 

r*  L 1 


(,  . fat.  7 

- e- 


It  thus  appears  that 


_,fr  _ks 

- fitf)  ^ ~ ft  }j£[±ltyrO  ~ e M 


t?zO 


As  far  as  the  second  integral  is  concerned,  noting 


L */,t 


- e 

~ ** 


and  evaluating  over  the  range  £ ^ 4 1£  ^ (2.  j ftzj  - *h.J 


yield: 


which  shows  that 


4' 


x.  / -K  - M j 

/e  («el+j  - £ 


” ^ I'v  yj"/^  • — /7 y* 

^ ftp  ft  - e (*e/*o] 


r-',^ 
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The  expression  of  the  coefficient  "P"  is  therefore 


* 
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APPENDIX  V 


Evaluation  of  the  Coefficient  "I,,"  of  Added  Moment  of  Inertia 

The  moment  of  inertia  of  the  water  displaced  by  the  buoy  with  respect 
to  the  buoy  c.g.  is  evaluated  with  the  help  of  the  parallel  axis  theorem, 


and  is  given  by 


- \ ^ 


where 


L+  ^ ( J ^ 3 ‘ 


!(*), = 


moment  of  inertia  of  cylinder  "t_"  with  respect 
to  its  own  c.g. 

49  c - £ 

wi  th  Mi,  - mass  of  water  displaced  by  cylinder  "C" 
= radius  of  cylinder  "6" 

' - height  of  cylinder  "C," 


- distance  of  c.g.  of  cylinder  " L"  to  keel 


- distance  of  buoy  c.  g.  to  keel. 
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APPENDIX  VI 


Computation  Method  for  Coefficients  1 * B1 1 # "D"  and  "P" 

The  actual  computation  of  the  roll  response  amplitude  operator  in  com- 
puter program  ROLLRAC  is  performed  using  different  forms  of  the 
expressions  for  gome  of  the  moments.  This  was  done  because  the  buoy 
configuration  is  input  as  geometri cal  "solids"  rather  than  surfaces.  It 
was,  therefore,  more  straightforward  to  implement  the  computation  of 
the  damping  forces  using  an  iterative  procedure  on  the  components.  The 
following  expressions  give  the  form  of  the  equations  used. 

For  the  buoy  damping  moment,  "B" 

cs*  l 


where  $(  - 


ei 


3k. 


as  defined  in  Appendix  II 


dt diameter  of  the  ith  buoy  component 
= depth  to  the  buoy  center  of  gravity 

depth  to  the  bottom  surface  of  the  ith  buoy  component 
Zn  - depth  to  the  top  surface  of  the  ith  buoy  component 

, . ^ z-~Zcfh 


4i  - number  of  buoy  components 


For  the  water  damping  moment,  "D" 


D 

4X3- 


/*)* 


1 -r  /**< 

L°1 

J Z^rc 


T 
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*P 


where 


and 


/3  a cj  4s. 7 a 

-3rt- 


as  defined  in  Appendix  III. 


For  the  water  particle  acceleration  inertia  moment  "P> 


JO* [(-*(*-*$)- M)e- 


P-.Xi. 
/r^ 


where  ^1.  as  defined  in  Appendix  IV. 


I 


i 


+ e 
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APPENDIX  VII 

Heave  Program  Listing 


1. 

C 

PROGRAM  H£RA0 

2. 

C 

3. 

c 

VERSION  2*0  JAN.  1977  R.  G8LDSMITH 

6. 

c 

VERSION  1.1  JUNE.  1976  R«  SBLDSMfTH 

5. 

c 

6 # 

c 

THIS  PROGRAM  IS  USED  19  COMPUTE  THE  HF A VE  RESPONSE 

7, 

c 

amplitude  operator.  and  associated  phase  angles,  pop 

8 • 

c 

SPAR  TYPE  dUOY  SYSTEMS 

9. 

c 

10. 

c 

VERSION  2 • 0 • MODIFIED  to  incorporate  a wave  drag  cbeff 

11. 

c 

*• 

12. 

c 

13. 

L0GICAL  IAMTERM 

1*. 

c 

15. 

dimension  DEPTHP(?s),AREAI251 

16. 

dimension  CFPTHI (25). aDDm5C( 251 .volume  125) 

17. 

dimension  DERTHO(2b)»DAMPC(2b).WDRAGC(25) 

1*. 

DIMENSION  maXv.aVN8(6)  . WVMAXCPF(a  ).HEVMAXHT(6) 

19. 

DIMENSION  FRaCamPS(5).AVRC0EFF(5). AVRESpNS(5) 

20. 

c 

?1. 

DATA  NCR. NLP/10’5. 108/ 

22. 

DATA  NMAX/25/ 

23. 

DAT*  Hl/3. 16159?/. RTPU/57.2958/ 

2*. 

DATA  RH8. a/1 *990351  32*176/ 

25. 

c 

26. 

DATA  FRACAMPS  /0.01a0.10i0.333iO»50a1.0  / 

?7. 

DATA  AVRCOEFF  /P.359.1 »8o0a1 *616. 1.256.0*886  / 

28. 

DATA  MAXRAVN8  /SO* 100.500* 1000# 100 00 i l 00000  / 

?9. 

DATA  WVMAXC9F  / 2. 12.2.2". 2. 61. 2.78. 3. 13*3. 67  / 

30* 

c 

31. 

c 

32. 

c 

INITIALISATION 

33. 

c 

36. 

NP  • 0 

35. 

M • a 

36. 

ND  ■ 0 

37. 

CARfcA^L  • 0*0 

38. 

VIRTMASS  ■ 0.0 

19. 

ISEASEL  ■ o 

*0« 

TIMtl  • 0.200 

61. 

TIM£2  ■ 50,0 

62. 

TImedel  ■ 0.200 

63. 

wtNUV  ■ 0.0 

66. 

R 1 2 • P I *?• 0 

65. 

RMOlj  • RH0*0 

66. 

c 

67. 

c 

CHECK  Fr,R  ON-LINE 

68. 

10NFLAQ  ■ 0 

69. 

IF  UAMTERM  (10UM)  ) IONFLAO  • 1 

50. 

c 

51. 

c 

input  data 

52. 

c 

53. 

o 

o 

CONTI NUF 

56. 

wRJTt  ( NLP*  9600 ) 

55. 

c 

56. 

c 

INPUT  NUMBER  of  PRESSURE  SURFACES 

57. 

IF  (10NFLAQ  .EQ.  1)  WRITE  (NLP.9M0) 

58. 

INPUT  NTEST 

59. 

IF  (NTEST  .GT.  NMAx)  WRITE  (NLP/9700)  NTES7#NMA*  $ STOP  100 

r — 

1 

00 

1 

? 

\ 

' 

60* 

IF  (NTEST  .LT.  0)  39  T»  173 

6 1 • 

J F (NTEST  < (i£  < 0)  NP  • NTEST 

63. 

1 F (np  .EQ.  0)  G8  T9  200 

63  . 

C 

INPUT  PRESSURE  TERMS 

64. 

I F UONFLAG  »F.G»  1>  WNJTE  (NLP.94?0> 

68. 

00  ISO  1 ■ l.NP 

66  . 

INPUT  DFPTHP(  I ). AR£A(l  ) 

67. 

150 

CONTINUE 

68  . 

c 

OUTPUT  TERMS 

69. 

1'5 

CONTINUE 

70. 

IF  UONFLAG  «EQ.  11  GO  T0  200 

71  . 

*f?ITE  ( NLP . 9430 ) NP 

72. 

IF  (NP  • QT . 0)  WHITE  (NLP. 9440)  ( DEPTHP ( I ) > ARE A ( I ) . 1 • 1 # NH ) 

73. 

c 

74. 

c 

INPUT  NUMBER  OF  INEHTIAL  COMPONENTS 

75. 

O 

o 

nj 

CONTINUE 

76. 

IF  UONFLAG  «EQ.  1)  WHITE  ( NLP/ 9450 1 

77. 

INPUT  NTEST 

78. 

IF  (NTEST  .GT.  NMAX)  WHITE  (NLP. 9700)  NTEST. NMAx  1 STOP  200 

79. 

IF  (NTEST  .LT.  n)  GO  TO  275 

HQ. 

IF  (NTEST  .GE*  0)  N I • NTEST 

HI. 

IF  (NJ  .EQ,  0)  GO  TO  300 

M2. 

c 

INPUT  INEHTIAL  TEHMS 

83* 

IF  UONFLAG  .EQ.  i>  WRJTE  (NLP.9460> 

84. 

08  2S0  I ■ 1 . N I 

i 

85. 

INPUT  UEpTH I ( I). AOOMSt ( I). VOLUME! I) 

86. 

250 

CONTINUE 

87. 

c 

OUTPUT  TERMS 

88. 

275 

CONTINUE 

89. 

IF  UONFLAG  .EQ.  1)  30  TO  300 

90. 

wRITt  (NLP. 9470)  *1 

91. 

IF  ( N I • GT • 0)  WHITE  (NLP. 9480)  ( DEPT* I ( I )/ ADDMSC ( I) . VOLUME  (!) . 

92* 

< 

► I ■ 1 . N I ) 

93. 

c 

94. 

c 

INPUT  NUMdEW  OF  DRAG  SUHFACES 

95. 

300 

CONTINUE 

96. 

IF  UONFLAG  .EQ.  1)  WHITE  (NLP.9490> 

97. 

INPUT  NTEST 

98. 

IF  (NTEST  ,ar.  NMAX)  WHITE  (NLP. 9700)  NTEST. NMAX  j STOP  300 

99c 

IF  (NTEST  .LT*  0)  00  TO  375 

100- 

IF  (NTEST  *G£ » 0)  NO  • NTEST 

101  • 

IF  (NO  .EQ.  Q)  GO  TO  400 

102. 

c 

INPUT  0«A G COMPONENTS 

► 

103* 

if  UONFLAG  »Eq.  1)  WHITE  (NLP. 9500) 

[ 

104. 

c 

SUM  DRAG  SURFACE  COEFF  F(DEPTN»0) 

10b. 

SUMUCO  • 0.0 

106. 

00  300  I • l.NO 

107. 

INPUT  OFPTHO( IJ.OAMPCI I ) . WOHAGC ( J) 

los. 

SUMJCo  ■ SUMOCO  ♦ DAMPCU) 

109. 

350 

CONTINUE 

110. 

c 

OUTPUT  TERMS 

111. 

373 

continue 

♦ 

112. 

IF  UONFLAG  .EQ.  1>  30  TO  400 

113* 

WHITE  (NLP. 9510)  NO 

114. 

IF  (NO  ,GT,  0)  WHITE  (NLP. 9520)  ( DEPTWD  { I ) . 0 AMPC  ( I ) . WDRAQC  ( I ) t 

115. 

« 

I • l.ND) 

116. 

c 

117. 

c 

INPUT  WATER  LEVEL  CROSS  SECTION  AREA 

1 1 a . 

c 

(TO  SIMPLIFY  THE  COMPUTATION  THIS  IS  ASSUMED 

! 

119. 

c 

constant  over  the  Range  of  vertical  motion  at 

85 


i 


120 

121 

1?2 

123 

1?4 

125 

126 
1?7 
128 

129 

130 
111 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 
1 4* 

145 

146 

147 

148 

149 
1 5q 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 
167 
16* 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 
179, 


C THE  WATER 

400  C8NT1NUE 

I F UCJNR.AQ  .EQ.  1)  WRITE  (NUP,9530> 

INPUT  CSATE8T 

If  (CSATEST  .UT*  0*01  U8  TO  450 
CAREAwt  • CSATEST 
450  CONI INUE 

IF  (18NFLAQ  «N£*  1»  WRITE  (NlP*9540)  CAREAWL 
HP  ■ CAREAW1.4RH83 
C 

C INPUT  VIRTUAL.  MASS 

IF  U8NFLAQ  .EQ.  D WRITE  (NUPi9550> 

INPUT  VMTEST 

IF  (VMTEST  *UT.  0»OJ  0»  T»  550 
V I RT^ASS  • VMTEST 
550  CONTINUE 

IF  (18NFUAQ  .N£„  1}  WRITE  INUP/9560)  VlRfMASS 
C 

C INPUT  TIME  RANGE 

IF  (18NFUAG  .EQ*  1 > WRITE  (NLP*9570> 

RpaU  ( NCR,  9025  ) Tl,T2,T3 
IF  ( T1  .UT.  0*0)  G»  T8  675 
TJMtl  ■ T 1 
T ! M£2  ■ T2 

IF  (TIME?  .UT.  TIME!)  TIME2  • TIME1 

IF  <T3  .UE*  0*0)  T3  4 TIMEOEL 
TIMtUEU  » T3 
675  CONTINUE 

IF  (I8NFLAG  .N£,  1)  WRITE  (NUP,9580>  T I M£i  , T I ME2,  T I MEDEU 
C 

C INPUT  W J NO  VEU8CITY  FBR  SEA  STATE 

700  CONTINUE 

IP  UBNFUAB  »EQ.  1»  WHITE  (N|_P#9590> 
we  A U < NCR#  9020 ) ISTEST* wAVEHT/ WAVEPER 
IF  (ISTE3T  .GT.  3»  JSTEST  ■ *1 
IF  USTEST  .UT.  0)  G»  TB  773 
ISEASEU  ■ ISTEST 

IF  USEASEU  *E0.  II  WINQV  ■ waVEHT 
WJNUVP4  ■ WIN0V«*4 
wAVEMTP?  • WAVEHT*WAVEHT 
*AVHfc.HP4  ■ WA  VEbER#*4 
775  CONI INUE 


IF 

( IBNFUAG 

• EQ. 

1 > 

U6  TB 

800 

IF 

( iseasfl 

• EQ. 

01 

WRITE 

(NUP/9600) 

IF 

( iseaseu 

• EO. 

1 ) 

WHI  TE 

( NUP# 9601 ) 

WINOV 

IF 

( ISEASFU 

• EQ. 

2) 

WRITE 

(NUP.960P) 

waveht, wavepe* 

IF 

( ISEA5EU 

* Eq  * 

3) 

WRITE 

(N|.p*9603) 

wavehT, WAVEPER 

C 

C CHECK  OUTPUT  BPTIBNS 

80C  CONTINUE 

IF  IJONFUAG  .EQ*  1)  WRITE  (NUP#9605) 

H£AU  ( NCR*  9900  ) JU/ IP#RMIN#HHAX 
I U I S»T  • 0 

IF  (IU  «EQ.  1 hy ) IUIST  ■ 1 
IPUOT  • 0 

IF  (IP  .EQ.  1HY)  IPU8T  •! 

IF  UPUBT  .EO.  1»  C AUU  PU8TIN1T 
r 

c C9*PUTE  RESPONSE  AND  PHASE  COMPONENT  BN  TJM£  ITERATION 
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180. 
1*1. 
1*2. 
183. 
1 *4. 

1*5. 

1*6. 

1*7. 

1*8. 

189. 

190. 
1 9 1 . 

192. 

193. 

194. 

195. 

196. 

197. 

198. 

199. 

200. 
201  . 
202* 
203* 

204. 

205. 

206. 
2«7. 
20*. 
209. 
21C  • 
211. 
212. 
213. 
21*  • 

215. 

216. 

217. 
21*. 
219. 
??0« 
221. 
222. 

223. 

224. 

225. 
2?6. 

227. 

228. 

229. 

230. 
211. 
232. 
213. 
23*. 
215. 

236. 

237, 

218. 
239. 


900  CONTINUE 

IF  ( 1U  1ST  »EU.  1)  WRITE  (NLP.9610) 

9WS1NTG  • 0*0 
SINTU  ■ 0. 

09  2000  TIME  ■ T I *£1,  TIME2#  TIMFDEL 
FKtU  • PI?*10000000000«0 
I*  ( T I ME  .NE.  0.0)  PKEQ  ■ PI2/TlME 
FheuP*  ■ FPEQ»*4 
FKfcUPS  ■ FRE0*FREUP4 
E*PT£RM  ■ -FRfQ«FPeU/U 
C 

C StJM  PRESSURE  COMPONENTS 

SUMP  • Q.o 

REPEAT  1100.  P(5W  l ■ (l.NP) 

SUMP  • SUMP  4 HH*a**9EA( I )*EXP(EXPTEWM*DEPTHP( I) ) 

1100  CONTINUE 

C SUM  INEHTIAL  COMPONENTS 

SUMJ  . o.O 

PE  PEAT  1200.  FOR  I • tl.M) 

SUM  I « SUPJ  ♦ AOOMSC  ( I ) »V8LUPE  U > *Ex|a<E*pTERN,*DEPTMl  II  ) ) 

1200  CONTINUE 

SUNJ  • RW0»FPEU«EPEU*SUMI 

c 

C SUM  ORAQ  COMPONENTS 

SUMO  • 0.0 

PEPEAT  1300.  F8H  I ■ ( 1 # NO ) 

SUMO  a SUMO  ♦ wOPAGCm*EXp(EXPTEWM*OEPTHOm  ) 

1300  CONTINUE 

SUMO  • f»EU»FPEU*SUMO 


SUMOC  ■ FHEU*FWE*W*SUMDCO 
C 

C Commute  HESPBNSE  amplitude  OPERATOR 

RNUM  « I SUMP  . SUMJ)**p  ♦ SUMO**g 

RUEN6M  ■ <WF  - VIRtMasS#ERE0*FPeQ)*«2  ♦ SUMQC**2 

RAO  a SOBT ( RNUM/K0EN6M ) 

c 

C COMPUTE  phase  BETWEEN  FORCE  AND  REAVE 

pMJ  a RTB0*ATAN2(. SUMOC.  CPF  . VJPT*vASS*FRE0#FREC3)  ) 
C 

C COMPUTE  PHASE  BETWEEN  FORCE  AND  WAyE 

SJUMA  a RT9D«ATANgCSUMD. (SUMP  • SUMI)  ) 

C 

c COMPUTE  PHASE  BETWEEN  wave  and  HEAVE 

THETA  ■ SIGMA  ♦ PHI 

c 

c GET  sea  SpECTPA 

call  SEASPEC 
C 

C compute  response  and  integrate 

HKS  a RA8*PA8*S 

I*  (TIME  .LE.  T I M£ 1 ) qo  TO  1400 

OELF  a rPE  UL AST  » FH£U 

HRSINTG  • ( HRS  ♦ HPSLAST  > -*0  » 50*DELp  ♦ RHSINTQ 
SJNTG  • (S  ♦ SLAST)*0.^0*GELF  ♦ SINTG 
1400  CONTINUE 

RPSLAST  a R«S 
SLAST  • S 
FHEULAST  • FRE0 

c 

C 


8UTPUT  LIST  IF  IT  WAS  SELECTED 


8? 


# 


2*0. 

2*1. 

2*2. 

C 

2*3. 

c 

2**. 
2*5. 
2*6 . 

1500 

2*7. 

2000 

2«8. 

c 

2*9. 

c 

250. 

c 

251. 

252. 

253. 

c 

25*. 

255. 

256. 

c 

297. 

2300 

258. 

259. 
2*0  • 

C 

261 . 
262* 
263. 
26*. 

2*00 

265. 

♦ 

266. 

A 

267. 

C 

268. 

269. 

c 

1 

270. 

271. 

272. 

273. 
?7*. 
275. 
275. 
277. 
27*. 
279. 

290. 

29 1. 
2*2. 
283. 
2*“  • 
285. 
2**. 
2*7. 
28*. 
289. 

2 90. 

291. 

292. 

293. 
29*. 

295. 

296. 

297. 

298. 

299. 


C 

c 

c 


!►  < IUIS7  ,l£.  01  08  fB  1500 

**JTE  (NL°»N615>  TIME.freL.WaB. SIGH*. PHI. THETA, $ 
Cm£Ck  FUR  PLBT 

CB^TfNUF 

l¥  < lpL8T  »LE  • 01  38  TB  2000 

C*LU  **U*T3  ( * • • »RAB,  TIME.  1 > 


net  statistics 

- COMMUTE  Wtm  hEAM  SQUARE  BF  WAVF 
R“S  • S'JRT  ( SI  NTG  ) 

MR*  B AMp  a 0 * 707*RMS 

COMMUTE  AVERAGE  WAVE  HEIGHT 
o:i  2300  I a 1,5 

AVRESMNSm  • AVRC8EFFU  )»HMS 
CONTINUE 

n„  , C8HMUTE  maximum  wave  AMPLITUDES 

L>3  c*00  I ■ 1«6 

HtVM*XMT(I)  a wVMAXC0F( I 1*RMS 
CONTINUE 

I T t (NLP, 96351  RMS 
"HI  rt  (NLP, 96*0)  RH88AMP 

"RITE  {NLP.  96*51  (EHACAMPSdl^VRfRPNSm, 

MAXWAVNBll  )«  ME V“AXMT ( I 1 j 1,1,5), 
MaXWAVNt)  (61.HLVMAXHT16) 

COMPUTE  ROOT  MEAN  SQUARE  8*  RESPONSE 
HMS  • SORT ( HRS  INTO)  6 

MRUOAMP  a 0.707«RMS 

C COMPUTE  AVERAGE  RESP8NSE  8F  WEAVE 

U8  2500  I • 1*5 

AV*ESPN<5(I)  a AVKCBEFF  { n ARMS 

2500  CONTINUE 

C COMPUTE  MAXIMUM  AMPLITUDES  8F  HEAVE 

D8  2600  I * 1.6 

HEVMAXHT(I)  a WVMAXC8F{ I)»RMS 
2600  CONTINUE 

"HITt  (NLP. 9620)  RM5 
WRITE  (NLP. 9625)  PROUAMP 

"RITE  (Nt». 96301  (FHAC Amps ( I ) . AVRESPNSt I ). 

* MAXWAVNB ( n .HEVMAXMT { I ),  I a 1 4 5 | . 

* MAXwAVNB(6».heVmaxht(6) 


check  for  plot 

IF  (IHL8T  .LE»  01  GO  T8  3000 
Call  mlotweav 
3000  CONTINUE 

IF  (IBNFLAQ  *EQ.  II  WRITE  (NLP. 9655) 
RfaU  (NCR, 901b. ENO«8COO)  J£ND 
IF  (lENO  . EU * 1**1  QB  TO  100 


8000  CONTINUE 
ST8M 

m 

» 

9000  FORMAT  < *1.1*. AJ.lX.gC.Q) 
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300.  9015  FORMAT  { *1 ) 

301.  9o?o  format  (I.f.oaf.oi 

302.  9025  FORMAT  (JP.O) 

303.  C 

30*.  9*00  FORMAT  <1W1#*  W£AV£  RESPONSE  AN*|_VSIS  PROGRAM'/ 

305.  ♦ ' ALL  DEPTHS  ARE  POSITIVE'/) 

306.  9410  FORMAT  I/'  Input  ThE  NUMBER  OF  PRESSURE  SURFACES  9') 

307.  9420  F (IRWAJ  </i  ENTER  F OP  EACH  SURFACE'/ 

308.  ♦ ' OK pth  (Ft)#  area  (SO  FT)  (•  ARp A F8R  DOWNWARD  FORCE) ') 

309.  9430  FfiqMAT  (/•  NUWSrP  OF  PRESSURE  SURFACES 

310.  9440  FORMAT  (/•  DEPTH  * ► T » * A P£  A ( So  FT)  (•  AREA  FOR  DOWNWARD  FORCE)'/ 

311.  ♦ <Flo.3.2X#FM.?>  ) 

312.  9450  format  (/'  INPUT  THE  NUMBER  OF  INERTIAL  COMPONENTS  O') 

313.  9460  FORMAT  </'  ENTER  K OP  EACH  CBMPRS.ENT'/ 

314.  ♦ « DEPTH  (FT)#  ADDED  MASS  CBEFF#  VOLUME  (CU  FT)*  ) 

315.  94/0  FORMAT  (/'  NUM8FP  OF  INERTIAL  COMPONENTS  *'«I5) 

316.  9480  FORMAT  (/«  DEPTH  (FT)  ADDED  MASS  CBEFF  vOLUmE  (CU  FT)'/ 

31 7.  ♦ (Fl0*3*6XiFl0*3.7XiF12,3)  ) 

318.  9490  FORMAT  (/'  INPUT  TH£  NUMBER  OF  DRAG  SURFACES  9*> 

319.  9500  FORMAT  (/'  ENTER  FOP  EACH  SURFACE'/ 

320.  * ' DEPTH  (FT)#  DAMPING  CR£FF#  WAVE  DRAG  COEFF' 

321.  ♦ • (LHF/(FT/SEC)/(WAD/SFC)  ) ' ) 

322.  9510  FORMAT  (/«  NUMBER  OF  DRAG  SURFACES  »'#I5) 

323-  9520  FORMAT  (/'  DEPTH  (FT)  DAMPING  COpFF  WAVE  DRAG  COEFT' 

324.  ♦ < O.BF/(FT/SEC>/(RAD/SEC) ) * / 

3?5.  ♦ (FI  0*  3#  3X#F1 1*3*  5X#Fli.3)  ) 

3?6 . 9530  format  (/•  ENTER  CROSS  SECTION  Are*  *T  SURFACE  (SQ  FT)  O') 

327.  9540  format  (/•  CROSS  SECTIONAL  AREA  AT  SURFACE  ••#FiO»4#'  SO  FEET') 

328.  9550  FORMAT  (/'  ENTER  VIRTUAL  MASS  (SLUGS)  O') 

3?9.  9560  format  (/i  VlRTu*L  MASS  ■ '#  F11.2,  ' SLUGS') 

330*  9570  FORMAT  {/•  ENT FR  START#  END#  INCREMENT  OF  PERIOD  RANGE  (SEC)  9«) 

331.  9580  format  (/•  period  range,  in  seconds  start  end  delta*/ 

332.  ♦ 26X#F8»3#F8.3#2X#F8.3) 

333.  9590  format  ( / » ENTER  SEA  SPECTRUM  type  AND  PARAMETERS' 

334.  ♦ /»  • 1.0  0 ' 

335.  ♦ /'  PlE*S0N-M9SK»wITZ  l,  WIND  SPEED  (KNOTS)  • 

336.  ♦ /'  PRETSCHNEIDE*  2#  SIGNJF  WAVE  HT  (FT)#  SIGNJFf 

337.  ♦ ' WAVE  PERIOD  (SEC)* 

338.  ♦ /'  I.S.S.C*  3#  SIGN  IF  WAVE  HT  (FT)#  SIUnIF* 

339.  ♦ ' WAVE  PERIOD  (SEC)'  ) 

340.  9600  format  (/'  SEA  SPECTRUM  • 1,0') 

341.  9601  format  (/'  PIERSON. mbs<*wJ TZ  SE*  SPECTRUM  « 

342.  ♦ /•  WIND  SPEED  • '#Fio.3#'  KNOTS') 

343.  960?  format  (/•  bretschnejder  sea  spectrum 

344.  ♦ /'  significant  wave  wt  ■ *#F10»3#'  FEET' 

345.  ♦ /'  SIGNIFICANT  wave  PERIOD  ■ '#Fto#3#'  SEC') 

346.  9603  format  (/•  I.S.S.C*  SEA  SPECTruM* 

347.  4 /'  significant  wave  mT  ■ * #Flo»  3#  * FEET' 

348.  ♦ /'  SIGNIFICANT  wave  PFRIOD  •>  *,Fiot3#*  SEC) 

349.  9605  Format  (/«  ENTER  TORN  FOR  A LISTING#  PLOT  OF  RESPONSE  »/ 

350*  ♦ * FOR  PL"T#  you  may  ALSO  ENTER  RAO  m[n,max  q « ) 

351.  9610  FORMAT  ( 1 H 1 / ' PERIOD  AnG  FREW  RAO  W-F  PHASE  F-H  PHASE  • 

352.  4 ' w.w  PHASE  amp  SPEC  ') 

353.  9615  FORMAT  (Fjo.3#Eic*3'fIO'3*Fio.3#fi0.3#FI0.3»Fio»3) 

354.  9620  format  (///////•  RMs  OF  RESPONSE  SPECTRUM  •«#F1o.3#>  FEET') 

355.  9625  FORMAT  (/*  PR9BA8LE  AMPLITUDE  '/ 

356.  ♦ • OF  heavE  RESPONSE  »'#Fi0»3#'  FEET') 

357.  9630  FORMAT  (// 

358.  ♦ /'  FRACTION  OF  AVERAGE  ! EXPECTED  • 

359.  ♦ /'  LARGEST  WfcAVE  I WEAVE  ' 
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360. 

361. 
362  • 

363. 

364. 
36b. 

366. 

367. 

368. 

369. 

370. 

371. 

372. 

373. 

374. 

375. 

376. 

377. 

378. 

379. 
3®0. 
381. 
3*2. 

, 383. 

\ 384. 

385. 

386. 

387. 

388. 

389. 

390. 
391  . 

392. 

393. 

394. 

395. 

396. 

397. 

398. 

399. 

400. 

401. 

402. 

403. 

» 404. 

405  » 

406. 

407. 

408. 

409. 

410. 

411. 

412. 

413. 

414. 

415. 
416* 
417. 
411. 
419. 


♦ /'  AMPLITUDES  AMPLITUDE  I NUMBER  MAXtMUM 

4 /*  CONSIDERED  RESPONSE  I OF  WAVES  AMPLITUDE 

♦ I 

♦ 5</T4iF5«3»  T 1 6#  F 9 • ?.  TgS.'I*.  T33»l6.  T43.F9.3) 

4 / T?H»  t I • # T33.16.  T43.F9.3  //) 

9635  Format  (lMl////*  wms  of  WA Vfc  SPECTRUM  •SF10.3.*  FEET') 
9640  FORMAT  (/•  PROBABLE  AMPLITUDE  •/ 

♦ ' OF  WAVE  f »*Fl0»3# ' FEET*  ) 


9645  FORMAT  (// 

♦ /'  FRACTION  OF 

♦ /'  largest 

♦ /'  AMPLITUDES 

♦ /'  considered 

♦ 


I EXPECTED 

AVERAGE  I WAVE 

WAVF  I NUMBER  MAXIMUM 
AMPLITUDE  I OF  WAVES  AMPLITUDE 
I 


♦ b(/T4,F5.3,  T 1 6. A'9 • 3 . T2«,*J',  T33.I6,  T43.F9.3) 

♦ / T28.'H,  T33»  16.  T43iF9.3  I 

9655  FORMAT  i/i  00  YOU  WANT  ANOTHER  CASE  8*) 


t 

I 


C 

9700  FORMAT  ('  NUMBER  OF  ENTRIES  IS  GREATER  THAN  ARRAY  SIZE  ALLOWS' 

4 /14X. I8i?0X. 17. 

4 /«  THE  PROGRAM  TERMINATES'  J 

C 

c * 

c 

SUBROUTINE  SEASPEC 

IF  (ISEASEL  .EQ.  0>  S * 1»0  I RETURN 
GO  TO  (4tco. 4200. 4300).  ISEASEL 
C PIERSON  - MOSKOWITZ 

4100  CONTINUE 

S ■ 135.0/FRE0P5*EKP«  •97ooO'0/(FREf3P4*wINDVP4)  ) 

C CORRECT  For  DOUBLE  height  SPECTRUM 

S ■ s/8.0 
return 


C BRETSCHNEIDEH 

4200  continue 

s • 4203.c#^AVEHTP2/(wAVPERP4*FREQP5>* 

♦ EXB<-105O.0/<WAVHcRP4*FRFUP4) ) 

C CORRECT  FOR  DOUBLE  HEIGHT  SPECTRUM 

S • S/8.0 
return 

C I.S.S.C. 

4300  CONTINUE 

S ■ 2760.O44AVEHTP2/(wAvPERP4»FREQP5)* 


4 ExP('630*C7<waVPEPP4*FRE0P4) ) 

C CORRECT  FOR  DOUBLE  HEIGHT  SPECTRUM 

S ■ S/8«0 
RETURN 
C 
C 

C * * 

c 

SUBROUTINE  PLOTINJT 

C THIS  SUBROUTINE  IS  USED  T9  INITIALIZE  A LINE  PRINTER 

C PLBT  OF  THE  HEAVE  RESPONSE. 

C 

DIMENSION  JPLOTBUF ( 3300) 


C 

DATA  JbUFSIZE  /3300/ 


C 

C 


90 


*4 


430. 

431 . 

432. 

433. 

434. 
435  . 

436. 

437. 

438. 

439. 
430* 

431. 

432. 

433. 

434. 

435. 

436. 

437. 

438. 

439. 

440  * 

441  . 

442. 

443. 

444. 

445. 

446. 

447. 

448  . 

449. 

450. 
451  • 

452. 

453. 

454. 

455. 

456. 

457. 


IF  (NMAX  .or.  RMJN)  ae  T8  1490 
KM AX  • 100*RMIN 

IF  (NMAX  ,£U.  0.0)  HMAX  4 3.0 
1490  CONTINUE 

WOEL  • RRAX  • RMIN 

HMAX  • IFlX(RDEl./9*0  ♦ 0.999)  »5 .0  ♦ RMJN 

C 

NUNES  • <T!"E2  * TIMED/TIHeoei. 

NMAH8  ■ (NUNES  * 9.11/10*0 

TMAK  ■ TIM£v  . NBANS*10*0»TI^EOEL 

N'LINES  * (TMAX  . T I^El  )/T  IMgOEL  ♦ 1*0 

IF  (NUNES  •(.!•  IBUFsUE/13)  00  TO  5000 

JRLOT  ■ 0 

*R!Tt  (NLP.9710) 

RETURN 

C 

5000  CONTINUE 

CAUL  HLOTl  (NBARS* 10*5#10> 

CALL  PLOTS  ( jPLOTbUF.NMlN/RMAX.TMAX,  TIMED 
RETURN 
C 

9710  F0«"AT  (/•  THE  PLOT  BUFFER  IS  NOT  LARGE  ENOUGH  I 

♦ /'  FOR  the  period  range  specified  » 

♦ /'  THE  PLOT  IS  SUPPRESSED’) 

c 

c 

c ******* * * 

c 

SUBROUTINE  PLOTHEAV 

c 

C THIS  SUBROUTINE  IS  USED  TO  OUTPUT  THE  LINE  PRINTER 

C PLOT  OF  THE  HEAVE  RESPONSE. 

C 

CALL  BLOT  5 O#  33* ' WE *VE  RESPONSE  AMPLITUDE  OPERATOR') 

CALL  PLOT 4 <14,  ' PERIOD  (SEC)  ') 

CALL  PL9T7  (10) 

RETURN 

END 


f 
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APPENDIX  VIII 


t 


1. 

2. 

3. 

4. 

5. 
6 » 

7. 

8. 

9. 

10. 

11* 

12. 

13* 

14. 

15. 
15. 

17. 

18. 
19. 
SO. 
21. 
22. 

23. 

24. 

25. 

26. 

27. 

28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 

41. 
42  • 

43. 

44. 

45. 
4b. 

47. 

48. 

49. 

50. 

51. 

52. 

53. 

54. 

55. 

56. 

57. 

58. 

59. 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 


Roll  Program  Listing 

PRBgRAM  RBLLPAO 

VERSJ8M  1*0  SEP,  1976  R«  GBLOSMUH 

this  PROGRAM  is  USED  TB  COMPUTE  THE  RBLL  RpSPBNSE 

AMPLITUDE  OPERATOR,  AND  ASS0C I ATED  PHASE  ANGLES,  FOR 
SPAR  T THE  BUSY  SYSTE^S, 

CURRENT  VERSION  RESTRICTS  DESIGN  TO  CYLINDRICAL  AND 
TRIANGULAR  BBDIES  BN  END,  AND  RECTANGULAR  PLATES, 


LOGICAL  lAMTERM 

DIMLNSIBN  FRACAMRS  lg>,  A VRCOEFF I 5 ) , AVRESPNSI*) 

0 1 Mi.NS I f)N  HAXWAVN3<6),  *VMaxC0F(6),  R»LLMaXI6) 

CBMMhn  / IPD£V  / NCR, NLP 

COMMON  / TP  / T I M£ l , T I ME2, T I MeDFL, PR£0, WA VEN 
common  / SEASTATE  / ISEASEl< w INDV, wAVEHT, WAVEPER# 

♦ RINDVM4.WAVEHTP2#WAVPERP4 

CBMMPN  / BINS  / NPARTS, ISHAPE(50),WtDTH(50>.HEJGHT(50),THlCK(S0), 

♦ DfcNSITv{R0)iDISTcGK(50).FRArN8RM(50> 

COMMON  / BOUTS  / VBLUmE;50),height(50j 

CBMdON  / WATERDIS  * HO(50J,HO(50>,XD(5O>,V0(5O),F0(50), 

♦ DEPTwb  C 00 ) , DE-PTHT  ( 50 ) 

COMMON  / CONSTANT  / PI.RTBD.RHB.G 

CBMMON  / BUOY  / NPMAX,R4L.AYERQ*Mpf TWETABAR,PERieOo# 

♦ 0EPTHK,BU8YCUK,DEPTHCG,BUBYCBt<,0EPTHCB,H0ISPLAC 
CBmmBN  / CBEFS  / 0RA0I5) »c»EFM(5> 

COMMON  / MOMENTS  / BUMTMI, ADOMI, VIRTINRT,WATERIM,BU0YMR, 

■»  buoymdt,«uoymd,watermd,dampm 

common  / outputs  / ilist, iplbt,rmin,Rmax 


DATA  NCRfNLP  /105.108/ 

DATA  P I , RTBO  /3. 141592,57.2958/ 

DAT*  RMO,r,  / 1.99035,  32.174/ 

MAXIMUM  ARRAY  SUES 
DATA  NPMAX  /50/ 

DRAG  COEFFICIENTS  FBR  CYLINDER  AND  PLATE 
DATA  ORAtj  / 1.2,  1*2,  1.5,  0*0,  0*0  / 

ADDED  MASS  CBEt-FjcIENTS  FOR  CYLINDER  AND  PLATE 
DATA  CBFFM  / l.<3,  1.0*  1»Q,  0*0,  0.0  / 

statistical  cbefficients 

DAT*  fRACAmpS  /0«01, 0*10, 0*3333, 0*50, 1.0/ 

DATA  AVRCOEFF  /?. 359, 1*8 00,1*416, 1*256, 0*886/ 

DATA  MAXWAVNB  /SO* 1 00 * 50Q, 1 000, 1 0000, 1 00000/ 
data  wVmAXCBF  /?*12, 2* 28*2*61, 2. 78, 3*13, 3*4 7/ 


INITIalUATIBN 

HI2  • PI*2.0 
RHOG  • RHQ«Q 
TIMfcl  • 0.2 
TIMEZ  • 5q»0 

TIMEUEL  • 0.200 


92 


1SEASEL  ■ 0 
w l NOV  ■ 0*0 
W AVtHT  a 0»0 
wAvtMER  • 0*0 

W -L  a 0*0 

avE*GAMR  a 3*0 

THE  f ABAW  a 11.5/WT80 

NPAHTS  a 0 

CONTINUE 

CHECK  FOR  ON  LINE  AND  INPUT  MODE 
10NFUAG  • 0 
CONTINUE 

I?1  t IAKTEP*  (?DUN)  ) IONFLAG  a 1 i CALL  T INPUT 

IF  UONFtAG  »EQ*  0 ) CAUL  SINPUT 

IF  UPL8T  ,E U.  1)  CALL  PL8TINJT  - - 

CenPUTF.  TOTAL  BUOY  WEIGHT,  DISTANCE  FROM  KEEL  T9 
’center  of  gravity,  depth  OF  cg 

BUOYWGT  ■ 0*0 

SUIT  a 0*0  " 

DO  400  I ' liNPARTS 

BUBYWGT  a BU9vwGT  ♦ weiGHTtl) 

SUmT  • SUmT  a weiGHTm.DISTCGKm 
CONTINUE 

BUOYLGK  • SUMT/HUOYWGT 
DEPTHCG  * DEPTHK  • bUOYCGK 

COMPUTE  THE  PArT  BASIC  MOMENT  OF  INErTIA  CONTRIBUTION 
SUMT  ■ 0*0 

CALL  aOOYMl  ( ISHApE ( I) ,HEIqHT ( I) , uI°TM( I ) ,TH£CK( I ) aPINERT • 
F"R  THE  BODY  ABOUT  ITS  0WN  AXIS 
PM  1 a WEIGHT ( t )*PJNERT/G 
A*8UT  the  CG 

uMICOMr  a ( weIqHT ( I l/y) »IDISTCgkI I>  • bU0YCqk)««2 

gyMT  a SUMT  ♦ PMJ  ♦ BMICOMP 

I continue 

BUOynl  ■ SUmT 

get  displacement  contributions 
CALL  DISPLACE 

MUOYLGCO  a HU9YCBK  • BUOYCGK 

WRITE  ( NiP,  R4  05  I ByOYwGT,WDlSPl.A{; 

COMPUTE  DISTANCE  TO  M£TACENTER  FR©m  CB 

SURF INRT  a PJ aRWL**4/4f 0 
BU9YLBM  a SURFjNRT*RH8G/WDISPLaC 

COMPUTE  RIghTINq  ARM,  g* 

BUOYLGM  • 9U0YCGCB  ♦ BUOYCBM 
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t 


i 


* 


130. 

121. 

122. 

123. 

124. 

125. 

126. 
127. 
120. 

129. 

130. 

131. 

132. 

133. 
1 34. 

135. 

136. 

137. 

138. 

139. 
140- 
141* 

142. 

143. 

144. 

145. 

146. 

147. 

148. 

149. 
150« 
1*51  • 

152. 

153. 

154. 

155. 

156. 

157. 

158. 

159. 

160. 
161. 
162* 
163. 
164  . 

165. 

166. 

167. 

168. 

169. 

170. 

171. 

172. 

173. 

174. 

175. 

176. 

177. 

178. 

179. 


C 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


X 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

X 

c 

c 


c 

c 


oec*  for  stability 

IF  (dUfiVCfjM  «LT  • 0*0)  W»ITE  !NLP#9700>  BUOYCbk+bUOYCbM# 

♦ BU8YCQK  | 

♦ STOP  550 


CPlPyTE  WIGHT INQ  MOMENT  TERM 

HUOYMR  ■ BUM YWGT  »OUO YCQM 
0U8TMRG  ■ .BJ9YHR/G 

C4MPUTF  NATURAL  PERIOD  OF  ROLL 
FIRST  GET  VIRTUAL  MOMENT 
V J R f J NRT  ■ HU9YMJ  ♦ AODMJ 
PER1OQ0  ■ 2 • 0 »P I *SWRT  ( V I RT I NRT/BuO YM* ) 
wRIIt  ( NLP i 94 1 0 ) PERI  BOO 

ma«E  ASSUMPTION  8F  UNIT  AMPLITUDE 
AND  START  FREQUENCY  ANALYSIS 

IF  ULIST  .EG.  i>  WRITE  (NLP#945o> 
rrsintg  ■ 0.0 
SINTg  » 0.0 

OUTPUT  MUflYMI .BUOYMW.auOYMRB 


NOTE  I THE  CALL  TB  MUBYOAMD  IS  PLACED  HERE  TO 

SIMPLIFY  THE  COMPUTATION.  THE  CONSTANT  TERM 
IS  COMPUTFD  HtRE  and  THE  FR£’)  IS  MULTIPLIED 
IN  AT  THE  BEGINNING  OF  THE  FREQUENCY  ITERATION, 


CALL  BUBYHamP 

DO  *000  TIME  • T I wEl#  T ImE2#  T ImFdEl 

FREU  . PJ  2*1 .OE*lO 

IF  ( T I Mr  , N£ , 0*0*  FREQ  ■ Plg/TIME 

WA  VgN  » FREQ*FRtU/Q 

•••  COMPUTE  Damping  MOMENTS 

SEE  NOTE  ABOVE 
MUOYMO  a HUOYMQTaFHEQ 
CALL  WATEROAMP 

***  COMPUTE  WAVE  INERTIA  MOMgviTS 
CALL  hATRINRT 

OUTPUT  BUOYMO.wATERMDaWATERIH  4 


ETfcM  a BUOYMrg  ♦ WATERIP 
FT£M  a taATERMO 

SET  MAIN  COMPONENT 
E • ETEM*FR£Q*FR6U 
F » FTEM»FWeQ 

PHASE  BETWEEN  WAVE  AND  TORQUE 
SIGMA  ■ RTOO* ATANg  ( *E«F ) 

exciting  TORQUE 
T • SORT  (E*E  ♦ F«F> 


C 

C 


- 94 


180* 

IM. 

182. 

183. 

184. 

185. 

186. 

187. 

188. 
189. 
190* 

191. 

192. 

193. 

194. 

195. 

196. 

197. 

198. 

199. 

200. 
201. 
202. 

203. 

204. 

205. 

206. 

207. 

208. 

209. 

210. 
211. 
212. 
213. 
21*. 

215. 

216. 

217. 

218. 
219. 
220* 
221. 
222. 

223. 

224. 
225* 
226. 
227. 
2*8. 
229. 
230- 
231. 
212. 

233. 

234. 

235. 
236* 

237. 

238. 

239. 


C 

C 


c 

c 

c 

c 

c 

c 


PWASE  ANqLE  between  torque  AND  ROLL 
ATtM  ■ RUMVM0»f WEQ 

OTfcM  ■ BUOYMR  . V|9T!NRT*FReQ*PReG 
PHI  ■ RT90.ATAN2  (-ATEM.bTEM) 

ROLL  RAO 

RtfLLRAO  » T/SfJRT  (ATEM*AT£M  4 BTEm*STEH) 

PHASE  ANeLE  BETWEEN  WAVE  AND  ROLL 


c 

c 

c 

c 


1400 


c 

c 


c 

c 

1500 

c 


theta  » SIOMA  * phi 

GET  SEA  SPEcTHA 

CALL  SEASPEC 


COMPUTE  RESPONSE  AND  INTEGRATE 
PWS  ■ W4LLRA0*RBLLRA0*S 
I»-  (TIME  .LE.  TIMEll  00  TO  1400 
DELp  ■ pRPULAST  • pH£Q 

RRSJNTG  ■ (HRS  ♦ RRSLAST ) *0»50*0eLP  ♦ RRSINfG 
SINTQ  -(S-.  SLAST>#0.50*OEtr  * sintg 
continue 
RHSLAST  ■ RRS 

SLAST  • S __ 

FNfcULAST  • pHEQ 


OUTPUT  LIST  IP  IT  WAS  SELECTED 
!►  (ILIST  .LE.  0)  SO  TO  1500 

white  (NLP. 94551  TIME.^EU.^OLLRaO^RtOD^SIGMa^pHI.tNETA.S 


CHECK  FOR  PLOT 

CONTINUE 

IF  t I PLOT  .LE.  0 1 SO  TO  2000 

call  plots  < ••',rbllrao*Rtoo,time.1> 


2000  CONTINUE 
C 

c GET  STATISTICS 

c 

c COMPUTE  ROOT  MEAN  SQUARE  OF  WAVE 

RMS  ■ SQRt(SINTG! 
pR08AMp  • o«707*RmS 

c COMPUTE  AVERAGE  WAVE  HEIGHT 

08  2300  I ■ 1.5 

AVXESPNsm  • AVRC0EFF(  I>»RmS 
2300  CONTINUE 

C COMPUTE  MAXIMUM  WAVE  AMPLITUDES 

08  2400  I • 1.6 

ROu.mAX<!>  « *vMAXCBF( !>*RmS 
2400  CONTINUE 

wR I Tfc  (NL».9500»  RMS 
“RfTfc  (NLP, 95051  PR08AMP 

WRITE  ( NLP. q5 10  1 (FHACAmpS( i ), AvRESpNS( I ), 

♦ maxwAVN8(1>,RPLLMAX{ I),  I ■ i , 5 ) « 

♦.  MaxwavN»<6>#R0LLMAX(6> 

C 

C COMPUTE  ROOT  MEAN  SQUARE  8F  RESPONSE 

RM8  ■ SORT (RRSINTG) 


t 


BR80AMP  ■ 0 • 707*WMS 

C COMPUTE  AVERAGE  RESPONSE  OF  rBll 

ue  aa oo  I > i « 5 

AVWtSPNSff)  . AVHC9EFF ( I ) »NHS«RTB0 
2500  CHNUNUE 

c COMPUTE  maximum  AMPLITUDES  8F  RBlL 

UB  2600  I a 1*6 

ROLLM A X ( 1 ) . WVMAXCBF(  I > •PMS*RT80 
2600  CONTINUE 

-PITS.  <Nlp,P515)  HmS»RT8D 

*RlTfc  (MLP|P520)  PWBHAMP*RTOD 

want  (NLP/9525)  IFRACAMPbm,  AVPir^PNSlII* 

♦ MaXwavnBI p , PBLLMAX ( I ) t 

♦ mAXi*AvNB(6)/RBLLmAX(6) 

CHECK  F8P  PLOT 

IF  UPL8T  ,lE.  0)  GB  TO  3000 

call  plbtroll 

3000  CONTINUE 

IF  (I8NFLAG  *EQ«  1)  WRITE  (NLP/94B5) 

W£A 0 <NCP,90B5,FMO»8000>  lE^O 
IF  IltNO  .Erj.  1HY)  qO  T8  150 

’*000  CONTINUE 
STBP 

9n85  FORMAT  (au 


9400  F8P"AT  <1^1/ ' 


g405  FORHat  </'  CM^CK  F»R  BUOYANCY  HUBY  W£jQHt  » ',F10,1,'  LBSM ' t 

♦ ' WATER  DISPLACED  ■ ' / F 10 • 1/  ' LBSm') 

9410  FORMAT  (/.  NATURAL  PERIOD  *'/Fg.3,»  SECBNDSi) 

9450  FpiRMat  UHl/'  PaB  IS  IN  DeGpEES/FB8t  8F  *AVe  AMPLITUDE'/ 

♦ ' PERIOD  ANci  FREQ  RAB  *.T  PHASE  T-«  PHASE 

♦ <w.R  phase  amp  SPEC  ») 

9455  F9RMAT  < F j o . 3,  E 1 0 . 3,  F i 0 . 3, F 1 0. 3 , F 1 o. 3/ F 1 0 .3, Fl 0 . 3 > 

9485  fwwmat  (/I  DS  V*U  WANT  ANOTHER  CASE  B*J 

9500  FBRMAT  < IHi //// * RMS  BF  WAVE  SPECTRUM  •»/FlO*3»'  FEET  I ) 

9505  format  {/,  PRBBABLE  AMPLITUDE  '/ 

♦ ' BF  WAVE  • ' iF10*3a • FEET* > 

9510  F8RMAT  <// 

♦ /»  FRACT18N  BF  I EXPECTED  ' 

♦ /'  LARQFST  AVERAqE  I WAVE  » 

♦ /*  AMPLITUDES  wave  I NUMBER  MAXIMUM  i 

♦ /'  CONSIDERED  AMPLITUDE  I BF  WAVES  AMPLITUDE  * 

♦ / ' -----------  I ........  .........  i 

♦ b(/TA,F5-3.  T16#F9*3#  Tg8,'I'#  T33.I6/  T43»F9.3) 

♦ / T 2H/ 1 1 * i T33/I6,  T43.F9.3  ) 

951?  f8rmat  <///////*  kms  8F  response  spectpum  «',fio.3,'  deq  ») 
9520  FORMAT  I/*  PR8pAdLE  AMPLITUDE  */ 

♦ • BF  MOLL  RESPONSE  ■ ' »F1o*3#  » DEG  ') 

9525  F8rmat  <// 

♦ /•  AMPLITUDES  BF  RBLL  are  IN  DEGREES*/ 


9410  format 
9450  FPIRMAT 


9455  F9RMAT 
9485  fwwmat 
9500  FBRMAT 
9505  format 

A 

9510  FORMAT 
♦ / 


RBLL  RESPONSE  ANNALY3!,'S  PROuRAM*/ 

ALL  DFPTMS  ARE  POSITIVE  '/) 

CM^CK  F»R  BUMYANCY  HUBY  Weight  » ',F10,1,'  LBSM ’ / 

w*TER  DISPLACED  ■ '/FlO. lx'  LBSm') 


NUMBER 
BF  WAVES 


FEET* ) 


EXPECTED 

WAVE 

maximum 

AMPLITUDE 


T33.I6/  T43.F9.3) 
T33/I6/  T43.F9.3 


/'  FRACTION  BF 

/’  largest 

/'  AMPLITUDES 
/*  CBNSIOERED 


AVERAGE 

ROLL 

amplitude 

RESPONSE 


NUM9ER 
BF  WAVES 


EXPECTED 

R8LU 

maximum 

amplitude 


- 96  " 

# 

✓ 

300* 

♦ /'  1 ' 

301  • 

♦ 5I/T4.F5.3.  T16/F9.3*  T28.*I*.  T33.16#  T43.F9.31 

aoa* 

♦ / T25.  ' I * . T33.16.  T43*F3.3  //) 

303. 

C 

30**  • 

9700 

FORMAT  (//•  STB**  EVERYTHING  . THIS  BUOY  WILL  ROLL  OVE**/ 

305. 

♦ « THE  CENTER  BF  GRAVITY  IS  ABOVE  THE  "ET**/ 

306. 

♦ » - CENTER.  MC  ABOVE  KEEL  * * . F6 * 2. * FEET*/ 

307. 

♦ t CG  ABB VE  KEEL  « *.F6*2. ' FEET*) 

308. 

C 

309. 

c 

310. 

c 

311. 

c 

312. 

SUBKOUT  I Nfr  SEAS^EC 

313. 

I c USEASEL  .EQ.  oJ  S * 1*0  1 RETURN 

31**. 

FR£UM4  * CPirO**** 

315. 

*•  REUHb  « fREUR4*fRE0 

316. 

OB  TO  ( 4 100*  42oo*  **300  > . ISEASEL 

317. 

c 

PIERSON  . MBSKBWIT4 

3 1 8 • 

4100 

CONTINUE 

319. 

S • l35.0/FREOp5»E*p(-97000.0/tFHEQP***wIN[)VP**) ) 

320. 

c 

CORRECT  FBR  DBUBLE  HEIGHT  SPECTRUM 

321. 

S « S/8.0 

322. 

RETURN 

323. 

c 

gRETScHNEIOER 

32**. 

**200 

CONTINUE 

325. 

S • 4200.0#wavEhTP2/(WAVPERP4*FREOP5>* 

326. 

* EXP(.1050.0/(WAVH(rRP4*FP(rUP4)) 

1 

327. 

c 

CORRECT  FOR  DOUBLE  Hr  I qht  SPECTRUM 

328. 

S • S/8* 0 

329. 

RETURN 

330. 

c 

i • s « s • c • 

331* 

4300 

continue 

332. 

S ■ 2760.0#w*VFWTP2/{WAVPEHR4*FREQP5)<* 

333. 

♦ pXP(  »630.0/1  W*VPeRP4*FRe^F,4)  1 

33**. 

C 

CORRECT  F0°  DOUBLE  HEIGHT  SPECTRUM 

335. 

s ■ s/a.o 

336. 

return 

337. 

C 

338. 

c 

339. 

c 

340 . 

c 

3**1  * 

subroutine  plbtimt 

3**p. 

c 

this  subroutine  is  used  to  initialize  a line  printer 

3**3« 

c 

PLOT  of  THF.  ROLL  RESPONSE. 

3****. 

c 

SUE  LIMITED  fop  on-line  USE  ONLY 

3<*5» 

c 

3**6  * 

DIMENSION  IPL0TMUF1 1300» 

* 

3**7. 

c 

3**8  • 

DATA  IbUFSUE  / 1300/ 

349. 

c 

350* 

c 

351. 

JF  (HMAX  ,GT.  Rm I N ) GO  TO  1490 

352. 

rma*  ■ ioo*rmjn 

9 

353. 

If  IRmAX  .E0»  0»0)  RmAX  • 10*0 

354. 

1 **9Q 

continue 

355. 

ROEL  • Rm AX  - RMJN 

356. 

hmax  • IF JX<R0EL/5»0  ♦ 0.9991*5.0  ♦ RMJN 

357. 

c 

358. 

NLINfeS  • (TIME2  - TJMEl l/TIMEOEL 

359. 

NBaRS  • (NLINeS  + 9.11/10.0 

k 

| 

- 97  - 


TnAX  • TJmM  ♦ N0ARSMO*O*TImEdEL 

NLlNtS  • (TMAX  - T1ME.1  )/TIHEDEL  ♦ 1*0 

jjr  (NtlNE?  .LT,  IBUfSUE/131  GP  TO  5000 

Ip^ht  ■ o 

wRtTfc  (NLP/9710J 

RETURN 

CONTINUE 

CALL  HL^Tt  (N0ARS* 10*3* 10) 

CALL  plots  <I0L0TBUf>«*I*/WMAX,T"AX, TIMED 
RETURN 


FORMAT  (/i 
► / ' 


THE  HLtJT  BUFFER  IS  NOT  LARGE  ENOUGH  * 

fur  the  period  range  specified  ' 

THE  PLOT  IS  SUPPRESSED') 


subroutine  plotroll 

THIS  SUPR8UT I N'E  IS  USED  T9  9UTPUT  TM£  LINE  PRINTER 
PLOT  of  the  roll,  response. 

CALL  PLPTe;  ( 3, 3-*,  i ROLL  HFSPBNSE  AMPLITUDE  OPERATOR  i ) 

Call  PL0TA  ()*,'  PeRJBD  (SEC)  m 

CALL  PLOTS  <S,2*»»  DEGREES  PER  FOOT  AMPLITUDE  •) 

CALL  PL0T7  <1Q) 

RETURN 
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1. 

2. 

3* 

4 « 

5. 
6 « 

7. 

8. 

9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 
21. 
?2« 
23. 
2*. 

25. 

26. 
37. 
28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 

41. 

42. 

43. 

44. 

*5. 

46. 

47. 

48. 

49. 

50. 

51. 

52. 

53. 

54. 
55* 

56. 

57. 

58. 

59. 


C 

C 

c 

c 

c 

c 


c 

c 

c 


c 

c 

200 


c 

c 

300 


c 

c 

400 


c 

c 

500 


SUBROUTINE  TJNPUT 

VERSION  l.Q  S£P#  1976  R.  Q8LDSMITH 

this  routine  iwputs  data  fbr  the  pull  pao  in  on-line  hade 


I BOEV  / NCR, NLP 

TP  / TIHfcl,TtHE2<TIMEDEUJFPEU#WAVEN 
SE4STATE  / !5EASEL'*lNDViWAVEHT/wAVEPEP# 
w I NQVP4* * AVFHTP3# WAVPCPP4 

P ! NS  / NHARTS,  IPHAPE(656>.WIDTH<50)  # H£ I GHt ( 5o ) , TH J C* ( 50 > # 
DtNSITV(50>#0ISTcGK(50j#EPACN9PM(50) 

POUTS  / VOLUME ( Po  > » bE I GHT { 50 ) 

CONSTANT  / PI*RTaO.KH8,G 

BU^Y  7 NpnAX.WrfL. AVtRQAMP#THET*HAR.PEP!BO0# 

DEPTHK#bU9TCb«#DEPTHCQ#8UBYCBK,0EPTHC8,W0ISPUAC 
5UTPUTS  / IL1ST* IPUOT^RHIM^RMaX 


INPUT  time  AND  range 
aWITfc  (NLf-«9400) 

R£AU  ( NCR  1 9000  > Tl#T2#T3 
IP  < T l .UT.  0*0)  G®  T»  200 
T I HE  1 ■ T l 
TIMfc.2  • T2 

If  (TIME2  .UT.  TIMtl)  T IME2  • TIME! 

IP  < T 3 .LE.  0*0)  13  • TJMEDEL 

TI*tUEl  « T3 

SELECT  SEA  STATE  PARAMETERS 

CBNT1NUE 

k-  p?  I T fc.  ( NLP/  9405  ) 

PEAL'  (NCR, 9005)  1STEST,  WAVtHr,  WaVePE* 

IF  (1STEST  .qT.  3)  JSTEST  • *1 
IP  (1STEST  .LT«  0>  G®  TO  300 
ISPASEL  » ISTEST 

IP  (1SEASEL  «E0.  1>  WJNDV  « WAVEHT 

W I NO VM*  • W I NO V**4 

t»AVtHTP2  • wavemt*i«aveHT 
wAVMt«H4  • wAvEpEW»*4 

ENTEP  WATER  PLANE  RADIUS 

CONTINUE 
"RITt  ( NLP#  P4 1 0 ) 

«E4U  (NC«, 90105  RRLTEST 
IP  (HwLTEST  *LT.  0*0)  GO  TB  400 
RWL  • RwLTEST 

INPUT  DEPTH  TB  KEEL 

CONTINUE 

WRITE  (NLP#R415J 
REAU  (NCR, 9015)  4*TEST 
IP  UkTEST  .LT.  O'O)  Q®  TO  500 
DE»*Th<  ■ ZKTEST 

ENTER  ESTIMATED  AVERAGE  AMPLITUDE 

CONTINUE 
WRITE  (NL«»#9420I 


COMMON  / 
COMMO!*  / 
COMMON  / 

COMMON  / 

Common  / 

COMMON  / 
COMMON  / 

COMMON  / 


» 


/ 
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M 


* 


60. 

61. 

68. 

63, 

64. 

65  » 

66  . 
67. 
68  . 

69. 

70. 

71. 

72. 

73. 

74. 

75. 

76. 

77. 

78. 
79 1 
8o« 
81. 
82* 

83. 

84. 

85. 

86. 
87. 
B8. 
89. 
90* 

9 1 . 

92. 

93. 

94. 

95. 
96* 

97. 

98. 
39. 

100. 

101. 

102* 

103. 

104. 

105. 

106. 
107. 
10«* 

109. 

110. 
Ill* 
112. 
113* 

114. 

115. 

116. 
117. 
11®. 
119. 


R£  AD  (NCR, 9020)  AMPTEM 
IF  <Ai“IPTEM  "LT.  0*0)  QO  T9  600 
AVEKUAMP  ■ amptem 

C FMTER  ESTIMATED  AVERAGE  ROLL 

600  CONTINUE 

wR  I Tt  ( NLP , 94^5  > 

READ  (NCR, 9025)  T*ETATEM 
IF  (THETATEM  .LT.  0*0)  G®  T®  700 
THETABAR  > THETATtM/HTBO 

C 

c 

700  CONTINUE 
C 

c DEPINE  BUOY 

C 

N ■ 0 

1000  CONTINUE 

*Rirt  (N(_p,9455) 

HEAD  ( NCR, 9055  J NPTEM 
IF  (NPTEM  «LE • 0>  Q®  TB  2000 


wwjTt  ( NLpj 9460 ) 

N • 1 ' - 

c 

c L89P  ON  ENTRY 

1050  CONTINUE 

IF  (N  ,GT,  NPMAX)  WRITE  (NCP.9725)  NPmax  I 
WHJTE  (NLP, 9465)  N 

READ  (NCR, 9065)  IDUM  - - 

wRITt  ( NLP l 9A  7q ) 

read  (ncr, 9070)  <,w,h,t,D#x»p 

C _ . . _ U r i Alktlur,*  fill  n ( 

c 


5) 

.PR.  ( K 

•LT*  1)) 

OUTPUT  K 1 
GO  TO  1050 

0.0) 

WRITE 

l NLP* 9705’ 

) 

OUTPUT 

w 

0«0> 

WRITE 

( NLP*  9705 ) 

1 

OUTPUT 

M 

0.0) 

• 0 

WRjte 

(NLP, 9710) 

1 

OUTPUT 

o*  0 ) 

WRITE 

( NLp*  97j 5 ) 

) 

OUTPUT 

DjSfLGK(N)  » * 

(iO  TO  (1100*1100*1100/1400*1500)*  K 
cylinders 

1100  CONTINUE 

IF  (T  ,£U.  -l»o>  T ■ P*6«0 

IF  ( ( K ,eQ,  2)  .SR*  * * ‘E0*  3*  > T ■ w*5»0 

IF  (T  .Ct.  0.0)  -RITE  (NLP/9705)  1 OUTPUT 

TwlCMN)  • T/J2.0 
FR^CnbRMJM  ■ 1,0 
ISMApE(N)  • 1 

GO  Tt>  1800 

triangle 

1A00  CONTjNUEy^  WRITE  (NLP*9705)  I OUTPUT 

THJCMN)  ■ T/12,0 


GO  TO  200 


WRITE  (NLR/9700)  I 
w i GO  TO  1050 
M , GO  TO  1050 
D J Go  TO  1050 
X / GO  TO  1q5q 


T | GO  TO  1050 


T | GO  TO  1050 


^ . 

- 100  - 

• 

T 

120. 

IF  <(F  .Gf.  1.0)  * 8R»  (F  *LT.  0*0))  WRITE 

(NLP'9720)  j 

OUTPUT  F 

121. 

♦ I 08 

TO  1050 

122. 

FRACNSRM(N)  ■ F 

123. 

ISWAW£(N)  • 2 

124. 

G8  TtJ  1800 

125. 

C 

126. 

C 

RECTANGLE 

127. 

1500 

continue 

128. 

IF  (T  .LT.  0*0)  WRITE  <NLP*97o5>  ) OUTPUT  T 1 GO  TO 

1030 

129. 

TRICMN)  • T/12.0 

130* 

IF  (CF  .qT.  1.0)  .OR.  (F  .LT.  0*0))  WRITE 

( NLP*  9720 ) j 

OUTPUT  F 

131. 

♦ j GO 

TO  1050 

132. 

fracnorm(n)  • f 

133. 

ISHAHL(N)  a ? 

134. 

GO  TO  1800 

135* 

c 

136. 

c 

COMPUTE  VOLUME  AND  WEIGHT 

137. 

1800 

CONTINUE 

118. 

C*LL  SOoTVOi.  ( ISHApE(N)*h,w,THIc<(N)*V) 

139. 

weIGHT(N)  a V«D 

uo. 

V9LU«L(N)  • V 

141. 

c 

142. 

c 

check  number  of  entries 

143* 

IF  <N  ,GE.  NPTEM)  GO  TO  2000 

144. 

N • N ♦ 1 

145. 

GO  TO  1050 

_ 

- 

146. 

c 

147. 

c 

PART  CHANGE 

- .. 

148. 

c 

149. 

2000 

CONTINUE 

...  _ 

150. 

IF  (N'TEM  ,GT.  0 > NPARTS  ■ NPTEM 

151. 

NPTfc.M  ■ o 

... 

152. 

Np ah r s ■ max  <noarts*n> 

153. 

wRJTt  < NLP#  9475 ) 

154 . 

HEAD  (NCH.9075)  N 

155. 

IF  (IN  .Gr.  0)  .ANO.  IN  ,LE*  NPa«TS  ♦ l)> 

38  f 8 1050 

156. 

c 

157. 

3000 

CONTINUE 

158. 

WR J Tfc.  ( NLP#  9480 ) 

159. 

HEAD  (NCR, 9080)  IL# IP#«mIN,RMaX 

160. 

ILI5T  > 0 

161. 

IF  (IL  .EO.  IHYj  ILIST  » 1 

162. 

I PLOT  • 0 

163. 

IF  UP  *E0.  lHy)  IPLOT  ■ 1 

164. 

C 

< 165. 

RETURN 

166. 

C 

157. 

9000 

FORMAT  (3F.0) 

168. 

9005 

FORMAT  (J, F.O, F.O) 

169. 

9010 

FORMAT  (F.O) 

170* 

9015 

FORMAT  (F.q) 

171. 

9020 

format  (F.O) 

172. 

9025 

FORMAT  (F.O) 

173. 

9055 

FORMAT  (I) 

174. 

9065 

FORMAT  (A) 

• 

175. 

9070 

FORMAT  (I,6F»0> 

176. 

9075 

FORMAT  (I) 

177. 

90«0 

FORMAT  ( A1,1X,A1#1*#2F.0> 

178. 

C 

179. 

VO 

o 

o 

rORMAT  (/»  ENTER  START#  END#  INCREMENT  of 

PERIOD  RANGE 

(SEC)  8*) 
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ENTER  SEA  SPECTRUM  TYPE  AND  PARAMETERS' 

■ 1 * 0 O' 

p1F.RS0n»m*SkBwitz  1*  WIND  SPEED  (*N8TS)* 
HRETSChNEJDER  2$  SIGMIE  WAVE  WT  (FT)# 

wave  PERIOD  < SEC ) • 

t.S.S.C.  3,  SIQNJE  WAVE  HT  (ET)« 

WAVE  PERIOD  ( SEC  > ' ) 

ENTER  "ATE*  plane  radius  at  SUREACE  (ET)  8 • ) 
£MTr<5  DEpr*  TO  *EEL  <ET)  8»> 

ENT^p  EXPECTED  average  AMPLITUDE  (ET)  0 » ) 
ENTER  EXPECTED  average  W6UL  OE3>  8*> 

£\TEP  \Um«E*  OE-  BUftV  Pa«TS  8') 

• ••  raw  each  p aw t number  ybu  MUST  ENTER  '/ 
S8me  IDENTIFIER  BEFORE  YBU  return. »/ 

THEN  ENTER  k,w,h,T,D#X#E 


sign  IE » 
SIGN  IE* 


ft  • SHAPE  C9DE  1 • HOLLOW  CYLINDER  '/ 

? - SOLID  CYLINDER  •/ 

3 • DISC  '/ 

A • TRIANGULAR  (RT)  PLATE  '/ 

ft  - RECTANGULAR  PLATE'/ 

m • wjDTM  OR  OUTSIDE  DUMeTfR  ( P T » •/ 

h ■ HEIGHT  C pT  > '/ 

t . thickness  (In>  '/ 

A *x  ENTERED  EBR  Case  <»l  WILL  ASSUME  '/ 

A SBLlD  ,T  “ w/?/l?>  'y 

FOR  C»SES  k.2#3  ENTER  ANYTHING  '/.... 

D • DENSITY  ( LSM/FT **3 ) '/ 

* ■ DISTANCE  Er8m  KEEL  Tp  PApT  CG  C FEET > */ 

E . FOR  PLATES  ONLY#  FRACTIONAL  AREA*/ 

ME  TWE  PLATE  Normal  TO  MOTION)/ 

enter  i ebr  cylinders  '///) 

part  N8.  *1  13#  *0* ) 

ENTER  K# w#H# T#D# X#E* ) 

ENTER  PART  NIJMWEH  TO  CHANGE  (•*  TS  STOP)  0*) 

Enter  y *r  n e»r  list  and  plot  »ptibns'/ 

FOR  PLOT  YOU  MAY  ALSO  ENT£R  RAO  MIN  and  max  0*) 

•••  YfiUHE  KIDDING  • THE  CODES  ONLY  GO  FROM  1 TO  3 1/ 
try  again  */) 

•••  WHAT  KIND  be  SHApE  IS  THIS  O') 

•••  what  do  you  have  IN  THERE  0*> 

•••  HHtRE  IS  IT  O') 

• ••  R A NyE  »E  f ■ 0*0  TO  1 • 0 o*) 

• ••  m\l Y • # 1 3#  « COMPONENTS  ARE  ALLOWED  I / 

buoy  definition  terminates  '/ 
buOv  cH*NGES  WILL  proceed  '/> 


HOLLOW  CYLINDER 
SOLID  CYLINDER 
DISC 

TRIANGULAR  (RT) 


fl 
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4* 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

10 

15 

16 

17 

18 

19 

20 
21 
22 
23 
20 

25 

26 

27 

28 

29 

30 

31 

32 

33 
30 

35 

36 
17 
3« 
39 

00 

01 
02 

03 

04 

05 

06 

07 

08 
09 

50 

51 

52 

53 

50 

55 

56 

57 

51 
59 


SUBROUTINE  BJNPUT 

c 

C VERSION  1*0  SEP#  ] 976  R»  GOLDSMITH 

C TWJS  ROUTINE  INPUTS  DATA  FOR  TWE  ROLL  RA9  IN  BATCH  MOPE 

C 

C 

C 

C 

0 1 NfcNS 1 ON  I COMMENT ( o ) # 10(10) 

C 

Common  / {OOEv  / NCR# NLP 

COMMON  / TP  / T|ME|,  T!^E?jTIOEDFL,ERE0/WAVEN 
COMMON  / SEASTATE  / I SEASEL#  *•  INpV,  RAVFHT#  WAVf.PgR, 

♦ „,INUVPOiwAVENTPP#^AVPERPO 

COMMON  / BINS  / mharts#  ISHAUfc  (5O)#WIDTH,5O)»^I0HT(5O)#THIC<{5O)# 

♦ OENS1TV(50>#OISTCGK(50)#FRaCN9RM(50) 
common  / pOUTS  / vOLUMfc(50>#wEIr.HT(50) 

COMMON  / CONSTANT  / PI#RTB0iRW8#0 

COMMON  / BUM*  / NHMAX1H*t.iAvERGAMP,TMETAyARjpEliIBDO# 

♦ DtPTHK#My9rcC>,<#DEPTWcG#BuBYCB<#0EPTHcB#  wDISPl^C 
common  / outputs  / ilist#  ipi»t.rmjn#rmax 


DATA  ID  /’  H CTL  '#•  S CTU  »,»  DISC  •#*  TPJ  PUT'#'  RCT  PLT»/ 

C — 

C INPUT  TIME  AND  PANCIE 

RE AU  ( NCR* 9005 ) T 1# T2# T3 
IE  ! T 1 » LT » 0.0)  GO  TO  178 

TIMfil  • TJ  • - - 

TIME2  • T? 

IE  ( T I M£2  ,CT.  TIMED  TIM£2  • TIMcl 

IF  ( T 3 .LE*  0*0)  T 3 ■ TJmEDEU 
TIMEUEL  • T3 
175  CONTINUE 

WRITE  (Ntp.9005)  TIME1#TIME2#TIMEDEL 
C 

C SELECT  SE*  STATE  PARAMETERS 

200  CONTINUE 

REAP  (NCR. 9010)  JSTEST#  WAV£HT#  WAVEPE* 

IE  (JSTEST  »GT.  3)  JSTEST  • -1 
IE  (JSTEST  .LT,  0)  G»  TO  275 
ISEAPEL  ■ ISTEST 

IF  (JSEASEL  ,E0.  1)  WJNDV  ■ WAVEHT 
WJNUVP4  ■ W f NOV*** 

WAVEMTP2  ■ wAVEHT**AVEMT 
WAVHERP4  ■ WAV£PER«o* 

275  CONTINUE 

IF  USEASEc  •EG*  0)  WHITE  (NcP#9410> 

IF  (ISEASEL  «EQ.  1)  WRITE  <NLP#9<H1)  WINDV 

IF  (1SEASFL  »EQ.  2)  RHITE  (NLP#9*12)  WAVEHT.HAVEPER 

IF  (ISEASEL  * £Q»  3)  wHJTE  (NLP/9413)  w* VERT # W*VEPER 

C ENTER  WATER  PLANE  RAOIUS 

300  CONTINUE 

REAP  ( NCR. 9015 ) WwtTEST 
IF  (RwlTEST  .LT*  0*0)  GO  TO  375 
hwl  • RWLTEST 
375  CONTINUE 
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60* 
61. 
62. 
63. 
64  . 

65. 

66. 
67* 
66. 


wRITt  (NLP. 9415)  R*L 
C 

C INPUT  DEPTH  T9  KEEL 

400  C9NMNUE 

Rfc'AD  ( NCR.  9020  > MTEST 
IE  UKTEST  .LT.  0.0)  GB  T9  475 
0£bThk  ■ 2KTEST 

475  continue 

wRITt  (NLP.942QI  depthk 


69. 

70. 

71. 

72. 

73. 

74. 

75. 

76. 

77. 

78. 

79. 
80* 
*1. 
82. 
83  a 

84. 

85. 

86. 

87. 

88. 
99. 
90* 

91. 

92. 

93. 

94. 
95* 
96. 


C 

C ENTER  ESTIMATES  AVERAqF  AMPLITUDE 

500  CONTINUE 

MEAL)  (NCW.9025)  AMHTEM 
IF  ( AMpfjTM  ,LT,  0.0)  G9  TB  575 
avE*uAMp  ■ AMpTEH 
575  CONTINUE 

*WITt  (NLC.9425)  avehuamp 
C 

C ENTER  ESTIMATED  average  R9LL 

600  C0NTINUE 

READ  (NCR, 9030)  TWETATEM 
IE  (THETATE*  .LT,  0.0)  39  TO  675 

thetaoar  ■ thejatem/rtbo 

675  CBNT1NUE 

WRITt  ( NLP. 9430 ) THETA8AR*WTBD 

c - - 

c 

700  C9NTINUE 
C 

C DEFINE  BUBY 

C 

N ■ 0 

1000  CONTINUE 

REAU  (Nc«. 9055)  npTEm 
IF  (NPTEM  ,UE«  0)  39  T9  2000 

C 

wRITt  ( NLP>  9455  ) NPTEM 


97. 

98. 

99. 

100* 
10  1 • 
102* 

103, 

104. 

■ 


117. 

118. 
119. 


wRITt  (NLP. 9460) 

N • 1 

c 

C L9BR  9N  ENTRY 

1050  CONTINUE 

IE  (N  .or,  NPMAX)  WRITE  (NLP. 9725)  NPMAX  t 
Read  ( NCR  # 9070 ) 8,R.h,T,D,*.E, ICbMMEnt 


C 

C 


c 

r 


SET  INPUTS  F9«  CBRRECT  SHAPE 


106. 

IF  <(K  .GT.  5) 

.BR*  (« 

•LT.  l ) ) OUTPUT  8 l 

107, 

♦ 

STOp  1050 

108. 

IF  (w  .LE.  0*0) 

WRITE 

(NLp*9705)  i 

OUTPUT 

109. 

w J 0 T ► : ( N ) a W 

11C, 

IE  ("  ,LT , 0.0) 

WRJTE 

(NLP.9705)  1 

BUrPUy 

111. 

ME  I ljHT  ( N ) • M 

112. 

IE  (U  *LT , 0*0) 

WRITE 

( NLP. 9710  J , 

BUTPUT 

113. 

0EN5»ITY(N)  • 0 

m. 

IE  (A  ,LT.  0.0 

wwjte 

(NLP.97I5)  i 

BUtPUt 

115. 

OISTLgk(N)  a X 

116. 

QB  TO  (1100*1100/1100* 

1400*1500)*  K 

1100  CBNT1NUE 


STBP  1050 

WRITE  (NLP. 9700*  I 
W I STPp  1050 
H l SyBP  1050 
D | STBP  1050 
X i STB*  1050 


cylinders 


104 


120 

121 

122 

123 

12* 

125 

126 

127 

128 

129 

130 

131 

132 

133 
13* 

135 

136 

137 

138 

139 
1*0 
1*1 
1*2 
1*3 
1** 
1*5 
1*6 
1*7 
1*8 
1*9 

150 

151 

152 

153 
15* 

155 

156 

157 

158 

159 

160 
161 
162 
163 
16* 

165 

166 

167 

168 

169 

170 

171 

172 

173 
17* 

175 

176 

177 

178 

179 


IF  (T  .EG.  -1*0)  T ■ w»6*0 
IF  I ( K . EO«  g)  .OR*  (<  *EQ*  3»  1 T • W*6.o 
| F (T  .LT,  0.0)  WHITE  (NLP. 9705)  J OUTPUT  T i 
THJCMM  • T/12»0 

FHACN(}RM(N1  • 1*0 
I SHAKfc.  ( -g  ) ■ 1 

GO  TO  1800 


9T0p  1100 


C 

c TRIANGLE 

1*00  CONTINUE 

IF  (T  .LT , o»0)  “RITE  (NLP»97q5)  I 
THILMN)  » T/12.0 

IF  (IF  ,Gt,  1.0)  .OR.  (F  .LT.  0.0)) 

V 

fRACNBRM(N)  • f 
ISWawe(n)  ■ 2 
09  TO  1800 
C 

C RECTANGLE 

1500  continue 

IF  (T  .LT.  o*0)  WRITE  «NLP*97o5>  l 
THJCMN)  • t/12.0 

IF  («F  • GT . 1*05  »9R.  If  »LT * 0*0)) 

F«aCN8RM(N)  • F 
ISMA^E(N)  « 3 
GO  TO  1800 


OUTPUT  T i STOP  i*00 

WRITE  (NLP. 9720)  J OUTPUT  F 
I STOP  1*00 


OUTPUT  T l STOP  i»00 

WRITE  (NLP. 9720)  l OUTPUT  f 
I STOP  1500 - 


C 

c COMPUTE  VOLUME  AND  WEIGHT 

1800  CONTINUE  

c*ll  aooYVBL  ( ishape(N).h.w#thick(n)#V) 
wEIliHT{N)  ■ V*D  

V9LUM£JN)  ■ V 


" WRITE  (NLP. 9*65)  N.I0(2*k-1)»ID(2*k)*W.H.THIC<(N).D»x*F 

+ .icomment 

c check  number  of  ENTRIES 

If  (N  «3E.  NPTEM)  go  TO  2000  . 

N . N ♦ 1 
GO  TO  1050 
C 

c PART  CHANGE 

C 

2000  CONTINUE 

IF  (NPTEM  .GT.  0)  NPARTS  ■ NPJEM 
NPTEM  * 0 

NPARTS  • MAX  (NPARTS. N) 

REAU  (NC«. 9075)  N „ 

IF  ( ( N .GT.  0)  . AND . IN  «LE*  NPARTS  * D>  WRITE  (NLP. 9*70)  4 
+ GO  TO  1050 


C 

3000  CONTINUE 

WEAU  (NCR. 9080)  ILa JP.KmIN.RMAX 
ILIST  • 0 

IF  UL  .EQ.  IHy)  ILIST  ■ 1 
I PLOT  ■ 0 

IF  (IP  ,£Q,  IHY)  I PLOT  ■ 1 

C 


C 


RETURN 
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1*0. 

9005  FORMAT 

<3F. 

0) 

1*1. 

9010  FORMAT 

( I.F 

•0»F«0> 

182. 

9015  FORMAT 

(F.O 

» 

183. 

9020  FORMAT 

(F.O 

) 

1*4. 

*025  format 

(F.o 

) 

185. 

9030  FORMAT 

(F.O 

J 

1*6. 

9055  FORMAT 

(I) 

187. 

9070  FORMAT 

( t.6F*0#3A4#A2) 

1*8. 

9075  FORMAT 

(I) 

189. 

90*0  format 

(AI,1X#A1#1X#2F.0) 

190* 

c 

191. 

9405  FORMAT 

(/* 

PERIOD  RANGE#  IN  SECONDS  START 

192. 

4 

2AX,FB,3#F*.3#2X,F8.3> 

193. 

9410  FORMAT 

(✓* 

sea  spectrum  . i»o* ) 

194. 

94U  FORMAT 

(/* 

PIERSBN-MOSKBWI Ti  SF*  SPECTRUM  * 

Ig5. 

4 

/ • 

WIND  SMLED  ■ ',Ft0.3,'  KNOTS') 

196. 

9412  f8PmAT 

(/» 

bretschneider  st. a spectrum* 

197. 

4 

/ • 

SIGNIFICANT  wave  ht  4 *#Flo«3# ' 

1 98 . 

4 

/• 

significant  wave  Period  • *,fio 

199. 

9413  fORmAT 

</• 

I . S • S * C * St*  SpEcTRuM* 

200. 

4 

/ • 

SIGNIFICANT  WAVE  HT  • *>F10.3#« 

201  * 

4 

/* 

SIGNIFICANT  wave  PERIOD  ■ *iFi0 

202. 

9415  format 

(/• 

wateR  flane  Radius  at  su»faCe  ■ *,i 

233. 

9420  format 

(/* 

DEPTH  to  the  <EEL  ■ * #F6»2# ' FT') 

234. 

94??  FORMAT 

(/. 

ESTIMATED  avfraue  AMPLITUDE  » *#F6 

205. 

94 jo  format 

(/• 

ESTImAT£0  avErAqE  ROLL  ■ '*F6,2, ' 1 

236. 

9455  FBRMAT 

(/' 

NUMgER  8F  PARTS  • *#I3) 

207. 

9460  FORMAT 

(//* 

208. 

4 

* 

c.a.  •/ 

209. 

210. 
211. 
212. 

213. 

214. 

215. 

216. 
217. 
21*» 

219, 

220. 
221  • 
222. 
223* 

224. 

225. 

226. 
227. 
22*. 


end 


DELTA*/ 


FEET* 

*3,'  SEC*> 

FEET* 

>3#  * SEC*  ) 
■6.2,'  FT*> 

<2»  * FT*) 

DeQ') 


♦ 

4 

♦ 

* 

4 

4 

4 

4 


*# 


ABOVE  */ 

part  width  height 

KEEL  FHACT'/ 

N9  SHAPE  (FT)  (FT) 

(FT)  N(}KM  COMMENTS’/ 


THICK 

(FT) 


density*# 

(LBSH/  *, 

FT**3 ) ' # 


9465 
94  7(3 

*9700 

9705 

9710 

9715 

972c 

9725 


FORMAT 

FORMAT 

FORMAT 
FARMAT 
F HPMAT 
FORMAT 
F ORM  A T 
FORMAT 


END 


( 15# 
(/) 

(/' 

</* 

(/• 

(/' 

(/« 

(/' 

* 


2X#?A4#  F*»  2#  FB.g#  F I i • 4 # 


»•••*/) 
FlQ. 1# 


F7 » 2#  F7.2#  2x*4A4) 


• • • 
*•* 
• •• 
*** 
*•  • 


V8U«E  <IDDJNU  a THE  CODES  ONLY  Q0  FROM  l TB  S '/> 
PMAT  KINO  BF  SHAPE  IS  THIS  8 *) 
what  Ob  Tru  haV£  IN  TheSe  8 ') 

WHERE  IS  IT  0 • ) 

RANGE  OF  r « o*0  T®  1*0  n *> 

ONLY  * # 1 3#  * COMPONENTS  APE  ALLOWED  '/ 

THE  PKBGPAM  TFPMJNATES  */) 


* 
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1, 

SUBROUTINE  5B0YMI  ( J S*H, W, T* PM J > 

2. 

C 

3. 

C 

VERSION  i.Q  SEP*  1976  R. 

goldsmith 

9* 

c 

5* 

c 

this  routine  computes  the  basic  shape  moment  of 

INERT 

6 « 

c 

ASSUMING  UNIFORM  DENSITY 

7. 

c 

8, 

C8MM8N  / J80£V  / NCR*  NLP 

9. 

c 

* 

13. 

c 

11. 

IF  (US  .qT*  3)  .OR.  US  .LT.  1))  WRITE  (Nt.P.9700) 

IS  l 

12. 

♦ 

STOP 

10 

13. 

c 

19. 

GS  TO  (100*200*3001*  IS 

15. 

c 

16. 

c 

CYLINDER 

17. 

100  CONTINUE 

18. 

WT  • w - 2«0*T 

1 9. 

IF  <*T  .LT.  0,0)  WT  • 0,0 

20. 

HMI  • (W*w  ♦ WT*WT)/16.0  ♦ H*H/ig.o 

21. 

RETURN 

22. 

c 

23. 

c 

TRIANGLES 

29. 

200  CONTINUE 

25. 

c 

•••-MUST  PLATrS  ARE  SMALL  SO  THE  INERTIA  IS  IGNORED  FOR 

NOW 

26. 

PM!  • 0.0  

27. 

RETURN 

28. 

c 

29. 

c 

RECTANGULAR  PLATE 

30* 

300  CONTINUE  . 

31. 

PM  I • W«H/t2»0 

32. 

RETURN  ' 

33. 

c 

39. 

c 

35. 

9700  FORMAT  (//•  what  KINO  OF  SHAPE  IS  CODE  *#  16/ 

36. 

♦ * PROGRAM  STOPS  I*  RHUtjNe  BODYMI 

•1 

37. 

c 

38. 

ENO 
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1. 

subroutine  sooyvol  <is#h,w#t#V) 

S' 

C 

3 ' 

C 

VE4SJ9M  i'O  SEP* 

1976  R.  G8LDSMITN 

4. 

c 

S' 

c 

COMPUTE  THE  BODY  VOLUME 

6 ' 

c 

7. 

COMMON  / f ODE V / NCRjNLP 

, 

A. 

COMMON  / CONSTANT  / PjiHT8Di«HB#a 

9. 

c 

10. 

c 

• 

11' 

IP  ms  'GT.  3)  .0*.  (IS  'L.T.  1)) 

WRITE  <NLP, 97001  IS  J 

12. 

STOP  10 

13. 

00  10  ( 1 00/ 200. 300  1 i is 

1*. 

c 

15. 

c 

16. 

c 

CYLINDERS 

17. 

100 

CONTINUE 

■ 

18. 

wT  ■ w « ?.n*T 

19. 

IP  <*T  .LT.  0.01  HT  • 0.0 

20. 

V ■ «»Pl*(w«*W  . WT*WT)/4»Q 

SI. 

RETURN 

22. 

c 

23. 

c 

triangles 

24. 

200 

CONTINUE 

. 

25. 

V ■ 0'5*H*W*T 

26. 

RETURN  . ...  ...  .. 



27. 

c 

26  * 

c 

rectangular  plate 

29. 

300 

CONTINUE 

30. 

V ■ m*w*T 

. . . . 

31. 

return 

32' 

c 

. ...  ...  . 

„ ...  ..  lt  . k 

33« 

c 

34. 

9700 

PORMAT  (//t  HHAT  IN  THE  WPRLD  IS  SHAPE  CODE  •»  16/ 

35. 

* ' program  st*ps 

JN  ROUTINE  bodyvol* 1 

36' 

c 

37, 

ENO 
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SUBROUTINE  (DISPLACE 


2. 

C 

VERSION  1»0  SEP*  1976  R.  38L0SMITH 

3. 

c 

« • 

c 

this  routine  is  used  to  compute  those  parameters 

5 • 

c 

6* 

c 

ASSOCIATED  WITH  THE  BUOY  DISPLACEMENT 

7. 

c 

8. 

C8*M0N 

/ BINS  / NPABTS*ISHAHE(5h)*WIDTH(5o)*hE!3hT(501*TwICK(50J 

9, 

♦ 

DeNSITV<50),DIStCG<(50)^RaCN8RM(50) 

10. 

C»"MON 

/ r9UTS  / VOLUME(50)*HEIr.HT(50) 

11. 

COMMON 

/ WATERDIS  / wD(bM*ND(56>*X0(5o>*V0(5O)*FD(50>*_ 

12* 

♦ 

DtPTH8(bO)#DEPTHT(50> 

13. 

COMMON 

/ CONSTANT  / PI *HT80»NH8«G 

14. 

CB^MON 

/ <3U0Y  / NHMAX.M^L.AVEPOahp^tnetaBaR.RERIBDo* 

15. 

♦ 

DePTMK^HUBVCUK^F.PTHCO.BUBYCBKjDERTHCB^WDISPLAC 

16. 

COMMON 

/ mShENTS  / BU9ymI*AOOmI#VIRTINRT*wATEPIm#BU8YmR# 

17. 

♦ 

BU8YMOT*BUBYHO,WATERMO,OAMPN 

18. 

c 

19. 

c 

20. 

PhBU  ■ 

RHB»G 

21* 

«dishlac  ■ 0*0 

22* 

SU*CBK 

• 0»0 

23. 

Aooni  ■ 

0.0 

24. 

c 

25. 

c 

gEGlN  loop  on  each  part 

26. 

c 

27. 

08  1000 

I • IMPARTS 

28. 

IS  - 

ISHAPEUJ 

29. 

M ■ HE  I QMT ( I ) 

30* 

W • WIDTH<IJ 

31. 

T ■ THIC«<  I ) 

32. 

V • V9LUMEU) 

33. 

* ■ DISTcaMJ) 

14. 

F • WEIGHT(I) 

35* 

<40  ( I ) 

■ w 

36, 

HO«IJ 

• H 

17. 

xom 

• X 

38. 

c 

COMPUTE  DEPTH  BF  PART  BOTTOM  AND  TBP 

39. 

c 

40* 

c 

41  . 

DCti  • 

Deptmk  . X 

42. 

c 

check  shape 

43. 

GO  TO 

(100*200*1001*  IS 

44. 

c 

CYLINDER  AND  RECTANGULAR  PLATE 

45. 

100 

continue 

46. 

Ttn  • 

H/2.0 

47. 

06  ■ 

DCG  ♦ TEm 

48. 

or  • 

ocg  • ter 

49. 

38  T8 

300 

50* 

c 

triangle 

si. 

200 

CONTINUE 

52. 

08  * 

OCG  ♦ H/3.0 

53. 

OT  ■ 

DCG  « H*0«6667 

54. 

c 

• •• 

THE  ABOVE  IS  BNLy  A GUESS  • CORRECT 

55. 

38  T8 

300 

56. 

c 

57. 

c 

IF  PART  TOTALLY  BUT  BF  RATER  I3NBRE 

58. 

300 

continue 

59. 

IE  (08  .LE*  Q.0>  VCOR  • V i 08  T8  750 
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60* 
61. 
62. 
63. 
64  • 

65. 

66. 

67. 

68. 
69. 
70- 

71. 

72. 

73. 

74. 

75. 

76. 

77. 

78. 

79. 

80. 
81. 
82. 

83. 

84. 

85. 

86. 

87. 

88. 
99. 

90. 

91. 

92. 

93. 

94. 

95. 

96. 

97. 

98. 

99. 
100. 
101. 
102. 

103. 

104. 

105. 

106. 

107. 

108. 


C if  past  tbp  is  in  wate*  e* 

tK  (DT  ,QE«  0*0)  VCOK  • 0*0  | q9  t8  750 

c 

C COMPUTE  VBUU-^E  CORRECTION  EBR  OUT  IF  WATER 

IP  ((IS  »N£.  l)  • AND . <IS  »NE*  4))  Q8  T0  350 
c CYLINDER  AND  HfcCTANGULAW  PLATE  H6D 

*0(1)  • DEPTH*  . DB/2* 0 

H0( l)  • 08 

U»  Tfl  600 

C TRIANGLE  BBD 

350  CONTINUE 

IP  (IS  «N£ . 2)  GO  TB  4Q0 

C *•**  *0(1)  SMpULO  *LSB  ChaNQe  IP  Hj  is  TB  BE  DBRRfCt  AUS0  HD 

-0(1)  ■ H*w/AbS(OT) 

37  TB  600 
400  CONTINUE 
C 

600  CONTINUE 

CALL  BB0VV8L  ( I S. ABS ( DT ) t *D( | ) , T. VCBR ) 

DT  • 0.0 
C 

c VOLUME  in  WATER 

750  CONTINUE 

VO(I)  ■ V - VCBR 

oehthhk  i » ■ oa 

OfcPTHT( I ) • DT 

c 

c WE  I SWT  BF  WATER 

wiiTw  « VD(II«*RHB0 

hdisplac  ■ wdisPLaC  ♦ WQtW 
SUHCuk  ■ SUHCbk  ♦ wqTw*XD( ( ) 

PU(I)  ■ WGTW 

C ■ 

c COMPUTE  THE  added  "I  BE  THE  WATER  B8DY 

C ABOUT  ITS  BWN  AXIS 

CALL  HBOVHl  ( IS»hd( I ). WD( U.T.WB9DVHI) 
aUwhi  ■ VO ( j ) *RWB*  ( XD ( i ) . B'J0YC3K>**2 
AUUMI  • AOQH I 4.  wqTW*wbBDYMI/q  + ADWMI 

1000  CONTINUE 
C 

C C9HPUTE  the  CENTER  Be  Bu»YANCY 

c 

HuOTLttK  • SUHCB*/WDISPLAC 
D£»THCB  ■ De^THK  . BUOTCBK 

C v 

RETURN 

c 

END 
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1. 

2. 

3* 

4. 

5. 
6 • 
7. 
8 » 
9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 
21. 
22* 
23* 

24. 

25. 

26. 

27. 

28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 

41. 

42. 

43. 

44. 

45. 

46. 

47. 

48. 

49. 

50. 
91. 

52. 

53. 

54. 

55. 

56. 

57. 

58. 

59. 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 


c 


X 

c 

c 

c 

c 

c 

c 


SUBROUT ! Np  HUOVDAMP 


VERSION  1*0  SEP/  1976  R.  GOLDSMITH 

this  routine  computes  the  buoy  damping  moment 

ASSUMPTIONS  t-  THAT  the  ANGLE  §F  ROLL  IS  'SMALL’ 

so  that  the  horizontal  Component  ®f 
BU0T  MOTION  and  water  VELOCITY  have 
PERPENDICULAR  EPPECT. 

. DAMPING  FORCE  IS  LINEAR  AND  PROPORTION 
TO  SPEED 


Common  / TP  / T I me  1 / T I ME 2. T I M£D EL / FREQ/ HAVEN 

COMMON  / BINS  / NHARTS# ISNAPfc<5r>,WlDTHI50)»wEIGMT(50)»TMICK(5D), 

♦ DENSITY(50J/DISTCG<(50)/FRaCN0RM(50J 
COMMON  / WATERDIS  / HDC3n>/HD(5D)/XP(5o>/VD(5o).FD(5o)/ 

♦ DePThH ( 50 ) i DEPtHT ( 50) 

COMMON  / CONSTANT  / Pi /RTOD/WHO,^ 

COMMON  / BUOY  / NPMAX,RrfLzAVfcRGAMplTHETABAR, PERIOD*)/ 

♦ DEHTHK#«ut)YCGK,OEPTHCGf  8U0YCBK.DEPTHCB.  WDISPLAC 
COMMON  / c^FFS  / DRAU(S) /c9EFM(5) 

COMMON  / MOMENTS  / BUM YM f . ADDMI , V I RT INRT/ WaTERI M, 8U0YMR, 

♦ BUOYMDT/BUOYMD,WATERMD/DAMPM 


ALPHAT«  4.o*RH0*THETABAR/(3.o*pn 

THIS  SECTION  COMpUTES  JUST  THE  BUOy  DAMpJNo 

BUOYMOT  ■ 0*0 
00  50  I • l.NPARTS 

CHECK  IF  ITS  but  of  WATER 
IF  (VD< I)  ,LE.  0*0)  SO  TO  50 
CHECK  SHAPE 

(10  TO  (20/30*30)/  ISMAPE(I) 

CYLINDER 

20  CONTINUE 

XB  ■ D£PTHB( I ) • OEPTHCG 
XT  ■ OePTHT<I)  « DfcPTHCG 

pmU  ■ O.P5»wIDTH( DMSIqN  (Xfl**4/XB)  • SIGN  (XT**4/XT)  ) 
30  TO  4o 

TRIANGULAR  RECTANGULAR  PLATES 

30  CONTINUE 

PAHEA  ■ VD( I )/THJCK( I) 

XC  • BUpTCGK  . XQ( I ) 

PMU  ■ PaREA*XC«XC«FRACN8RM( 1 1 
30  TO  40 
40  CONTINUE 

ALPHA  ■ ALPHAT*DRAU( ISWAPE( I ) ) 

8U0YMUT  • BU8VMUT  ♦ ALPHA*PMQ 
OUTPUT  ALP«A*PmD 
50  CONTINUE 


NOTE)  to  simplify  the  COMPUTATION  this  TERM  HAS  BEEN 
computed  as  a constant*  the  frequency 
CONTRIBUTION  IS  MULTIPLIED  in  at  the  BEGINNING 
OF  the  main  FREQUENCY  ITERATION  IN  TH£  MAIN 
PR03RAM. 


60. 
61. 
62. 
63* 
64. 
65* 
66. 
67. 
68  • 

69. 

70. 

71. 

72. 

73. 

74. 

75. 

76. 

77. 

78. 

79. 
80* 
81. 
82. 
S3. 
8*  • 
85. 
96  . 

87. 

88. 
89. 
90* 

91. 

92. 

93. 

94. 

95. 

96. 

97. 

98. 

99. 
100. 
101. 
102* 

103. 

104. 

105. 

106. 

107. 

108. 

109. 

110. 
ill* 

1 IS* 

1 1 3 • 
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C 

return 

c 

c 

ENTRY  WATERDAMP 

C 

C WATER  MOMENT  OF  DAMPING 

C 

HF.TAT  ■ 4.0*RHB*AVERGAMP*FPEW/(3»0«Pn 
waTEHMD  « 0.0 
C 

C COMPUTE  OAMRINU 

08  1000  I • l.NPAHTS 

C CWEC*  IF  ITS  BUT  BF  WATER 

IF  ( VD ( I)  »LE.  0*0)  (iO  TB  1000 
C CHECK  SHAKE 

OB  T9  ( l cn*  300' 300) * ISHAKE(I) 

C CYLINDERS 

100  CBNT1WUE 

Xb  • OEPTHBIJ)  - OEKTHCQ 
*T  • OEpTHT<I)  - OEKTHCQ 
wAvEXR  ■ wAvEN*Xa*2»0 
wavexT  ■ WAVE^*XT#2*0 

'.AV£NP2  a WAVrN»WAVEN 

Tt,HMH  . (<-wAVEXB  - 1 *0) #EXp  { "i* AVEXB ) * 1*0). 

TfcKMT  a ((-WAVEXT  . 1*0)#EXP  (-WAVEXT)  ♦ 1*0) 

wMU  a <te«*9  - TER'’,T>*WIt)TH(I  )*EXP  ( -2.0*waVeN*DepTHCQ) 

♦ / ( 4 . 0*wA VtNp2 ) 

OB  T9  900 

C 

C TRIANGULAR  AND  RECTANGULAR  PLATES  , 

300  CBNTINuc 

RAKE A a VO ( I ) /THICK ( I ) 

XC  ■ BUflVCUK  . XD ( i ) 
lQ  • OEpTh<  - xDU) 

C WATER  MHKENT 

WHO  a PAREA»XC.FHACNBRH(  I ),»EXP  (*WAVER#ZC) 

38  TO  900 
900  CBNTINUE 
c 

C GET  total  Dahring 

BETA  a BETAT.ORAGI ISHAPEl I) ) 
who  ■ BfTA#«MD 
wATERmD".  wATERmD  * WMO 
1000  CONTINUE 
c 

C TOTAL  DAMRING  MOMENT 

c 

OAMKM  a BU9TM0  ♦ WATERED 

c 

RETyRN 

C 

END 
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C 

c 

c 

c 

c 

c 

c 


c 

c 


SUBROUTINE  watrinkt 

VERSION  1*0 


SEP*  1976 


R.  GOLDSMITH 


THIS  R4UTINE  COMPUTES  TH£  COEFFICIENT  FOR  THE  INERTIA 
MOMENT  DUE  TO  WATER  PARTICLE  ACCELERATION, 


COMMON 

COMMON 


TP  / TTHEi»TIME?>TI*EDELiFREQ,wAVEN 

BINS  / N-ARTS, ISHAHE(50),w!DTH(50) ,HEIGHT(50) ,THICK(50)» 
DENSITY! 50).DISTCgk(50)#eRACN0Rm 150) 

COMMON  / WATEHOIS  / wD(5o).w9(5o).XD(50).VD(5o).FD(50). 

DfcHTwB(sO),DEPTHTC50> 

CONSTANT  / Hl.HTBO.HHO.Q 

HUHY  / NPHaX,K*L»AVERGAHP,THETABaR.P£RIODo, 

OEPTh<#BUBTCGK| depthcg, buoycok*  depthcb. wdjsplaC 

CPfFs  DHAQI5) iCOEFH(S) 

MPPENTb  / UUOv^I/AIJOMI.vlRTlNRT.wATERIMfBUOYMR,  . 
BUOTHQT  j buoyhd * watermd.damph 


common 

COMMON 

COMMON 

COMMON 


c 

C 

c 

c 

c 

c 


OETAT  . p I *RH9/4 • 0 

CONSI  a BrTAT/(WAVEN»WAVEN)»EXP 

wATEHIm  • 0*0 


(.wAVEN*OEPTHCQj 


COMPUTE  WATER  PARTICLE 
DO  1000  I ■ l.NPARTS 

CHECK  IF  ITS  OUT  OF 
IF  <VO(I)  «lE.  0*0)  GO  TO  1000 
CHECK  SHAPE 

00  TO  ( 100'300. 300>*  ISHAPE(I) 


inertia  moment 
water 


CYLINDERS 

ioo  continue 

XO  ■ DE»ThB( I ) - DEPTMCQ 

xr  ■ OE°THT ( I ) - depthco 

wAVENXH  ■ wAVEN**a 

wAYENXT  . WAVEN*XT 

TtRM  J H a (.WAVENXB  • 1.0>  *£XP 

TfeRMlB  • * TEOMlB  ♦ 1*0 

TERMlT  , (-WAVENXT  - 1«0)*EXP 

TfeRMIT  a ♦ TEPM1T  ♦ 1.0 

*Ml  ■ c0^^1'#<wl'3THU  )*wIDTH(  I ) 

00  TB  900 


•WAVENXB) 

•WAVENXT, 

)*(TERM1B 


terhiti 


300  CONTINUE 
••••CURHLNTLV  set 
wM I a 0.0 

30  TO  900 


TRIANGULAR  AND  RECTANGULAR  PLATES 
TB  ZERO 


900  CONTINUE 

WM | a •WMI*C0EFM<  J SHAPE I I ) ) 
WATtRJM  a WATER! M ♦ WM  I 

1000  CONTINUE 
•RETURN 
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